return split_radix_permutation(i, m) * 4 - 1;
}
-/** 1 / sqrt(2). */
-#define SQRTHALF (float)0.70710678118654752440
-
#define BF(x, y, a, b) {\
x = a - b;\
y = a + b;\
BF(z[6].re, z[2].re, z[2].re, t7);
BF(z[6].im, z[2].im, z[2].im, t8);
- TRANSFORM(z[1], z[3], z[5], z[7], SQRTHALF, SQRTHALF);
+ TRANSFORM(z[1], z[3], z[5], z[7], M_SQRT1_2, M_SQRT1_2);
}
static void fft16(struct fft_complex *z)
fft4(z + 12);
TRANSFORM_ZERO(z[0], z[4], z[8], z[12]);
- TRANSFORM(z[2], z[6], z[10], z[14], SQRTHALF, SQRTHALF);
+ TRANSFORM(z[2], z[6], z[10], z[14], M_SQRT1_2, M_SQRT1_2);
TRANSFORM(z[1], z[5], z[9], z[13], cos_16[1], cos_16[3]);
TRANSFORM(z[3], z[7], z[11], z[15], cos_16[3], cos_16[1]);
}