From b92dce32703aa0a7899db22c4d28e8cab6b5fdcf Mon Sep 17 00:00:00 2001
From: Gregory Maxwell <greg@xiph.org>
Date: Sun, 28 Jun 2009 19:51:30 -0400
Subject: [PATCH] Additional mathops.h tests.
---
libcelt/mathops.h | 4 +--
tests/mathops-test.c | 65 ++++++++++++++++++++++++++++++++++++++++++++
2 files changed, 67 insertions(+), 2 deletions(-)
diff --git a/libcelt/mathops.h b/libcelt/mathops.h
index 62e8e4aa6..72a4fbf4f 100644
--- a/libcelt/mathops.h
+++ b/libcelt/mathops.h
@@ -155,8 +155,8 @@ static inline float celt_exp2(float x)
}
#else
-#define celt_log2(x) (1.442695*log(x))
-#define celt_exp2(x) (exp(0.69315*(x)))
+#define celt_log2(x) (1.442695040888963387*log(x))
+#define celt_exp2(x) (exp(0.6931471805599453094*(x)))
#endif
#endif
diff --git a/tests/mathops-test.c b/tests/mathops-test.c
index 214e2c212..cb3c1325a 100644
--- a/tests/mathops-test.c
+++ b/tests/mathops-test.c
@@ -71,10 +71,75 @@ void testrsqrt(void)
}
}
+#ifndef FIXED_POINT
+void testlog2(void)
+{
+ float x;
+ for (x=0.001;x<1677700.0;x+=(x/8.0))
+ {
+ float error = fabs((1.442695040888963387*log(x))-celt_log2(x));
+ if (error>0.001)
+ {
+ fprintf (stderr, "celt_log2 failed: fabs((1.442695040888963387*log(x))-celt_log2(x))>0.001 (x = %f, error = %f)\n", x,error);
+ ret = 1;
+ }
+ }
+}
+
+void testexp2(void)
+{
+ float x;
+ for (x=-11.0;x<24.0;x+=0.0007)
+ {
+ float error = fabs(x-(1.442695040888963387*log(celt_exp2(x))));
+ if (error>0.0005)
+ {
+ fprintf (stderr, "celt_exp2 failed: fabs(x-(1.442695040888963387*log(celt_exp2(x))))>0.0005 (x = %f, error = %f)\n", x,error);
+ ret = 1;
+ }
+ }
+}
+
+void testexp2log2(void)
+{
+ float x;
+ for (x=-11.0;x<24.0;x+=0.0007)
+ {
+ float error = fabs(x-(celt_log2(celt_exp2(x))));
+ if (error>0.001)
+ {
+ fprintf (stderr, "celt_log2/celt_exp2 failed: fabs(x-(celt_log2(celt_exp2(x))))>0.001 (x = %f, error = %f)\n", x,error);
+ ret = 1;
+ }
+ }
+}
+#else
+void testilog2(void)
+{
+ celt_word32_t x;
+ for (x=1;x<=268435455;x+=127)
+ {
+ celt_word32_t error = abs(celt_ilog2(x)-(int)floor(log2(x)));
+ if (error!=0)
+ {
+ printf("celt_ilog2 failed: celt_ilog2(x)!=floor(log2(x)) (x = %d, error = %d)\n",x,error);
+ ret = 1;
+ }
+ }
+}
+#endif
+
int main(void)
{
testdiv();
testsqrt();
testrsqrt();
+#ifndef FIXED_POINT
+ testlog2();
+ testexp2();
+ testexp2log2();
+#else
+ testilog2();
+#endif
return ret;
}
--
GitLab