From 9e5622841f97cb424f98ba13d418bc4514f51340 Mon Sep 17 00:00:00 2001
From: Paul Brossier <piem@piem.org>
Date: Sat, 2 Mar 2013 22:48:49 -0500
Subject: [PATCH] src/mathutils.c: fix computation of gauss window

---
 src/mathutils.c | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

diff --git a/src/mathutils.c b/src/mathutils.c
index 38431c97..69a79c60 100644
--- a/src/mathutils.c
+++ b/src/mathutils.c
@@ -104,8 +104,16 @@ new_aubio_window (char_t * window_type, uint_t size)
           - 0.01168 * COS(3.0*TWO_PI*i/(size-1.0));
       break;
     case aubio_win_gaussian:
-      for (i=0;i<size;i++)
-        w[i] = EXP(- 1.0 / SQR(size) * SQR(2.0*i-size));
+      {
+        lsmp_t a, b, c = 0.5;
+        uint_t n;
+        for (n = 0; n < size; n++)
+        {
+          a = (n-c*(size-1))/(SQR(c)*(size-1));
+          b = -c*SQR(a);
+          w[n] = EXP(b);
+        }
+      }
       break;
     case aubio_win_welch:
       for (i=0;i<size;i++)
-- 
2.11.0