diff options
author | Mark Borgerding <mark@borgerding.net> | 2009-10-30 19:50:11 -0400 |
---|---|---|
committer | Mark Borgerding <mark@borgerding.net> | 2009-10-30 19:50:11 -0400 |
commit | a26b729cc9784eba82a5951be0e7738c1d209484 (patch) | |
tree | 05fee9b471125c9baf9808523f3166f9af86c6db /unsupported/Eigen/src/FFT | |
parent | 0fa68b9e50ba8b9e7eca9127ba5b10262e2e25d9 (diff) |
moved scaling to Eigen::FFT
Diffstat (limited to 'unsupported/Eigen/src/FFT')
-rw-r--r-- | unsupported/Eigen/src/FFT/ei_fftw_impl.h | 11 | ||||
-rw-r--r-- | unsupported/Eigen/src/FFT/ei_kissfft_impl.h | 10 |
2 files changed, 0 insertions, 21 deletions
diff --git a/unsupported/Eigen/src/FFT/ei_fftw_impl.h b/unsupported/Eigen/src/FFT/ei_fftw_impl.h index e1f67f334..cdb5861e8 100644 --- a/unsupported/Eigen/src/FFT/ei_fftw_impl.h +++ b/unsupported/Eigen/src/FFT/ei_fftw_impl.h @@ -187,12 +187,6 @@ void inv(Complex * dst,const Complex *src,int nfft) { get_plan(nfft,true,dst,src).inv(ei_fftw_cast(dst), ei_fftw_cast(src),nfft ); - - //TODO move scaling to Eigen::FFT - // scaling - Scalar s = Scalar(1.)/nfft; - for (int k=0;k<nfft;++k) - dst[k] *= s; } // half-complex to scalar @@ -200,11 +194,6 @@ void inv( Scalar * dst,const Complex * src,int nfft) { get_plan(nfft,true,dst,src).inv(ei_fftw_cast(dst), ei_fftw_cast(src),nfft ); - - //TODO move scaling to Eigen::FFT - Scalar s = Scalar(1.)/nfft; - for (int k=0;k<nfft;++k) - dst[k] *= s; } protected: diff --git a/unsupported/Eigen/src/FFT/ei_kissfft_impl.h b/unsupported/Eigen/src/FFT/ei_kissfft_impl.h index c068d8765..86ec7a1ca 100644 --- a/unsupported/Eigen/src/FFT/ei_kissfft_impl.h +++ b/unsupported/Eigen/src/FFT/ei_kissfft_impl.h @@ -334,7 +334,6 @@ void inv(Complex * dst,const Complex *src,int nfft) { get_plan(nfft,true).work(0, dst, src, 1,1); - scale(dst, nfft, Scalar(1)/nfft ); } // half-complex to scalar @@ -362,7 +361,6 @@ m_tmpBuf[k] = fek + fok; m_tmpBuf[ncfft-k] = conj(fek - fok); } - scale(&m_tmpBuf[0], ncfft, Scalar(1)/nfft ); get_plan(ncfft,true).work(0, reinterpret_cast<Complex*>(dst), &m_tmpBuf[0], 1,1); } } @@ -403,12 +401,4 @@ } return &twidref[0]; } - - // TODO move scaling up into Eigen::FFT - inline - void scale(Complex *dst,int n,Scalar s) - { - for (int k=0;k<n;++k) - dst[k] *= s; - } }; |