47 # define SAMPPROD int64_t 48 #define SAMP_MAX 2147483647 51 # define SAMPPROD int32_t 52 #define SAMP_MAX 32767 55 #define SAMP_MIN -SAMP_MAX 57 #if defined(CHECK_OVERFLOW) 58 # define CHECK_OVERFLOW_OP(a,op,b) \ 59 if ( (SAMPPROD)(a) op (SAMPPROD)(b) > SAMP_MAX || (SAMPPROD)(a) op (SAMPPROD)(b) < SAMP_MIN ) { \ 60 fprintf(stderr,"WARNING:overflow @ " __FILE__ "(%d): (%d " #op" %d) = %ld\n",__LINE__,(a),(b),(SAMPPROD)(a) op (SAMPPROD)(b) ); } 64 # define smul(a,b) ( (SAMPPROD)(a)*(b) ) 65 # define sround( x ) (kiss_fft_scalar)( ( (x) + (1<<(FRACBITS-1)) ) >> FRACBITS ) 67 # define S_MUL(a,b) sround( smul(a,b) ) 69 # define C_MUL(m,a,b) \ 70 do{ (m).r = sround( smul((a).r,(b).r) - smul((a).i,(b).i) ); \ 71 (m).i = sround( smul((a).r,(b).i) + smul((a).i,(b).r) ); }while(0) 73 # define DIVSCALAR(x,k) \ 74 (x) = sround( smul( x, SAMP_MAX/k ) ) 76 # define C_FIXDIV(c,div) \ 77 do { DIVSCALAR( (c).r , div); \ 78 DIVSCALAR( (c).i , div); }while (0) 80 # define C_MULBYSCALAR( c, s ) \ 81 do{ (c).r = sround( smul( (c).r , s ) ) ;\ 82 (c).i = sround( smul( (c).i , s ) ) ; }while(0) 86 # define S_MUL(a,b) ( (a)*(b) ) 87 #define C_MUL(m,a,b) \ 88 do{ (m).r = (a).r*(b).r - (a).i*(b).i;\ 89 (m).i = (a).r*(b).i + (a).i*(b).r; }while(0) 90 # define C_FIXDIV(c,div) 91 # define C_MULBYSCALAR( c, s ) \ 93 (c).i *= (s); }while(0) 96 #ifndef CHECK_OVERFLOW_OP 97 # define CHECK_OVERFLOW_OP(a,op,b) 100 #define C_ADD( res, a,b)\ 102 CHECK_OVERFLOW_OP((a).r,+,(b).r)\ 103 CHECK_OVERFLOW_OP((a).i,+,(b).i)\ 104 (res).r=(a).r+(b).r; (res).i=(a).i+(b).i; \ 106 #define C_SUB( res, a,b)\ 108 CHECK_OVERFLOW_OP((a).r,-,(b).r)\ 109 CHECK_OVERFLOW_OP((a).i,-,(b).i)\ 110 (res).r=(a).r-(b).r; (res).i=(a).i-(b).i; \ 112 #define C_ADDTO( res , a)\ 114 CHECK_OVERFLOW_OP((res).r,+,(a).r)\ 115 CHECK_OVERFLOW_OP((res).i,+,(a).i)\ 116 (res).r += (a).r; (res).i += (a).i;\ 119 #define C_SUBFROM( res , a)\ 121 CHECK_OVERFLOW_OP((res).r,-,(a).r)\ 122 CHECK_OVERFLOW_OP((res).i,-,(a).i)\ 123 (res).r -= (a).r; (res).i -= (a).i; \ 128 # define KISS_FFT_COS(phase) floor(.5+SAMP_MAX * cos (phase)) 129 # define KISS_FFT_SIN(phase) floor(.5+SAMP_MAX * sin (phase)) 130 # define HALF_OF(x) ((x)>>1) 131 #elif defined(USE_SIMD) 132 # define KISS_FFT_COS(phase) _mm_set1_ps( cos(phase) ) 133 # define KISS_FFT_SIN(phase) _mm_set1_ps( sin(phase) ) 134 # define HALF_OF(x) ((x)*_mm_set1_ps(.5)) 136 # define KISS_FFT_COS(phase) (kiss_fft_scalar) cos(phase) 137 # define KISS_FFT_SIN(phase) (kiss_fft_scalar) sin(phase) 138 # define HALF_OF(x) ((x)*.5) 141 #define kf_cexp(x,phase) \ 143 (x)->r = KISS_FFT_COS(phase);\ 144 (x)->i = KISS_FFT_SIN(phase);\ 150 fprintf(stderr,"%g + %gi\n",(double)((c)->r),(double)((c)->i) ) 153 #ifdef KISS_FFT_USE_ALLOCA 159 #define KISS_FFT_TMP_ALLOC(nbytes) alloca(nbytes) 160 #define KISS_FFT_TMP_FREE(ptr) 162 #define KISS_FFT_TMP_ALLOC(nbytes) KISS_FFT_MALLOC(nbytes) 163 #define KISS_FFT_TMP_FREE(ptr) KISS_FFT_FREE(ptr)
int factors[2 *MAXFACTORS]