comparison dsp/transforms/FFT.cpp @ 74:769da847732b

* A few fixes prompted by vamp-plugin-tester

 
 FFT::FFT(unsigned int n) :
     m_n(n),
     m_private(0)
 {
+    if (m_n & 1) {
+        std::cerr << "Error: Odd FFT size " << m_n
+                  << " not supported in this implementation"
+                  << std::endl;
+        return;
+    }
     KissFFTRec *rec = new KissFFTRec;
     rec->forward = kiss_fft_alloc(m_n, 0, 0, 0);
     rec->inverse = kiss_fft_alloc(m_n, 1, 0, 0);
     rec->in = new kiss_fft_cpx[m_n];
     rec->out = new kiss_fft_cpx[m_n];
     m_private = rec;
 }
 
 FFT::~FFT()
 {
+    if (!m_private) return;
     KissFFTRec *rec = (KissFFTRec *)m_private;
     kiss_fft_free(rec->forward);
     kiss_fft_free(rec->inverse);
     delete[] rec->in;
     delete[] rec->out;

 void
 FFT::process(bool inverse,
              const double *rin, const double *iin,
              double *rout, double *iout)
 {
+    if (!m_private) return;
     KissFFTRec *rec = (KissFFTRec *)m_private;
     for (int i = 0; i < m_n; ++i) {
         rec->in[i].r = rin[i];
     }
     if (iin) {

 
 FFTReal::FFTReal(unsigned int n) :
     m_n(n),
     m_private(0)
 {
+    if (m_n & 1) {
+        std::cerr << "Error: Odd FFT size " << m_n
+                  << " not supported in this implementation"
+                  << std::endl;
+        return;
+    }
     KissFFTRealRec *rec = new KissFFTRealRec;
     rec->forward = kiss_fftr_alloc(m_n, 0, 0, 0);
     rec->inverse = kiss_fftr_alloc(m_n, 1, 0, 0);
     rec->out = new kiss_fft_cpx[m_n];
     m_private = rec;
 }
 
 FFTReal::~FFTReal()
 {
+    if (!m_private) return;
     KissFFTRealRec *rec = (KissFFTRealRec *)m_private;
     kiss_fftr_free(rec->forward);
     kiss_fftr_free(rec->inverse);
     delete[] rec->out;
 }

 void
 FFTReal::process(bool inverse,
                  const double *rin,
                  double *rout, double *iout)
 {
+    if (!m_private) return;
     KissFFTRealRec *rec = (KissFFTRealRec *)m_private;
     if (inverse) {
         kiss_fftr(rec->inverse, rin, rec->out);
         for (int i = 0; i < m_n; ++i) {
             rout[i] = rec->out[i].r;