comparison qm-dsp-decimate/decimate.cpp @ 28:69ee50c19c0c tip

Add decimate-b

 /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
 
+#ifdef DECIMATE_B
+#include "qm-dsp/dsp/rateconversion/DecimatorB.h"
+#else
 #include "qm-dsp/dsp/rateconversion/Decimator.h"
+#endif
+
 #include "qm-dsp/maths/MathUtilities.h"
 
 #include <iostream>
 #include <sndfile.h>
 #include <time.h>

 
     if (!MathUtilities::isPowerOfTwo(factor)) {
         cerr << "ERROR: Factor must be a power of two" << endl;
         return 1;
     }
+#ifdef DECIMATE_B
+    if (factor < 2) {
+        cerr << "ERROR: Factor must be at least 2" << endl;
+        return 1;
+    }
+#else    
     if (factor < 2 || factor > Decimator::getHighestSupportedFactor()) {
         cerr << "ERROR: Only factors between 2 and "
              << Decimator::getHighestSupportedFactor() 
              << " inclusive are supported" <<endl;
         return 1;
     }
+#endif
     
     timeval tv;
     (void)gettimeofday(&tv, 0);
 
     char *fileName = strdup(argv[optind++]);

 	     << sf_strerror(sndfileOut) << endl;
 	return 1;
     }
 
     int channels = sfinfo.channels;
-    vector<Decimator *> decimators; // one per channel
-
     int ibs = 1024;
     int obs = ibs / factor;
 
+#ifdef DECIMATE_B
+    vector<DecimatorB *> decimators; // one per channel
+    for (int c = 0; c < channels; ++c) {
+	decimators.push_back(new DecimatorB(ibs, factor));
+    }
+#else
+    vector<Decimator *> decimators; // one per channel
     for (int c = 0; c < channels; ++c) {
 	decimators.push_back(new Decimator(ibs, factor));
     }
+#endif
 
     float *ibuf = new float[channels * ibs];
     float *obuf = new float[channels * obs];
     
     double **prebuf = new double*[channels];