comparison dsp/transforms/FFT.cpp @ 170:8c86761a5533

Fix overrun in reading inverse complex-to-real FFT input (contrary to docs)

         delete[] io;
     }
 
     void inverse(const double *ri, const double *ii, double *ro) {
 
-        for (int i = 0; i < m_n; ++i) {
+        // kiss_fftr.h says
+        // "input freqdata has nfft/2+1 complex points"
+
+        for (int i = 0; i < m_n/2 + 1; ++i) {
             m_c[i].r = ri[i];
             m_c[i].i = ii[i];
         }
         
         kiss_fftri(m_plani, m_c, ro);