1 files changed, 97 insertions, 0 deletions
diff --git a/unittest/tst_codec2_fft_init.c b/unittest/tst_codec2_fft_init.c
new file mode 100644
index 0000000..71ddc38
--- /dev/null
+++ b/unittest/tst_codec2_fft_init.c
@@ -0,0 +1,97 @@
+/*---------------------------------------------------------------------------*\
+
+  FILE........: tst_codec2_fft_init.c,
+  AUTHOR......: David Rowe, Don Reid
+  DATE CREATED: 30 May 2013, Oct 2018, Feb 2018
+
+  Test FFT Window initialization in Codec2_create
+
+\*---------------------------------------------------------------------------*/
+
+/*
+  Copyright (C) 2014 David Rowe
+
+  All rights reserved.
+
+  This program is free software; you can redistribute it and/or modify
+  it under the terms of the GNU Lesser General Public License version 2.1, as
+  published by the Free Software Foundation.  This program is
+  distributed in the hope that it will be useful, but WITHOUT ANY
+  WARRANTY; without even the implied warranty of MERCHANTABILITY or
+  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public
+  License for more details.
+
+  You should have received a copy of the GNU Lesser General Public License
+  along with this program; if not, see <http://www.gnu.org/licenses/>.
+*/
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <stdint.h>
+#include <math.h>
+
+#include "codec2.h"
+#include "codec2_internal.h"
+#include "defines.h"
+
+#include "machdep.h"
+
+static const float expect_w[] = {
+  0.004293, 0.004301, 0.004309, 0.004315,
+  0.004320, 0.004323, 0.004326, 0.004328,
+  0.004328, 0.004328, 0.004326, 0.004323,
+  0.004320, 0.004315, 0.004309, 0.004301};
+
+
+static const float expect_W[] = {
+ -0.002176,  0.002195,  0.004429, -0.008645,
+ -0.012196,  0.065359,  0.262390,  0.495616, 
+  0.601647,  0.495616,  0.262390,  0.065359, 
+ -0.012196, -0.008645,  0.004429,  0.002195};
+
+
+int float_cmp(float a, float b) {
+    if ( fabsf(a - b) < 1e-6f ) return 1;
+    else return 0;
+    }
+
+int main(int argc, char *argv[]) {
+
+    struct CODEC2 *codec2;
+    int i, j;
+
+    ////////
+    codec2 = codec2_create(CODEC2_MODE_700C);
+
+    j = (codec2->c2const.m_pitch / 2) - 8;
+    for (i=0; i<16; i++) {
+        printf("w[%d] = %f", j+i, 
+                (double)codec2->w[j+i]);
+        if (!float_cmp(codec2->w[j+i], expect_w[i])) {
+            printf(" Error, expected %f", (double)expect_w[i]);
+            }
+        printf("\n");
+        }
+
+    printf("\n");
+
+    j = (FFT_ENC / 2) - 8;
+    for (i=0; i<16; i++) {
+        printf("W[%d] = %f", j+i, 
+                (double)codec2->W[j+i]);
+        if (!float_cmp(codec2->W[j+i], expect_W[i])) {
+            printf(" Error, expected %f", (double)expect_W[i]);
+            }
+        printf("\n");
+        }
+
+    codec2_destroy(codec2);
+
+    printf("\nEnd of Test\n");
+    fclose(stdout);
+    fclose(stderr);
+
+    return(0);
+}
+
+/* vi:set ts=4 et sts=4: */