mirror/kissfft
1
0
Fork 0
mirror of https://github.com/mborgerding/kissfft.git synced 2025-05-27 21:20:27 -04:00
Code Issues Packages Projects Releases Wiki Activity

*** empty log message ***

Browse Source
This commit is contained in:
Mark Borgerding 2003-08-09 13:12:56 +00:00
parent ccbc48dc63
commit db556661ed
2 changed files with 8 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Makefile
View File

@ -1,5 +1,5 @@
all: kiss_fft_s kiss_fft_f kiss_fft_d freqpeak tones
all: kiss_fft_s kiss_fft_f kiss_fft_d freqpeak tones testsig
kiss_fft_s: kiss_fft.h kiss_fft.c
gcc -Wall -O3 -o kiss_fft_s -DFIXED_POINT -DFFT_UTIL kiss_fft.c -lm

16
freqpeak.c
View File

@ -28,18 +28,14 @@ int main(int argc, char ** argv)
memset( rmag2 , 0 , sizeof(rmag2) );
while ( fread( sampsin , sizeof(short) * 2*NFFT, 1 , stdin ) == 1 ) {
//perform two ffts in parallel by packing the channels into the real and imaginary
//
//perform two ffts in parallel by packing the left&right channels into the real and imaginary
for (k=0;k<NFFT;++k) {
//cbuf[k].r = 0;
//cbuf[k].i = 0;
cbuf[k].r = sampsin[2*k];
cbuf[k].i = sampsin[2*k+1];
//cbuf[k].i = sampsin[2*k+1];
}
if (removedc){
if (removedc) {
float dcr=0,dci=0;
for (k=0;k<NFFT;++k){
dcr += cbuf[k].r;
@ -48,14 +44,16 @@ int main(int argc, char ** argv)
dcr /= NFFT;
dci /= NFFT;
for (k=0;k<NFFT;++k){
for (k=0;k<NFFT;++k) {
cbuf[k].r -= dcr;
cbuf[k].i -= dci;
}
}
// perform the fft on the L+R packed buffer
kiss_fft( st , cbuf );
// get the half-symmetric FFTs for the Left and Right Channels
for (k=0;k<NFFT/2;++k) {
int k2 = (NFFT-k)%NFFT;
kiss_fft_cpx r,l;