Changeset 8979752 for src/spectral/fft.c

Timestamp:

Oct 11, 2009, 12:38:43 PM (15 years ago)

Author:

Paul Brossier <piem@piem.org>

Branches:

feature/autosink, feature/cnn, feature/cnn_org, feature/constantq, feature/crepe, feature/crepe_org, feature/pitchshift, feature/pydocstrings, feature/timestretch, fix/ffmpeg5, master, pitchshift, sampler, timestretch, yinfft+

Children:

4f4299d

Parents:

2fde783

Message:

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

File:

• src/spectral/fft.c (modified) (1 diff)

src/spectral/fft.c

@@ -184 +184 @@
   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.;
+    }
   }
 }