mirror/kissfft
1
0
Fork 0
mirror of https://github.com/mborgerding/kissfft.git synced 2025-06-04 01:28:23 -04:00
Code Issues Packages Projects Releases Wiki Activity

fixed warnings

Browse Source
This commit is contained in:
Mark Borgerding 2005-06-24 01:27:31 +00:00
parent 2f70f1172d
commit 24be1c5850
1 changed files with 10 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
test/twotonetest.c
View File

@ -20,19 +20,17 @@ double two_tone_test( int nfft, int bin1,int bin2)
double f2 = bin2*2*M_PI/nfft; double f2 = bin2*2*M_PI/nfft;
double sigpow=0; double sigpow=0;
double noisepow=0; double noisepow=0;
#if FIXED_POINT==32
long maxrange = LONG_MAX;
#else
long maxrange = SHRT_MAX;/* works fine for float too*/
#endif
cfg = kiss_fftr_alloc(nfft , 0, NULL, NULL); cfg = kiss_fftr_alloc(nfft , 0, NULL, NULL);
tbuf = malloc(nfft * sizeof(kiss_fft_scalar)); tbuf = malloc(nfft * sizeof(kiss_fft_scalar));
kout = malloc(nfft * sizeof(kiss_fft_cpx)); kout = malloc(nfft * sizeof(kiss_fft_cpx));
//assert(bin1<=nfft/2);
//assert(bin2<=nfft/2);
#if FIXED_POINT==32 /* generate a signal with two tones*/
long maxrange = LONG_MAX;
#else
long maxrange = SHRT_MAX;// works fine for float too
#endif
// generate a signal with two tones
for (i = 0; i < nfft; i++) { for (i = 0; i < nfft; i++) {
tbuf[i] = (maxrange>>1)*cos(f1*i) tbuf[i] = (maxrange>>1)*cos(f1*i)
+ (maxrange>>1)*cos(f2*i); + (maxrange>>1)*cos(f2*i);
@ -45,19 +43,19 @@ double two_tone_test( int nfft, int bin1,int bin2)
double tmpi = (double)kout[i].i / (double)maxrange; double tmpi = (double)kout[i].i / (double)maxrange;
double mag2 = tmpr*tmpr + tmpi*tmpi; double mag2 = tmpr*tmpr + tmpi*tmpi;
if (i!=0 && i!= nfft/2) if (i!=0 && i!= nfft/2)
mag2 *= 2; // all bins except DC and Nyquist have symmetric counterparts implied mag2 *= 2; /* all bins except DC and Nyquist have symmetric counterparts implied*/
// if there is power in one of the expected bins, it is signal, otherwise noise /* if there is power in one of the expected bins, it is signal, otherwise noise*/
if ( i!=bin1 && i != bin2 ) if ( i!=bin1 && i != bin2 )
noisepow += mag2; noisepow += mag2;
else else
sigpow += mag2; sigpow += mag2;
} }
//printf("TEST %d,%d,%d noise @ %fdB\n",nfft,bin1,bin2,10*log10(noisepow/sigpow +1e-30) ); /*printf("TEST %d,%d,%d noise @ %fdB\n",nfft,bin1,bin2,10*log10(noisepow/sigpow +1e-30) );*/
return 10*log10(sigpow/(noisepow+1e-50) ); return 10*log10(sigpow/(noisepow+1e-50) );
} }
int main(int argc,char ** argv) int main(void)
{ {
int nfft = 4*2*2*3*5; int nfft = 4*2*2*3*5;
int i,j; int i,j;