ChTheo/FFT_and_FP_math
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

now printf is disabled if twiddles arrays are external. Now while used for emdedded there is no prfintf garbage in compiled code.

Browse Source
This commit is contained in:
Theodor Chikin
2025-10-09 15:38:10 +03:00
parent 86f488e456
commit 5e36a64acb
2 changed files with 10 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
C/FFT_FP_realisation.c
View File

@ -6,8 +6,11 @@
#ifndef FFT_FP_EXTERNAL_TWIDDLES #ifndef FFT_FP_EXTERNAL_TWIDDLES
static int64_t twiddle_re[TWIDDLE_L] = {0,}; static int64_t twiddle_re[TWIDDLE_L] = {0,};
static int64_t twiddle_im[TWIDDLE_L] = {0,}; static int64_t twiddle_im[TWIDDLE_L] = {0,};
#define PRINTF
#endif #endif
void fft_twiddle_gen(int64_t* tw_re, int64_t* tw_im){ void fft_twiddle_gen(int64_t* tw_re, int64_t* tw_im){
for (uint32_t k = 0; k < TWIDDLE_L; ++k){ for (uint32_t k = 0; k < TWIDDLE_L; ++k){
double angle = 2.0 * PI * k / DATA_L; double angle = 2.0 * PI * k / DATA_L;
@ -19,7 +22,9 @@ void fft_twiddle_gen(int64_t* tw_re, int64_t* tw_im){
void fft_fp_prepare(void){ void fft_fp_prepare(void){
fft_twiddle_gen(twiddle_re, twiddle_im); fft_twiddle_gen(twiddle_re, twiddle_im);
for (uint32_t k = 0; k < TWIDDLE_L; ++k){ for (uint32_t k = 0; k < TWIDDLE_L; ++k){
#ifdef PRINTF
printf("k, angle, tw_re, tw_im: %u %g %lld %lld\n", k, 2.0 * PI * k / DATA_L, (long long)twiddle_re[k], (long long)twiddle_im[k]); printf("k, angle, tw_re, tw_im: %u %g %lld %lld\n", k, 2.0 * PI * k / DATA_L, (long long)twiddle_re[k], (long long)twiddle_im[k]);
#endif
} }
} }
@ -90,7 +95,10 @@ void FFT_fp(int64_t* inp, uint32_t inp_L, int64_t* buf){
for (uint32_t ii = 0; ii < inp_L; ++ii){ for (uint32_t ii = 0; ii < inp_L; ++ii){
buf[ii*2] /= inp_L; buf[ii*2] /= inp_L;
buf[ii*2+1] /= inp_L; buf[ii*2+1] /= inp_L;
#ifdef PRINTF
printf("re,im: %lld, %lld\n", (long long)buf[ii*2], (long long)buf[ii*2 +1]); printf("re,im: %lld, %lld\n", (long long)buf[ii*2], (long long)buf[ii*2 +1]);
#endif
} }
} }

4
C/FFT_FP_realisation.h
View File

@ -20,8 +20,8 @@
#endif #endif
#ifdef FFT_FP_EXTERNAL_TWIDDLES #ifdef FFT_FP_EXTERNAL_TWIDDLES
extern int64_t twiddle_re[TWIDDLE_L]; extern volatile int64_t twiddle_re[TWIDDLE_L];
extern int64_t twiddle_im[TWIDDLE_L]; extern volatile int64_t twiddle_im[TWIDDLE_L];
#endif #endif
void fft_twiddle_gen(int64_t* twiddle_re, int64_t* twiddle_im); void fft_twiddle_gen(int64_t* twiddle_re, int64_t* twiddle_im);