mirror/kissfft
1
0
Fork 0
mirror of https://github.com/mborgerding/kissfft.git synced 2025-07-16 03:54:22 -04:00
Code Issues Packages Projects Releases Wiki Activity

Create long double pi constant

Browse Source
This commit is contained in:
Chris Mumford
2025-07-11 13:59:41 -04:00
parent d12c0e1206
commit aff158e922
1 changed files with 7 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
test/testcpp.cc
View File

@ -7,6 +7,7 @@
*/ */
#include "kissfft.hh" #include "kissfft.hh"
#include <iostream> #include <iostream>
#include <cmath>
#include <cstdlib> #include <cstdlib>
#include <typeinfo> #include <typeinfo>
@ -41,9 +42,14 @@ void dotest(int nfft)
long double totalpower=0; long double totalpower=0;
long double difpower=0; long double difpower=0;
// Create long double constant for pi because M_PIl is not defined by
// all toolchains.
const long double pi = std::acosl(-1);
for (int k0=0;k0<nfft;++k0) { for (int k0=0;k0<nfft;++k0) {
complex<long double> acc = 0; complex<long double> acc = 0;
long double phinc = 2*k0* M_PI / nfft; long double phinc = 2*k0* pi / nfft;
for (int k1=0;k1<nfft;++k1) { for (int k1=0;k1<nfft;++k1) {
complex<long double> x(inbuf[k1].real(),inbuf[k1].imag()); complex<long double> x(inbuf[k1].real(),inbuf[k1].imag());
acc += x * exp( complex<long double>(0,-k1*phinc) ); acc += x * exp( complex<long double>(0,-k1*phinc) );