diff options
author | Mark Borgerding <mark@borgerding.net> | 2009-10-30 20:26:30 -0400 |
---|---|---|
committer | Mark Borgerding <mark@borgerding.net> | 2009-10-30 20:26:30 -0400 |
commit | d659fd9b148870e3c1e367139bec388142f2818e (patch) | |
tree | 19aa8beb5a2f572f2ea73178d28719c2477b1b6b /unsupported/Eigen/src/FFT | |
parent | a26b729cc9784eba82a5951be0e7738c1d209484 (diff) |
moved real-half-spectrum reflection into Eigen::FFT
Diffstat (limited to 'unsupported/Eigen/src/FFT')
-rw-r--r-- | unsupported/Eigen/src/FFT/ei_fftw_impl.h | 3 | ||||
-rw-r--r-- | unsupported/Eigen/src/FFT/ei_kissfft_impl.h | 2 |
2 files changed, 0 insertions, 5 deletions
diff --git a/unsupported/Eigen/src/FFT/ei_fftw_impl.h b/unsupported/Eigen/src/FFT/ei_fftw_impl.h index cdb5861e8..18473a29b 100644 --- a/unsupported/Eigen/src/FFT/ei_fftw_impl.h +++ b/unsupported/Eigen/src/FFT/ei_fftw_impl.h @@ -177,9 +177,6 @@ void fwd( Complex * dst,const Scalar * src,int nfft) { get_plan(nfft,false,dst,src).fwd(ei_fftw_cast(dst), ei_fftw_cast(src) ,nfft); - int nhbins=(nfft>>1)+1; - for (int k=nhbins;k < nfft; ++k ) - dst[k] = conj(dst[nfft-k]); } // inverse complex-to-complex diff --git a/unsupported/Eigen/src/FFT/ei_kissfft_impl.h b/unsupported/Eigen/src/FFT/ei_kissfft_impl.h index 86ec7a1ca..859e7e6c9 100644 --- a/unsupported/Eigen/src/FFT/ei_kissfft_impl.h +++ b/unsupported/Eigen/src/FFT/ei_kissfft_impl.h @@ -322,8 +322,6 @@ // place conjugate-symmetric half at the end for completeness // TODO: make this configurable ( opt-out ) - for ( k=1;k < ncfft ; ++k ) - dst[nfft-k] = conj(dst[k]); dst[0] = dc; dst[ncfft] = nyquist; } |