GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-06-04 18:54:00 +08:00
Code Issues Packages Projects Releases Wiki Activity

fix pi in kissfft

Browse Source
This commit is contained in:
Alexey Korepanov 2024-08-02 22:57:47 +00:00 committed by Rasmus Munk Larsen
parent 2b7b7aac57
commit ec18dd09c8
1 changed files with 7 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
unsupported/Eigen/src/FFT/ei_kissfft_impl.h
View File

@ -27,12 +27,14 @@ struct kiss_cpx_fft {
std::vector<Complex> m_scratchBuf; std::vector<Complex> m_scratchBuf;
bool m_inverse; bool m_inverse;
static const Scalar m_pi4; // constant pi / 4
inline void make_twiddles(int nfft, bool inverse) { inline void make_twiddles(int nfft, bool inverse) {
using numext::cos; using numext::cos;
using numext::sin; using numext::sin;
m_inverse = inverse; m_inverse = inverse;
m_twiddles.resize(nfft); m_twiddles.resize(nfft);
double phinc = 0.25 * double(EIGEN_PI) / nfft; Scalar phinc = m_pi4 / nfft;
Scalar flip = inverse ? Scalar(1) : Scalar(-1); Scalar flip = inverse ? Scalar(1) : Scalar(-1);
m_twiddles[0] = Complex(Scalar(1), Scalar(0)); m_twiddles[0] = Complex(Scalar(1), Scalar(0));
if ((nfft & 1) == 0) m_twiddles[nfft / 2] = Complex(Scalar(-1), Scalar(0)); if ((nfft & 1) == 0) m_twiddles[nfft / 2] = Complex(Scalar(-1), Scalar(0));
@ -279,6 +281,10 @@ struct kiss_cpx_fft {
} }
}; };
template <typename _Scalar>
const typename kiss_cpx_fft<_Scalar>::Scalar kiss_cpx_fft<_Scalar>::m_pi4 =
numext::atan(kiss_cpx_fft<_Scalar>::Scalar(1));
template <typename Scalar_> template <typename Scalar_>
struct kissfft_impl { struct kissfft_impl {
typedef Scalar_ Scalar; typedef Scalar_ Scalar;