# HG changeset patch
# User Sean Enderby <sean.enderby@gmail.com>
# Date 1393337549 0
# Node ID a2f66a305e53ab243bb1e127a6a2de64799717de
# Parent  cccd314664ac61066312fcde9ec67d520182ad7a
tristimulus3 correction

diff -r cccd314664ac -r a2f66a305e53 src/scalar.c
--- a/src/scalar.c	Tue Feb 25 13:53:39 2014 +0000
+++ b/src/scalar.c	Tue Feb 25 14:12:29 2014 +0000
@@ -486,21 +486,21 @@
 
 int xtract_tristimulus_3(const double *data, const int N, const void *argv, double *result)
 {
+    int n = N >> 1, h = 0, i;
+    double den = 0.0, num = 0.0, fund = 0.0, temp = 0.0;
+    const double *freqs;
 
-    int n = N, count = 0;
+    fund = *(double *)argv;
+    freqs = data + n;
 
-    double den, num, temp;
-
-    den = num = temp = 0.0;
-
-    for(n = 0; n < N; n++)
+    for(i = 0; i < n; i++)
     {
-        if((temp = data[n]))
+        if((temp = data[i]))
         {
             den += temp;
-            if(count >= 5)
+            h = round(freqs[i] / fund);
+            if(h >= 5)
                 num += temp;
-            count++;
         }
     }