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

removed troublesome M_PI and M_PIl constants

Browse Source
This commit is contained in:
Mark Borgerding 2009-12-02 17:23:09 -05:00
parent e12f5adbde
commit e5c91b4e95
2 changed files with 4 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
unsupported/test/FFT.cpp
View File

@ -45,13 +45,10 @@ complex<long double> promote(long double x) { return complex<long double>( x);
long double totalpower=0; long double totalpower=0;
long double difpower=0; long double difpower=0;
cerr <<"idx\ttruth\t\tvalue\t|dif|=\n"; cerr <<"idx\ttruth\t\tvalue\t|dif|=\n";
long double pi = acos(-1);
for (size_t k0=0;k0<size_t(fftbuf.size());++k0) { for (size_t k0=0;k0<size_t(fftbuf.size());++k0) {
complex<long double> acc = 0; complex<long double> acc = 0;
#ifdef _GNU_SOURCE long double phinc = -2.*k0* pi / timebuf.size();
long double phinc = -2.*k0* M_PIl / timebuf.size();
#else
long double phinc = -2.*k0* M_PI / timebuf.size();
#endif
for (size_t k1=0;k1<size_t(timebuf.size());++k1) { for (size_t k1=0;k1<size_t(timebuf.size());++k1) {
acc += promote( timebuf[k1] ) * exp( complex<long double>(0,k1*phinc) ); acc += promote( timebuf[k1] ) * exp( complex<long double>(0,k1*phinc) );
} }

3
unsupported/test/FFTW.cpp
View File

@ -45,10 +45,11 @@ complex<long double> promote(long double x) { return complex<long double>( x);
{ {
long double totalpower=0; long double totalpower=0;
long double difpower=0; long double difpower=0;
long double pi = acos(-1);
cerr <<"idx\ttruth\t\tvalue\t|dif|=\n"; cerr <<"idx\ttruth\t\tvalue\t|dif|=\n";
for (size_t k0=0;k0<fftbuf.size();++k0) { for (size_t k0=0;k0<fftbuf.size();++k0) {
complex<long double> acc = 0; complex<long double> acc = 0;
long double phinc = -2.*k0* M_PIl / timebuf.size(); long double phinc = -2.*k0* pi / timebuf.size();
for (size_t k1=0;k1<timebuf.size();++k1) { for (size_t k1=0;k1<timebuf.size();++k1) {
acc += promote( timebuf[k1] ) * exp( complex<long double>(0,k1*phinc) ); acc += promote( timebuf[k1] ) * exp( complex<long double>(0,k1*phinc) );
} }