Test cleanups and renaming.
[opus.git] / celt / tests / test_unit_mathops.c
1 #ifdef HAVE_CONFIG_H
2 #include "config.h"
3 #endif
4
5 #ifndef CUSTOM_MODES
6 #define CUSTOM_MODES
7 #endif
8
9 #define CELT_C
10
11 #include "mathops.c"
12 #include "entenc.c"
13 #include "entdec.c"
14 #include "entcode.c"
15 #include "bands.c"
16 #include "vq.c"
17 #include "cwrs.c"
18 #include <stdio.h>
19 #include <math.h>
20
21 #ifdef FIXED_POINT
22 #define WORD "%d"
23 #else
24 #define WORD "%f"
25 #endif
26
27 #ifdef FIXED_DEBUG
28 long long celt_mips=0;
29 #endif
30 int ret = 0;
31
32 void testdiv(void)
33 {
34    opus_int32 i;
35    for (i=1;i<=327670;i++)
36    {
37       double prod;
38       opus_val32 val;
39       val = celt_rcp(i);
40 #ifdef FIXED_POINT
41       prod = (1./32768./65526.)*val*i;
42 #else
43       prod = val*i;
44 #endif
45       if (fabs(prod-1) > .00025)
46       {
47          fprintf (stderr, "div failed: 1/%d="WORD" (product = %f)\n", i, val, prod);
48          ret = 1;
49       }
50    }
51 }
52
53 void testsqrt(void)
54 {
55    opus_int32 i;
56    for (i=1;i<=1000000000;i++)
57    {
58       double ratio;
59       opus_val16 val;
60       val = celt_sqrt(i);
61       ratio = val/sqrt(i);
62       if (fabs(ratio - 1) > .0005 && fabs(val-sqrt(i)) > 2)
63       {
64          fprintf (stderr, "sqrt failed: sqrt(%d)="WORD" (ratio = %f)\n", i, val, ratio);
65          ret = 1;
66       }
67       i+= i>>10;
68    }
69 }
70
71 void testbitexactcos(void)
72 {
73    int i;
74    opus_int32 min_d,max_d,last,chk;
75    chk=max_d=0;
76    last=min_d=32767;
77    for(i=64;i<=16320;i++)
78    {
79       opus_int32 d;
80       opus_int32 q=bitexact_cos(i);
81       chk ^= q*i;
82       d = last - q;
83       if (d>max_d)max_d=d;
84       if (d<min_d)min_d=d;
85       last = q;
86    }
87    if ((chk!=89408644)||(max_d!=5)||(min_d!=0)||(bitexact_cos(64)!=32767)||
88        (bitexact_cos(16320)!=200)||(bitexact_cos(8192)!=23171))
89    {
90       fprintf (stderr, "bitexact_cos failed\n");
91       ret = 1;
92    }
93 }
94
95 void testbitexactlog2tan(void)
96 {
97    int i,fail;
98    opus_int32 min_d,max_d,last,chk;
99    fail=chk=max_d=0;
100    last=min_d=15059;
101    for(i=64;i<8193;i++)
102    {
103       opus_int32 d;
104       opus_int32 mid=bitexact_cos(i);
105       opus_int32 side=bitexact_cos(16384-i);
106       opus_int32 q=bitexact_log2tan(mid,side);
107       chk ^= q*i;
108       d = last - q;
109       if (q!=-1*bitexact_log2tan(side,mid))
110         fail = 1;
111       if (d>max_d)max_d=d;
112       if (d<min_d)min_d=d;
113       last = q;
114    }
115    if ((chk!=15821257)||(max_d!=61)||(min_d!=-2)||fail||
116        (bitexact_log2tan(32767,200)!=15059)||(bitexact_log2tan(30274,12540)!=2611)||
117        (bitexact_log2tan(23171,23171)!=0))
118    {
119       fprintf (stderr, "bitexact_log2tan failed\n");
120       ret = 1;
121    }
122 }
123
124 #ifndef FIXED_POINT
125 void testlog2(void)
126 {
127    float x;
128    for (x=0.001;x<1677700.0;x+=(x/8.0))
129    {
130       float error = fabs((1.442695040888963387*log(x))-celt_log2(x));
131       if (error>0.0009)
132       {
133          fprintf (stderr, "celt_log2 failed: fabs((1.442695040888963387*log(x))-celt_log2(x))>0.001 (x = %f, error = %f)\n", x,error);
134          ret = 1;
135       }
136    }
137 }
138
139 void testexp2(void)
140 {
141    float x;
142    for (x=-11.0;x<24.0;x+=0.0007)
143    {
144       float error = fabs(x-(1.442695040888963387*log(celt_exp2(x))));
145       if (error>0.0002)
146       {
147          fprintf (stderr, "celt_exp2 failed: fabs(x-(1.442695040888963387*log(celt_exp2(x))))>0.0005 (x = %f, error = %f)\n", x,error);
148          ret = 1;
149       }
150    }
151 }
152
153 void testexp2log2(void)
154 {
155    float x;
156    for (x=-11.0;x<24.0;x+=0.0007)
157    {
158       float error = fabs(x-(celt_log2(celt_exp2(x))));
159       if (error>0.001)
160       {
161          fprintf (stderr, "celt_log2/celt_exp2 failed: fabs(x-(celt_log2(celt_exp2(x))))>0.001 (x = %f, error = %f)\n", x,error);
162          ret = 1;
163       }
164    }
165 }
166 #else
167 void testlog2(void)
168 {
169    opus_val32 x;
170    for (x=8;x<1073741824;x+=(x>>3))
171    {
172       float error = fabs((1.442695040888963387*log(x/16384.0))-celt_log2(x)/1024.0);
173       if (error>0.003)
174       {
175          fprintf (stderr, "celt_log2 failed: x = %ld, error = %f\n", (long)x,error);
176          ret = 1;
177       }
178    }
179 }
180
181 void testexp2(void)
182 {
183    opus_val16 x;
184    for (x=-32768;x<15360;x++)
185    {
186       float error1 = fabs(x/1024.0-(1.442695040888963387*log(celt_exp2(x)/65536.0)));
187       float error2 = fabs(exp(0.6931471805599453094*x/1024.0)-celt_exp2(x)/65536.0);
188       if (error1>0.0002&&error2>0.00004)
189       {
190          fprintf (stderr, "celt_exp2 failed: x = "WORD", error1 = %f, error2 = %f\n", x,error1,error2);
191          ret = 1;
192       }
193    }
194 }
195
196 void testexp2log2(void)
197 {
198    opus_val32 x;
199    for (x=8;x<65536;x+=(x>>3))
200    {
201       float error = fabs(x-0.25*celt_exp2(celt_log2(x)))/16384;
202       if (error>0.004)
203       {
204          fprintf (stderr, "celt_log2/celt_exp2 failed: fabs(x-(celt_exp2(celt_log2(x))))>0.001 (x = %ld, error = %f)\n", (long)x,error);
205          ret = 1;
206       }
207    }
208 }
209
210 void testilog2(void)
211 {
212    opus_val32 x;
213    for (x=1;x<=268435455;x+=127)
214    {
215       opus_val32 error = abs(celt_ilog2(x)-(int)floor(log2(x)));
216       if (error!=0)
217       {
218          printf("celt_ilog2 failed: celt_ilog2(x)!=floor(log2(x)) (x = %d, error = %d)\n",x,error);
219          ret = 1;
220       }
221    }
222 }
223 #endif
224
225 int main(void)
226 {
227    testbitexactcos();
228    testbitexactlog2tan();
229    testdiv();
230    testsqrt();
231    testlog2();
232    testexp2();
233    testexp2log2();
234 #ifdef FIXED_POINT
235    testilog2();
236 #endif
237    return ret;
238 }