54 void check(kiss_fft_scalar  * in,kiss_fft_scalar  * out,int nfft,int isinverse)  in check()  argument
59     for (bin=0;bin<nfft/2;++bin) {  in check()
63         for (k=0;k<nfft;++k) {  in check()
64            double phase = 2*M_PI*(k+.5+.25*nfft)*(bin+.5)/nfft;  in check()
67            re /= nfft/4;  in check()
77     printf("nfft=%d inverse=%d,snr = %f\n",nfft,isinverse,snr );  in check()
84 void check_inv(kiss_fft_scalar  * in,kiss_fft_scalar  * out,int nfft,int isinverse)  in check_inv()  argument
89    for (bin=0;bin<nfft;++bin) {  in check_inv()
93       for (k=0;k<nfft/2;++k) {  in check_inv()
94          double phase = 2*M_PI*(bin+.5+.25*nfft)*(k+.5)/nfft;  in check_inv()
107    printf("nfft=%d inverse=%d,snr = %f\n",nfft,isinverse,snr );  in check_inv()
115 void test1d(int nfft,int isinverse)  in test1d()  argument
118     size_t buflen = sizeof(kiss_fft_scalar)*nfft;  in test1d()
123     opus_val16  * window= (opus_val16*)malloc(sizeof(opus_val16)*nfft/2);  in test1d()
126     clt_mdct_init(&cfg, nfft, 0);  in test1d()
127     for (k=0;k<nfft;++k) {  in test1d()
131     for (k=0;k<nfft/2;++k) {  in test1d()
134     for (k=0;k<nfft;++k) {  in test1d()
140        for (k=0;k<nfft;++k) {  in test1d()
141           in[k] /= nfft;  in test1d()
145     for (k=0;k<nfft;++k)  in test1d()
151        for (k=0;k<nfft;++k)  in test1d()
153        clt_mdct_backward(&cfg,in,out, window, nfft/2, 0, 1);  in test1d()
155        for (k=0;k<nfft/4;++k)  in test1d()
156           out[nfft-k-1] = out[nfft/2+k];  in test1d()
157        check_inv(in,out,nfft,isinverse);  in test1d()
159        clt_mdct_forward(&cfg,in,out,window, nfft/2, 0, 1);  in test1d()
160        check(in_copy,out,nfft,isinverse);  in test1d()