+ complex<double> complexOutput[m_nSignalPoints];
+ complex<double> filteredSignal[m_nSignalPoints];
+ finiteFourierTransform (input, fftSignal, m_nSignalPoints, -1);
+ if (m_traceLevel >= TRACE_PLOT) {
+ double test[m_nSignalPoints];
+ for (int i = 0; i < m_nSignalPoints; i++)
+ test[i] = abs(fftSignal[i]);
+ ezplot_1d(test, m_nSignalPoints);
+ cio_kb_getc();
+ }
+ dotProduct (m_vecFilter, fftSignal, filteredSignal, m_nSignalPoints);
+ if (m_traceLevel >= TRACE_PLOT) {
+ double test[m_nSignalPoints];
+ for (int i = 0; i < m_nSignalPoints; i++)
+ test[i] = abs(filteredSignal[i]);
+ ezplot_1d(test, m_nSignalPoints);
+ cio_kb_getc();
+ }
+ 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);
+ if (m_traceLevel >= TRACE_PLOT) {
+ double test[m_nFilterPoints];
+ for (int i = 0; i < m_nFilterPoints; i++)
+ test[i] = sqrt(out[i].re * out[i].re + out[i].im * out[i].im);
+ ezplot_1d(test, m_nFilterPoints);
+ cio_kb_getc();
+ }
+ 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;
+ }
+ if (m_traceLevel >= TRACE_PLOT) {
+ double test[m_nFilterPoints];
+ for (int i = 0; i < m_nFilterPoints; i++)
+ test[i] = sqrt(out[i].re * out[i].re + out[i].im * out[i].im);
+ ezplot_1d(test, m_nFilterPoints);
+ cio_kb_getc();
+ }
+ fftw_one(m_planBackward, out, in);
+ if (m_traceLevel >= TRACE_PLOT) {
+ double test[m_nFilterPoints];
+ for (int i = 0; i < m_nFilterPoints; i++)
+ test[i] = sqrt(in[i].re * in[i].re + in[i].im * in[i].im);
+ ezplot_1d(test, m_nFilterPoints);
+ cio_kb_getc();
+ }
+ for (int i = 0; i < m_nSignalPoints; i++)
+ output[i] = sqrt (in[i].re * in[i].re + in[i].im * in[i].im);