mirror/kissfft
1
0
Fork 0
mirror of https://github.com/mborgerding/kissfft.git synced 2025-06-03 09:08:10 -04:00
Code Issues Packages Projects Releases Wiki Activity

going to bed

Browse Source
This commit is contained in:
Mark Borgerding 2003-11-07 03:42:14 +00:00
parent ae305ca400
commit e98f9ff29a
1 changed files with 7 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
kiss_fft.c
View File

@ -292,7 +292,7 @@ void bfly_generic(
kiss_fft_cpx * scratch2 = scratch + p; kiss_fft_cpx * scratch2 = scratch + p;
kiss_fft_cpx * tw = st->twiddles; kiss_fft_cpx * tw = st->twiddles;
kiss_fft_cpx tlo,t3,t4; kiss_fft_cpx tlo,t3,t4;
kiss_fft_cpx *Foutlo,*Fouthi,*tw2; kiss_fft_cpx *Foutlo,*Fouthi;
fsm = fstride*m; fsm = fstride*m;
halfp=p/2; halfp=p/2;
@ -315,7 +315,6 @@ void bfly_generic(
scratch2 = scratch + p; scratch2 = scratch + p;
uf=u*fstride; uf=u*fstride;
/* d==0 */
Foutlo=Fout; Foutlo=Fout;
for ( q=1 ; q<p ; ++q ) { for ( q=1 ; q<p ; ++q ) {
Foutlo += m; Foutlo += m;
@ -325,14 +324,12 @@ void bfly_generic(
for ( q=1; q<p ; ++q ) { for ( q=1; q<p ; ++q ) {
C_ADDTO(*Fout , scratch[q] ); C_ADDTO(*Fout , scratch[q] );
Foutlo=Fout; Foutlo=Fout+m;
Fouthi=Fout + mp; Fouthi=Fout+mp;
for ( d=1; d<=halfp;++d) { do{
tlo = *scratch2++; tlo = *scratch2++;
Foutlo += m;
t3.r = scratch[q].r * tlo.r; t3.r = scratch[q].r * tlo.r;
t3.i = scratch[q].r * tlo.i; t3.i = scratch[q].r * tlo.i;
t4.r = scratch[q].i * tlo.i; t4.r = scratch[q].i * tlo.i;
@ -341,10 +338,10 @@ void bfly_generic(
Fouthi -= m; Fouthi -= m;
Foutlo->r += t3.r - t4.r; Foutlo->r += t3.r - t4.r;
Fouthi->r += t3.r + t4.r;
Foutlo->i += t3.i - t4.i; Foutlo->i += t3.i - t4.i;
Fouthi->r += t3.r + t4.r;
Fouthi->i -= t3.i + t4.i; Fouthi->i -= t3.i + t4.i;
} }while( (Foutlo += m) < Fouthi );
} }
++Fout; ++Fout;
} }
@ -372,7 +369,7 @@ void fft_work(
switch (p) { switch (p) {
case 2: bfly2(Fout,fstride,st,m); break; case 2: bfly2(Fout,fstride,st,m); break;
#if 0 #if 1
case 3: bfly3(Fout,fstride,st,m); break; case 3: bfly3(Fout,fstride,st,m); break;
#endif #endif
case 4: bfly4(Fout,fstride,st,m); break; case 4: bfly4(Fout,fstride,st,m); break;