Additional mathops.h tests.

1 files changed, 65 insertions, 0 deletions

diff --git a/tests/mathops-test.c b/tests/mathops-test.c
index 214e2c2..cb3c132 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;
 }