src/fft.c: fix norm computation when HAVE_COMPLEX_H is undefined
authorPaul Brossier <piem@piem.org>
Sat, 3 Nov 2007 18:03:58 +0000 (19:03 +0100)
committerPaul Brossier <piem@piem.org>
Sat, 3 Nov 2007 18:03:58 +0000 (19:03 +0100)
src/fft.c

index 5cbc47807a34f5a872111a900c7d08f95e5f2cd3..e5c1a03770241f25c6de3b8a2eab5d963f4983f2 100644 (file)
--- a/src/fft.c
+++ b/src/fft.c
@@ -120,9 +120,9 @@ void aubio_fft_getspectrum(fft_data_t * spectrum, smpl_t *norm, smpl_t * phas, u
 
 void aubio_fft_getnorm(smpl_t * norm, fft_data_t * spectrum, uint_t size) {
   uint_t i;
-  norm[0] = SQR(spectrum[0]);
-  for (i=1;i<size/2;i++) norm[i] = (SQR(spectrum[i]) + SQR(spectrum[size-i]));
-  norm[size/2] = SQR(spectrum[size/2]);
+  norm[0] = spectrum[0];
+  for (i=1;i<size/2;i++) norm[i] = SQRT((SQR(spectrum[i]) + SQR(spectrum[size-i])));
+  norm[size/2] = spectrum[size/2];
 }
 
 void aubio_fft_getphas(smpl_t * phas, fft_data_t * spectrum, uint_t size) {