shithub: aubio

Info  •  Files  •  Log  •  Branches

src/spectral/fft.c: do not zero phase for no good reason, use ATAN2 macro

--- a/src/spectral/fft.c

+++ b/src/spectral/fft.c

@@ -185,9 +185,7 @@

   for (i = 0; i < spectrum->channels; i++) {

     spectrum->phas[i][0] = 0.;

     for (j=1; j < spectrum->length - 1; j++) {

-      if (compspec->data[i][j] == 0.) spectrum->phas[i][j] = 0;

-      else

-      spectrum->phas[i][j] = atan2f(compspec->data[i][compspec->length-j],

+      spectrum->phas[i][j] = ATAN2(compspec->data[i][compspec->length-j],

           compspec->data[i][j]);

     }

     spectrum->phas[i][spectrum->length-1] = 0.;