Changeset 52c1de9 for tests/src/spectral

Timestamp:

Aug 22, 2018, 3:56:50 PM (6 years ago)

Author:

Paul Brossier <piem@piem.org>

Branches:

feature/autosink, feature/cnn, feature/cnn_org, feature/constantq, feature/crepe, feature/crepe_org, feature/pitchshift, feature/pydocstrings, feature/timestretch, fix/ffmpeg5, master

Children:

49b19d4

Parents:

3beacb3

Message:

tests/src/spectral/test-dct.c: check reconstruction works

File:

• tests/src/spectral/test-dct.c (modified) (4 diffs)

tests/src/spectral/test-dct.c

@@ -1 @@
-#include <aubio.h>
+#include <math.h>
+#include "aubio.h"
+#include "utils_tests.h"
 
 int main (void)
@@ -5 +7 @@
   int return_code = 0;
   uint_t win_s = 32; // window size
-  uint_t i, n_iters = 10; // number of iterations
+  uint_t i, j, n_iters = 10; // number of iterations
   // create dct object
   aubio_dct_t * dct = new_aubio_dct(win_s);
@@ -11 +13 @@
   fvec_t * in = new_fvec (win_s); // input buffer
   fvec_t * dctout = new_fvec (win_s); // output buffer
+  fvec_t * out = new_fvec (win_s); // input buffer
 
   if (!dct || !in || !dctout) {
@@ -20 +23 @@
   for (i = 0; i < n_iters; i++) {
     aubio_dct_do (dct, in, dctout);
-    aubio_dct_rdo (dct, dctout, in);
+    aubio_dct_rdo (dct, dctout, out);
+    for (j = 0; j < in->length; j++) {
+      if (fabsf(in->data[j] - out->data[j]) > 10.e-4) {
+        fprintf(stderr, "dct reconstruction failed\n");
+      }
+    }
   }
+
+  fvec_print(in);
   fvec_print(dctout);
-  fvec_print(in);
+  fvec_print(out);
+
   del_fvec(dctout);
   del_fvec(in);
+  del_fvec(out);
+  del_aubio_dct(dct);
 
-  del_aubio_dct(dct);
   return return_code;
 }