diff options
| author | Jamie Bullock <jamie@postlude.co.uk> | 2008-05-04 11:02:40 +0000 |
|---|---|---|
| committer | Jamie Bullock <jamie@postlude.co.uk> | 2008-05-04 11:02:40 +0000 |
| commit | 6a4be601746f6a84fc70434a61b16bea86b6c147 (patch) | |
| tree | fd8e8bfac0e6846f848a38490d583f4bfa7f6a6e /src/vector.c | |
| parent | d817fc9b42b0c04d6452159ca63f5a30834f5489 (diff) | |
| download | LibXtract-6a4be601746f6a84fc70434a61b16bea86b6c147.tar.gz LibXtract-6a4be601746f6a84fc70434a61b16bea86b6c147.tar.bz2 LibXtract-6a4be601746f6a84fc70434a61b16bea86b6c147.zip | |
Fixed bug in peak interpolation algorithm in xtract_peak_spectrum()
Diffstat (limited to 'src/vector.c')
| -rw-r--r-- | src/vector.c | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/vector.c b/src/vector.c index 449cd0d..b891111 100644 --- a/src/vector.c +++ b/src/vector.c @@ -436,7 +436,7 @@ int xtract_peak_spectrum(const float *data, const int N, const void *argv, float for(n = 1; n < N; n++){ if(input[n] >= threshold){ if(input[n] > input[n - 1] && input[n] > input[n + 1]){ - result[N + n] = q * (n + (p = .5 * (y = input[n-1] - + result[N + n] = q * (n + (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; |