comparison src/scalar.c @ 223:fdbef1474be9

correction to xtract_spectral_inharmonicicty

 }
 
 int xtract_spectral_inharmonicity(const double *data, const int N, const void *argv, double *result)
 {
 
-    int n = N >> 1;
+    int n = N >> 1, h = 0;
     double num = 0.0, den = 0.0, fund;
     const double *freqs, *amps;
 
     fund = *(double *)argv;
     amps = data;

 
     while(n--)
     {
         if(amps[n])
         {
-            num += fabs(freqs[n] - n * fund) * XTRACT_SQ(amps[n]);
+            h = round(freqs[n] / fund);
+            num += fabs(freqs[n] - h * fund) * XTRACT_SQ(amps[n]);
             den += XTRACT_SQ(amps[n]);
         }
     }
 
     *result = (2 * num) / (fund * den);