Chris@366: #include <stdio.h>
Chris@366: #include <stdlib.h>
Chris@366: #include <sys/times.h>
Chris@366: #include <unistd.h>
Chris@366: #include "kiss_fft.h"
Chris@366: #include "kiss_fftr.h"
Chris@366: #include "kiss_fftnd.h"
Chris@366: #include "kiss_fftndr.h"
Chris@366: 
Chris@366: #include "pstats.h"
Chris@366: 
Chris@366: static
Chris@366: int getdims(int * dims, char * arg)
Chris@366: {
Chris@366:     char *s;
Chris@366:     int ndims=0;
Chris@366:     while ( (s=strtok( arg , ",") ) ) {
Chris@366:         dims[ndims++] = atoi(s);
Chris@366:         //printf("%s=%d\n",s,dims[ndims-1]);
Chris@366:         arg=NULL;
Chris@366:     }
Chris@366:     return ndims;
Chris@366: }
Chris@366: 
Chris@366: int main(int argc,char ** argv)
Chris@366: {
Chris@366:     int k;
Chris@366:     int nfft[32];
Chris@366:     int ndims = 1;
Chris@366:     int isinverse=0;
Chris@366:     int numffts=1000,i;
Chris@366:     kiss_fft_cpx * buf;
Chris@366:     kiss_fft_cpx * bufout;
Chris@366:     int real = 0;
Chris@366: 
Chris@366:     nfft[0] = 1024;// default
Chris@366: 
Chris@366:     while (1) {
Chris@366:         int c = getopt (argc, argv, "n:ix:r");
Chris@366:         if (c == -1)
Chris@366:             break;
Chris@366:         switch (c) {
Chris@366:             case 'r':
Chris@366:                 real = 1;
Chris@366:                 break;
Chris@366:             case 'n':
Chris@366:                 ndims = getdims(nfft, optarg );
Chris@366:                 if (nfft[0] != kiss_fft_next_fast_size(nfft[0]) ) {
Chris@366:                     int ng = kiss_fft_next_fast_size(nfft[0]);
Chris@366:                     fprintf(stderr,"warning: %d might be a better choice for speed than %d\n",ng,nfft[0]);
Chris@366:                 }
Chris@366:                 break;
Chris@366:             case 'x':
Chris@366:                 numffts = atoi (optarg);
Chris@366:                 break;
Chris@366:             case 'i':
Chris@366:                 isinverse = 1;
Chris@366:                 break;
Chris@366:         }
Chris@366:     }
Chris@366:     int nbytes = sizeof(kiss_fft_cpx);
Chris@366:     for (k=0;k<ndims;++k)
Chris@366:         nbytes *= nfft[k];
Chris@366: 
Chris@366: #ifdef USE_SIMD        
Chris@366:     numffts /= 4;
Chris@366:     fprintf(stderr,"since SIMD implementation does 4 ffts at a time, numffts is being reduced to %d\n",numffts);
Chris@366: #endif
Chris@366: 
Chris@366:     buf=(kiss_fft_cpx*)KISS_FFT_MALLOC(nbytes);
Chris@366:     bufout=(kiss_fft_cpx*)KISS_FFT_MALLOC(nbytes);
Chris@366:     memset(buf,0,nbytes);
Chris@366: 
Chris@366:     pstats_init();
Chris@366: 
Chris@366:     if (ndims==1) {
Chris@366:         if (real) {
Chris@366:             kiss_fftr_cfg st = kiss_fftr_alloc( nfft[0] ,isinverse ,0,0);
Chris@366:             if (isinverse)
Chris@366:                 for (i=0;i<numffts;++i)
Chris@366:                     kiss_fftri( st ,(kiss_fft_cpx*)buf,(kiss_fft_scalar*)bufout );
Chris@366:             else
Chris@366:                 for (i=0;i<numffts;++i)
Chris@366:                     kiss_fftr( st ,(kiss_fft_scalar*)buf,(kiss_fft_cpx*)bufout );
Chris@366:             free(st);
Chris@366:         }else{
Chris@366:             kiss_fft_cfg st = kiss_fft_alloc( nfft[0] ,isinverse ,0,0);
Chris@366:             for (i=0;i<numffts;++i)
Chris@366:                 kiss_fft( st ,buf,bufout );
Chris@366:             free(st);
Chris@366:         }
Chris@366:     }else{
Chris@366:         if (real) {
Chris@366:             kiss_fftndr_cfg st = kiss_fftndr_alloc( nfft,ndims ,isinverse ,0,0);
Chris@366:             if (isinverse)
Chris@366:                 for (i=0;i<numffts;++i)
Chris@366:                     kiss_fftndri( st ,(kiss_fft_cpx*)buf,(kiss_fft_scalar*)bufout );
Chris@366:             else
Chris@366:                 for (i=0;i<numffts;++i)
Chris@366:                     kiss_fftndr( st ,(kiss_fft_scalar*)buf,(kiss_fft_cpx*)bufout );
Chris@366:             free(st);
Chris@366:         }else{
Chris@366:             kiss_fftnd_cfg st= kiss_fftnd_alloc(nfft,ndims,isinverse ,0,0);
Chris@366:             for (i=0;i<numffts;++i)
Chris@366:                 kiss_fftnd( st ,buf,bufout );
Chris@366:             free(st);
Chris@366:         }
Chris@366:     }
Chris@366: 
Chris@366:     free(buf); free(bufout);
Chris@366: 
Chris@366:     fprintf(stderr,"KISS\tnfft=");
Chris@366:     for (k=0;k<ndims;++k)
Chris@366:         fprintf(stderr, "%d,",nfft[k]);
Chris@366:     fprintf(stderr,"\tnumffts=%d\n" ,numffts);
Chris@366:     pstats_report();
Chris@366: 
Chris@366:     kiss_fft_cleanup();
Chris@366: 
Chris@366:     return 0;
Chris@366: }
Chris@366: