- complex<double> fftSignal[m_nSignalPoints];
- complex<double> complexOutput[m_nSignalPoints];
- complex<double> filteredSignal[m_nSignalPoints];
- finiteFourierTransform (input, fftSignal, m_nSignalPoints, -1);
- dotProduct (m_vecFilter, fftSignal, filteredSignal, m_nSignalPoints);
- finiteFourierTransform (filteredSignal, complexOutput, m_nSignalPoints, 1);
- for (int i = 0; i < m_nSignalPoints; i++) {
- output[i] = abs( complexOutput[i] );
- }
- } else if (m_idFilterMethod == FILTER_METHOD_FFT || FILTER_METHOD_FFT_ZEROPAD_2 || FILTER_METHOD_FFT_ZEROPAD_4) {
- fftw_complex in[m_nFilterPoints], out[m_nFilterPoints];
- for (int i = 0; i < m_nSignalPoints; i++) {
- in[i].re = input[i];
- in[i].im = 0;
- }
- for (int i = m_nSignalPoints; i < m_nFilterPoints; i++) {
- in[i].re = in[i].im = 0; // ZeroPad
- }
- fftw_one(m_planForward, in, out);
- for (int i = 0; i < m_nFilterPoints; i++) {
- out[i].re = m_vecFilter[i] * out[i].re / m_nSignalPoints;
- out[i].im = m_vecFilter[i] * out[i].im / m_nSignalPoints;
- }
- fftw_one(m_planBackward, out, in);
+ double inputSignal[m_nFilterPoints];
+ for (int i = 0; i < m_nSignalPoints; i++)
+ inputSignal[i] = input[i];
+ for (int i = m_nSignalPoints; i < m_nFilterPoints; i++)
+ inputSignal[i] = 0; // zeropad
+ complex<double> fftSignal[m_nFilterPoints];
+ finiteFourierTransform (inputSignal, fftSignal, m_nFilterPoints, -1);
+ for (int i = 0; i < m_nFilterPoints; i++)
+ fftSignal[i] *= m_vecFilter[i];
+ double inverseFourier[m_nFilterPoints];
+ finiteFourierTransform (fftSignal, inverseFourier, m_nFilterPoints, 1);
+ for (int i = 0; i < m_nSignalPoints; i++)
+ output[i] = inverseFourier[i];
+ } else if (m_idFilterMethod == FILTER_METHOD_FOURIER_TABLE) {
+ double inputSignal[m_nFilterPoints];
+ for (int i = 0; i < m_nSignalPoints; i++)
+ inputSignal[i] = input[i];
+ for (int i = m_nSignalPoints; i < m_nFilterPoints; i++)
+ inputSignal[i] = 0; // zeropad
+ complex<double> fftSignal[m_nFilterPoints];
+ finiteFourierTransform (inputSignal, fftSignal, -1);
+ for (int i = 0; i < m_nFilterPoints; i++)
+ fftSignal[i] *= m_vecFilter[i];
+ double inverseFourier[m_nFilterPoints];
+ finiteFourierTransform (fftSignal, inverseFourier, 1);