30 #define cexpf(x) (expf(crealf(x))*(cosf(cimagf(x))+sinf(cimagf(x))*I))
31 #warning Your C library does not have cexpf(). Please update it.
48 for (
int n =
LOGN; n --; )
50 y = (y << 1) | (x & 1);
64 for (
int n = 0; n <
N; n ++)
65 hamming[n] = 1 - 0.85 * cosf (2 * M_PI * n / N);
66 for (
int n = 0; n <
N; n ++)
68 for (
int n = 0; n < N / 2; n ++)
88 for (
int g = 0; g <
N; g += half << 1)
91 for (
int b = 0, r = 0;
b < half;
b ++, r += inv)
93 float complex even =
a[g +
b];
94 float complex odd =
roots[r] *
a[g + half +
b];
95 a[g +
b] = even + odd;
96 a[g + half +
b] = even - odd;
115 for (
int n = 0; n <
N; n ++)
122 for (
int n = 0; n < N / 2 - 1; n ++)
123 freq[n] = 2 * cabsf (a[1 + n]) /
N;
126 freq[N / 2 - 1] = cabsf (a[N / 2]) /
N;