# HG changeset patch
# User Jamie Bullock <jamie@jamiebullock.com>
# Date 1371684167 25200
# Node ID e678a849f2519530f2cd2f8c1f55ef0228961ea2
# Parent  4d6483986506f8a5c32b727164396beab1a11c76
Offset peak interpolation by one spectral bin. Fixes #20

diff -r 4d6483986506 -r e678a849f251 src/vector.c
--- a/src/vector.c	Wed Jun 19 16:05:59 2013 -0700
+++ b/src/vector.c	Wed Jun 19 16:22:47 2013 -0700
@@ -520,10 +520,11 @@
         {
             if(input[n] > input[n - 1] && n + 1 < N && input[n] > input[n + 1])
             {
-                result[N + n] = q * (n + (p = .5 * ((y = input[n-1]) -
+                result[N + n] = q * (n + 1 + (p = .5 * ((y = input[n-1]) -
                                                     (y3 = input[n+1])) / (input[n - 1] - 2 *
                                                             (y2 = input[n]) + input[n + 1])));
                 result[n] = y2 - .25 * (y - y3) * p;
+                printf("result[%d] = %f\tresult[N+n] = %f\n", n, result[n], result[n+N]);
             }
             else
             {