From: Paul Brossier <piem@piem.org>
Date: Sun, 11 Oct 2009 10:47:43 +0000 (+0200)
Subject: src/spectral/fft.c: fix imag boundaries, including for odd fft sizes
X-Git-Tag: bzr2git~140
X-Git-Url: http://git.tremily.us/?a=commitdiff_plain;h=74a4865ea20987d9adf612ca88f4c832905fffb8;p=aubio.git

src/spectral/fft.c: fix imag boundaries, including for odd fft sizes
---

diff --git a/src/spectral/fft.c b/src/spectral/fft.c
index 71acca64..b79fc064 100644
--- a/src/spectral/fft.c
+++ b/src/spectral/fft.c
@@ -224,7 +224,7 @@ void aubio_fft_get_norm(fvec_t * compspec, cvec_t * spectrum) {
 void aubio_fft_get_imag(cvec_t * spectrum, fvec_t * compspec) {
   uint_t i, j;
   for (i = 0; i < compspec->channels; i++) {
-    for (j = 1; j < compspec->length / 2 + 1; j++) {
+    for (j = 1; j < ( compspec->length + 1 ) / 2 /*- 1 + 1*/; j++) {
       compspec->data[i][compspec->length - j] =
         spectrum->norm[i][j]*SIN(spectrum->phas[i][j]);
     }
@@ -234,7 +234,7 @@ void aubio_fft_get_imag(cvec_t * spectrum, fvec_t * compspec) {
 void aubio_fft_get_real(cvec_t * spectrum, fvec_t * compspec) {
   uint_t i, j;
   for (i = 0; i < compspec->channels; i++) {
-    for (j = 0; j< compspec->length / 2 + 1; j++) {
+    for (j = 0; j < compspec->length / 2 + 1; j++) {
       compspec->data[i][j] = 
         spectrum->norm[i][j]*COS(spectrum->phas[i][j]);
     }