Changeset 6481c0c for src

Timestamp:

Sep 29, 2009, 7:58:43 PM (15 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, pitchshift, sampler, timestretch, yinfft+

Children:

a54502c

Parents:

78d14d4

Message:

src/temporal/filter.c: fix aubio_filter_do and aubio_filter_do_outplace

File:

• src/temporal/filter.c (modified) (2 diffs)

src/temporal/filter.c

@@ -37 +37 @@
 };
 
-void aubio_filter_do(aubio_filter_t * f, fvec_t * in) {
-  aubio_filter_do_outplace(f, in, in);
+void aubio_filter_do_outplace(aubio_filter_t * f, fvec_t * in, fvec_t * out) {
+  fvec_copy(in, out);
+  aubio_filter_do (f, out);
 }
 
-void aubio_filter_do_outplace(aubio_filter_t * f, fvec_t * in, fvec_t * out) {
+void aubio_filter_do(aubio_filter_t * f, fvec_t * in) {
   uint_t i,j,l, order = f->order;
   lsmp_t *x;
@@ -53 +54 @@
     for (j = 0; j < in->length; j++) {
       /* new input */
-      if (IS_DENORMAL(in->data[i][j])) {
-        x[0] = y[0] = 0.;
-      } else {
-        x[0] = in->data[i][j];
-        y[0] = b[0] * x[0];
-        for (l=1;l<order; l++) {
-          y[0] += b[l] * x[l];
-          y[0] -= a[l] * y[l];
-        }
+      x[0] = KILL_DENORMAL(in->data[i][j]);
+      y[0] = b[0] * x[0];
+      for (l=1;l<order; l++) {
+        y[0] += b[l] * x[l];
+        y[0] -= a[l] * y[l];
       }
       /* new output */
-      out->data[i][j] = y[0];
+      in->data[i][j] = y[0];
       /* store for next sample */
       for (l=order-1; l>0; l--){