c@409: #include "kiss_fft.h"
c@409: 
c@409: 
c@409: void check(kiss_fft_cpx  * in,kiss_fft_cpx  * out,int nfft,int isinverse)
c@409: {
c@409:     int bin,k;
c@409:     double errpow=0,sigpow=0;
c@409:     
c@409:     for (bin=0;bin<nfft;++bin) {
c@409:         double ansr = 0;
c@409:         double ansi = 0;
c@409:         double difr;
c@409:         double difi;
c@409: 
c@409:         for (k=0;k<nfft;++k) {
c@409:             double phase = -2*M_PI*bin*k/nfft;
c@409:             double re = cos(phase);
c@409:             double im = sin(phase);
c@409:             if (isinverse)
c@409:                 im = -im;
c@409: 
c@409: #ifdef FIXED_POINT
c@409:             re /= nfft;
c@409:             im /= nfft;
c@409: #endif            
c@409: 
c@409:             ansr += in[k].r * re - in[k].i * im;
c@409:             ansi += in[k].r * im + in[k].i * re;
c@409:         }
c@409:         difr = ansr - out[bin].r;
c@409:         difi = ansi - out[bin].i;
c@409:         errpow += difr*difr + difi*difi;
c@409:         sigpow += ansr*ansr+ansi*ansi;
c@409:     }
c@409:     printf("nfft=%d inverse=%d,snr = %f\n",nfft,isinverse,10*log10(sigpow/errpow) );
c@409: }
c@409: 
c@409: void test1d(int nfft,int isinverse)
c@409: {
c@409:     size_t buflen = sizeof(kiss_fft_cpx)*nfft;
c@409: 
c@409:     kiss_fft_cpx  * in = (kiss_fft_cpx*)malloc(buflen);
c@409:     kiss_fft_cpx  * out= (kiss_fft_cpx*)malloc(buflen);
c@409:     kiss_fft_cfg  cfg = kiss_fft_alloc(nfft,isinverse,0,0);
c@409:     int k;
c@409: 
c@409:     for (k=0;k<nfft;++k) {
c@409:         in[k].r = (rand() % 65536) - 32768;
c@409:         in[k].i = (rand() % 65536) - 32768;
c@409:     }
c@409: 
c@409:     kiss_fft(cfg,in,out);
c@409: 
c@409:     check(in,out,nfft,isinverse);
c@409: 
c@409:     free(in);
c@409:     free(out);
c@409:     free(cfg);
c@409: }
c@409: 
c@409: int main(int argc,char ** argv)
c@409: {
c@409:     if (argc>1) {
c@409:         int k;
c@409:         for (k=1;k<argc;++k) {
c@409:             test1d(atoi(argv[k]),0);
c@409:             test1d(atoi(argv[k]),1);
c@409:         }
c@409:     }else{
c@409:         test1d(32,0);
c@409:         test1d(32,1);
c@409:     }
c@409:     return 0;
c@409: }