From 78a4b5718ecfa7cb91457bc9b4deff543ec10373 Mon Sep 17 00:00:00 2001
From: Jamie Bullock <jamie@jamiebullock.com>
Date: Tue, 3 Jun 2014 21:17:51 +0100
Subject: Add xtract_peak() determines whether the 'current' value is a peak

---
 src/scalar.c | 33 +++++++++++++++++++++++++++++++++
 1 file changed, 33 insertions(+)

(limited to 'src/scalar.c')

diff --git a/src/scalar.c b/src/scalar.c
index e0d2578..54cb459 100644
--- a/src/scalar.c
+++ b/src/scalar.c
@@ -1002,3 +1002,36 @@ int xtract_midicent(const double *data, const int N, const void *argv, double *r
     return XTRACT_SUCCESS;
 }
 
+int xtract_peak_picker(const double *data, const int N, const void *argv, double *result)
+{
+    double threshold = *(double *)argv;
+    double current = data[N - 1];
+    double average = 0.0;
+    double maximum = -DBL_MAX;
+    
+    for (uint32_t n = 0; n < N; ++n)
+    {
+        average += data[n];
+        if (data[n] > maximum)
+        {
+            maximum = data[n];
+        }
+    }
+    
+    average /= (double)N;
+    
+    if (current != maximum)
+    {
+        return XTRACT_NO_RESULT;
+    }
+    
+    if (current < average + threshold)
+    {
+        return XTRACT_NO_RESULT;
+    }
+    
+    return XTRACT_SUCCESS;
+    
+}
+
+
-- 
cgit v1.2.3