From: Paul Brossier <piem@piem.org>
Date: Sat, 3 Nov 2007 18:03:58 +0000 (+0100)
Subject: src/fft.c: fix norm computation when HAVE_COMPLEX_H is undefined
X-Git-Tag: bzr2git~481
X-Git-Url: http://git.tremily.us/?a=commitdiff_plain;h=ac5f22ca55c59c732ac71f7c256fa21cf2cf349a;p=aubio.git

src/fft.c: fix norm computation when HAVE_COMPLEX_H is undefined
---

diff --git a/src/fft.c b/src/fft.c
index 5cbc4780..e5c1a037 100644
--- 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) {