/* fftd.h */ int bitrev(int a, int k); int ilog2(int n); void fft (double b[], double a1[], double a2[], int n, int sgn);