comparison audioio/PhaseVocoderTimeStretcher.cpp @ 19:f17798a555df

...

     else m_prevMag = 0;
     m_prevPercussiveCount = new int[m_channels];
     m_prevPercussive = false;
 
     m_dbuf = (float *)fftwf_malloc(sizeof(float) * m_wlen);
-    m_time = (fftwf_complex *)fftwf_malloc(sizeof(fftwf_complex) * m_wlen);
+    m_time = (float *)fftwf_malloc(sizeof(float) * m_wlen);
     m_freq = (fftwf_complex *)fftwf_malloc(sizeof(fftwf_complex) * m_wlen);
         
-    m_plan = fftwf_plan_dft_1d(m_wlen, m_time, m_freq, FFTW_FORWARD, FFTW_ESTIMATE);
-    m_iplan = fftwf_plan_dft_c2r_1d(m_wlen, m_freq, m_dbuf, FFTW_ESTIMATE);
+    m_plan = fftwf_plan_dft_r2c_1d(m_wlen, m_time, m_freq, FFTW_ESTIMATE);
+    m_iplan = fftwf_plan_dft_c2r_1d(m_wlen, m_freq, m_time, FFTW_ESTIMATE);
 
     m_inbuf = new RingBuffer<float> *[m_channels];
     m_outbuf = new RingBuffer<float> *[m_channels];
     m_mashbuf = new float *[m_channels];
 

     for (i = 0; i < m_wlen/2; ++i) {
 	float temp = buf[i];
 	buf[i] = buf[i + m_wlen/2];
 	buf[i + m_wlen/2] = temp;
     }
-    
+
     for (i = 0; i < m_wlen; ++i) {
-	m_time[i][0] = buf[i];
-	m_time[i][1] = 0.0;
+	m_time[i] = buf[i];
     }
 
     fftwf_execute(m_plan); // m_time -> m_freq
 
     if (m_sharpen && c == 0) { //!!!

         }
 
         m_prevPercussiveCount[c] = count;
     }
 
-    for (i = 0; i < m_wlen; ++i) {
+    for (i = 0; i < m_wlen; ++i) { //!!! /2
 
         float mag;
 
         if (m_sharpen && c == 0) {
             mag = m_prevMag[c][i]; // can reuse this

 	
         float expectedPhase = m_prevPhase[c][i] + omega;
 
         float phaseError = princargf(phase - expectedPhase);
 
-        float phaseIncrement = (omega + phaseError) / m_n1;
-
-        float adjustedPhase = m_prevAdjustedPhase[c][i] +
-            lastStep * phaseIncrement;
-
-        if (isPercussive) adjustedPhase = phase;
+        float adjustedPhase = phase;
+
+        if (!isPercussive) {
+//            if (fabsf(phaseError) < (1.1f * (lastStep * M_PI) / m_wlen)) {
+
+                float phaseIncrement = (omega + phaseError) / m_n1;
+
+                adjustedPhase = m_prevAdjustedPhase[c][i] +
+                    lastStep * phaseIncrement;
+//            }
+        }
+
+//        if (isPercussive) adjustedPhase = phase;
 	
 	float real = mag * cosf(adjustedPhase);
 	float imag = mag * sinf(adjustedPhase);
 	m_freq[i][0] = real;
 	m_freq[i][1] = imag;
 
         m_prevPhase[c][i] = phase;
         m_prevAdjustedPhase[c][i] = adjustedPhase;
     }
     
-    fftwf_execute(m_iplan); // m_freq -> in, inverse fft
-    
+    fftwf_execute(m_iplan); // m_freq -> m_time, inverse fft
+
     for (i = 0; i < m_wlen/2; ++i) {
-	float temp = buf[i] / m_wlen;
-	buf[i] = buf[i + m_wlen/2] / m_wlen;
-	buf[i + m_wlen/2] = temp;
-    }
-    
-    m_window->cut(buf);
+	float temp = m_time[i];
+	m_time[i] = m_time[i + m_wlen/2];
+	m_time[i + m_wlen/2] = temp;
+    }
 
     for (i = 0; i < m_wlen; ++i) {
-        out[i] += buf[i];
+	m_time[i] = m_time[i] / m_wlen;
+    }
+
+    m_window->cut(m_time);
+
+    for (i = 0; i < m_wlen; ++i) {
+        out[i] += m_time[i];
     }
 
     if (modulation) {
 
         float area = m_window->getArea();