X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=libctsim%2Ffourier.cpp;h=149b4de234b5fc2c4738295ce04df5c51e792dd6;hp=a3a62d984ef80966bfa847445be2597ecc78be9f;hb=1a050c98763fbbc0662731b0b76953acede6f5d7;hpb=c8b19dfaffba9f06d8b6c40cb1bb83a8964867f7 diff --git a/libctsim/fourier.cpp b/libctsim/fourier.cpp index a3a62d9..149b4de 100644 --- a/libctsim/fourier.cpp +++ b/libctsim/fourier.cpp @@ -100,14 +100,14 @@ Fourier::shuffleNaturalToFourierOrder (ImageFile& im) } #ifdef HAVE_FFTW -void Fourier::shuffleNaturalToFourierOrder (fftw_complex* pVector, const int n) +void Fourier::shuffleNaturalToFourierOrder (fftw_complex* pVector, const int n) { fftw_complex* pTemp = static_cast(fftw_malloc(sizeof(fftw_complex) * n)); int i; - + if (isOdd(n)) { // Odd int iHalfN = (n - 1) / 2; - + pTemp[0][0] = pVector[iHalfN][0]; pTemp[0][1] = pVector[iHalfN][1]; for (i = 0; i < iHalfN; i++) { @@ -131,7 +131,7 @@ void Fourier::shuffleNaturalToFourierOrder (fftw_complex* pVector, const int n) pTemp[i + iHalfN][1] = pVector[i][1]; } } - + for (i = 0; i < n; i++) { pVector[i][0] = pTemp[i][0]; pVector[i][1] = pTemp[i][1]; @@ -145,7 +145,7 @@ void Fourier::shuffleFourierToNaturalOrder (fftw_complex* pVector, const int n) int i; if (isOdd(n)) { // Odd int iHalfN = (n - 1) / 2; - + pTemp[iHalfN][0] = pVector[0][0]; pTemp[iHalfN][1] = pVector[0][1]; for (i = 0; i < iHalfN; i++) { @@ -169,12 +169,12 @@ void Fourier::shuffleFourierToNaturalOrder (fftw_complex* pVector, const int n) pTemp[i + iHalfN + 1][1] = pVector[i+1][1]; } } - + for (i = 0; i < n; i++) { pVector[i][0] = pTemp[i][0]; pVector[i][1] = pTemp[i][1]; } - + fftw_free(pTemp); } #endif