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

*** empty log message ***

Browse Source
This commit is contained in:
Mark Borgerding 2003-10-18 00:32:54 +00:00
parent 6f8bcedc24
commit a3d3217ae6
2 changed files with 7 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
.cvsignore Normal file
View File

@ -0,0 +1,2 @@
*.pyc
*.o

11
kiss_fft.c
View File

@ -180,16 +180,14 @@ void bfly3(
int u; int u;
kiss_fft_cpx * scratch = st->scratch; kiss_fft_cpx * scratch = st->scratch;
kiss_fft_cpx * twiddles = st->twiddles; kiss_fft_cpx * twiddles = st->twiddles;
kiss_fft_cpx t[6]; kiss_fft_cpx t[2];
kiss_fft_cpx epi3; kiss_fft_cpx epi3;
//kiss_fft_cpx epi3 = { -0.5 , -0.866025403784439 }; epi3 = twiddles[fstride*m];
epi3 = twiddles[fstride*m];
Fout0=Fout; Fout0=Fout;
Fout1=Fout0+m; Fout1=Fout0+m;
Fout2=Fout0+2*m; Fout2=Fout0+2*m;
//printf("epi3=%e,%ei\n",epi3.r,epi3.i);
for ( u=0; u<m; ++u ) { for ( u=0; u<m; ++u ) {
kiss_fft_cpx sum23,t0pt1,t0mt1; kiss_fft_cpx sum23,t0pt1,t0mt1;
@ -198,7 +196,6 @@ void bfly3(
Fout1->r /= 3; Fout1->i /= 3; Fout1->r /= 3; Fout1->i /= 3;
Fout2->r /= 3; Fout2->i /= 3; Fout2->r /= 3; Fout2->i /= 3;
#endif #endif
scratch[0] = *Fout0; scratch[0] = *Fout0;
C_MUL(t[0],*Fout1 , twiddles[fstride*u ] ); C_MUL(t[0],*Fout1 , twiddles[fstride*u ] );
@ -214,8 +211,8 @@ void bfly3(
t0mt1.r *= epi3.i; t0mt1.r *= epi3.i;
t0mt1.i *= epi3.i; t0mt1.i *= epi3.i;
sum23.r = t0pt1.r - t0mt1.i ; sum23.r = t0pt1.r - t0mt1.i;
sum23.i = t0pt1.i + t0mt1.r ; sum23.i = t0pt1.i + t0mt1.r;
C_ADD( *Fout1, scratch[0] , sum23 ); C_ADD( *Fout1, scratch[0] , sum23 );