# HG changeset patch
# User Sean Enderby <sean.enderby@gmail.com>
# Date 1393246825 0
# Node ID fdbef1474be94ea425266b095ead1163ba6f3003
# Parent  15fde6ea8a99171ebfa74ddc0068c3ca9eb43ff0
correction to xtract_spectral_inharmonicicty

diff -r 15fde6ea8a99 -r fdbef1474be9 src/scalar.c
--- a/src/scalar.c	Tue Feb 18 20:18:44 2014 +0000
+++ b/src/scalar.c	Mon Feb 24 13:00:25 2014 +0000
@@ -623,7 +623,7 @@
 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;
 
@@ -635,7 +635,8 @@
     {
         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]);
         }
     }