void FFT_inverse(int N, double *data) { int n = N/2; double norm = 0.0; int i=0; FFT_transform_internal(N, data, +1); /* Normalize */ norm=1/((double) n); for(i=0; i<N; i++) data[i] *= norm; }
void FFT_transform(int N, double *data) { FFT_transform_internal(N, data, -1); }
// MODIFIED FOR TEST CODE void FFT_transform(int N, double *data, int print) { FFT_transform_internal(N, data, -1, print); }