28 #define cexpf(x) (expf(crealf(x))*(cosf(cimagf(x))+sinf(cimagf(x))*I))
29 #warning Your C library does not have cexpf(). Please update it.
46 for (
int n =
LOGN; n --; )
48 y = (y << 1) | (x & 1);
62 for (
int n = 0; n <
N; n ++)
63 hamming[n] = 1 - 0.85 * cosf (2 * M_PI * n / N);
64 for (
int n = 0; n <
N; n ++)
66 for (
int n = 0; n < N / 2; n ++)
86 for (
int g = 0; g <
N; g += half << 1)
89 for (
int b = 0, r = 0;
b < half;
b ++, r += inv)
91 float complex even =
a[g +
b];
92 float complex odd =
roots[r] *
a[g + half +
b];
93 a[g +
b] = even + odd;
94 a[g + half +
b] = even - odd;
113 for (
int n = 0; n <
N; n ++)
120 for (
int n = 0; n < N / 2 - 1; n ++)
121 freq[n] = 2 * cabsf (a[1 + n]) /
N;
124 freq[N / 2 - 1] = cabsf (a[N / 2]) /
N;