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]