-
- ImageFileArrayConst vLHS = getArray();
- ImageFileArrayConst vLHSImag = getImaginaryArray();
- ImageFileArray vRealResult = result.getArray();
- ImageFileArray vImagResult = result.getImaginaryArray();
-
- unsigned int ix, iy;
-
- // alloc output matrix
- CTSimComplex** complexOut = new CTSimComplex* [m_nx];
- for (ix = 0; ix < m_nx; ix++)
- complexOut[ix] = new CTSimComplex [m_ny];
-
- // fourier each x column
- CTSimComplex* pY = new CTSimComplex [m_ny];
- for (ix = 0; ix < m_nx; ix++) {
- for (iy = 0; iy < m_ny; iy++) {
- double dImag = 0;
- if (isComplex())
- dImag = vLHSImag[ix][iy];
- pY[iy] = std::complex<double>(vLHS[ix][iy], dImag);
- }
- ProcessSignal::finiteFourierTransform (pY, complexOut[ix], m_ny, ProcessSignal::FORWARD);