shithub: aubio

Download patch

ref: 897975259f6fdea0ea011c8f3f3f77652b968682
parent: 2fde783b5a4aff7e18c9dcabbb6cdef73d5a4b12
author: Paul Brossier <piem@piem.org>
date: Sun Oct 11 08:38:43 EDT 2009 

src/spectral/fft.c: 0 and N/2 + 1 phase must carry the sign of respective real components


--- a/src/spectral/fft.c
+++ b/src/spectral/fft.c
@@ -183,12 +183,20 @@
 void aubio_fft_get_phas(fvec_t * compspec, cvec_t * spectrum) {
   uint_t i, j;
   for (i = 0; i < spectrum->channels; i++) {
-    spectrum->phas[i][0] = 0.;
+    if (compspec->data[i][0] < 0) {
+      spectrum->phas[i][0] = PI;
+    } else {
+      spectrum->phas[i][0] = 0.;
+    }
     for (j=1; j < spectrum->length - 1; j++) {
       spectrum->phas[i][j] = ATAN2(compspec->data[i][compspec->length-j],
           compspec->data[i][j]);
     }
-    spectrum->phas[i][spectrum->length-1] = 0.;
+    if (compspec->data[i][compspec->length/2] < 0) {
+      spectrum->phas[i][spectrum->length - 1] = PI;
+    } else {
+      spectrum->phas[i][spectrum->length - 1] = 0.;
+    }
   }
 }