27 #define cexpf(x) (expf(crealf(x))*(cosf(cimagf(x))+sinf(cimagf(x))*I))
28 #warning Your C library does not have cexpf(). Please update it.
45 for (
int n =
LOGN; n --; )
47 y = (y << 1) | (x & 1);
61 for (
int n = 0; n <
N; n ++)
62 hamming[n] = 1 - 0.85 * cosf (2 * M_PI * n / N);
63 for (
int n = 0; n <
N; n ++)
65 for (
int n = 0; n < N / 2; n ++)
85 for (
int g = 0; g <
N; g += half << 1)
88 for (
int b = 0, r = 0;
b < half;
b ++, r += inv)
90 float complex even =
a[g +
b];
91 float complex odd =
roots[r] *
a[g + half +
b];
92 a[g +
b] = even + odd;
93 a[g + half +
b] = even - odd;
112 for (
int n = 0; n <
N; n ++)
119 for (
int n = 0; n < N / 2 - 1; n ++)
120 freq[n] = 2 * cabsf (a[1 + n]) /
N;
123 freq[N / 2 - 1] = cabsf (a[N / 2]) /
N;