mirror/kissfft
1
0
Fork 0
mirror of https://github.com/mborgerding/kissfft.git synced 2025-05-27 21:20:27 -04:00
Code Issues Packages Projects Releases Wiki Activity

Add global pi constant as define

Browse Source
This commit is contained in:
Andreas Beeck 2019-04-12 10:49:00 +02:00
parent 36dbc05760
commit b59132b60f
3 changed files with 6 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
kiss_fft.c
View File

@ -349,8 +349,7 @@ kiss_fft_cfg kiss_fft_alloc(int nfft,int inverse_fft,void * mem,size_t * lenmem
st->inverse = inverse_fft; st->inverse = inverse_fft;
for (i=0;i<nfft;++i) { for (i=0;i<nfft;++i) {
const double pi=3.141592653589793238462643383279502884197169399375105820974944; double phase = -2*kiss_fft_pi*i / nfft;
double phase = -2*pi*i / nfft;
if (st->inverse) if (st->inverse)
phase *= -1; phase *= -1;
kf_cexp(st->twiddles+i, phase ); kf_cexp(st->twiddles+i, phase );

4
kiss_fft.h
View File

@ -56,6 +56,10 @@ extern "C" {
# endif # endif
#endif #endif
#ifndef kiss_fft_pi
#define kiss_fft_pi 3.141592653589793238462643383279502884197169399375105820974944
#endif
typedef struct { typedef struct {
kiss_fft_scalar r; kiss_fft_scalar r;
kiss_fft_scalar i; kiss_fft_scalar i;

2
tools/kiss_fftr.c
View File

@ -50,7 +50,7 @@ kiss_fftr_cfg kiss_fftr_alloc(int nfft,int inverse_fft,void * mem,size_t * lenme
for (i = 0; i < nfft/2; ++i) { for (i = 0; i < nfft/2; ++i) {
double phase = double phase =
-3.14159265358979323846264338327 * ((double) (i+1) / nfft + .5); -kiss_fft_pi * ((double) (i+1) / nfft + .5);
if (inverse_fft) if (inverse_fft)
phase *= -1; phase *= -1;
kf_cexp (st->super_twiddles+i,phase); kf_cexp (st->super_twiddles+i,phase);