From: Paul Brossier <piem@piem.org>
Date: Sun, 3 Mar 2013 03:48:49 +0000 (-0500)
Subject: src/mathutils.c: fix computation of gauss window
X-Git-Tag: 0.4.0-beta1~302
X-Git-Url: https://git.aubio.org/?a=commitdiff_plain;h=9e5622841f97cb424f98ba13d418bc4514f51340;p=aubio.git

src/mathutils.c: fix computation of gauss window
---

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++)