comparison dsp/mfcc/MFCC.cpp @ 29:1c9c4d2c0592

...

 
     /*
      * We calculate now mfccDCT matrix 
      * NB: +1 because of the DC component
      */
+
+    const double pi = 3.14159265358979323846264338327950288;
   
     for (i=0; i<nceps+1; i++) {
         for (j=0; j<totalFilters; j++) {
             mfccDCTMatrix[i][j] = (1./sqrt((double) totalFilters / 2.))  
-                * cos((double) i * ((double) j + 0.5) / (double) totalFilters * M_PI);
+                * cos((double) i * ((double) j + 0.5) / (double) totalFilters * pi);
         }
     }
 
     for (j=0;j<totalFilters;j++){
         mfccDCTMatrix[0][j]     =  (sqrt(2.)/2.) * mfccDCTMatrix[0][j];

     fftMag    = (double*)calloc(fftSize, sizeof(double)); 
     
     /* Zero-pad if needed */
     memcpy(inputData, inframe, length*sizeof(double));
 
-    window->cut(inputData);
+//!!!    window->cut(inputData);
   
     /* Calculate the fft on the input frame */
     FFT::process(fftSize, 0, inputData, imagIn, realOut, imagOut);
 
     for (i = 0; i < fftSize; ++i) {