3 #include <gsl/gsl_errno.h>
4 #include <gsl/gsl_fft_complex.h>
6 #define REAL(z,i) ((z)[2*(i)])
7 #define IMAG(z,i) ((z)[2*(i)+1])
12 int i; double data[2*128];
14 for (i = 0; i < 128; i++)
16 REAL(data,i) = 0.0; IMAG(data,i) = 0.0;
21 for (i = 1; i <= 10; i++)
23 REAL(data,i) = REAL(data,128-i) = 1.0;
26 for (i = 0; i < 128; i++)
28 printf ("%d %e %e\n", i,
29 REAL(data,i), IMAG(data,i));
33 gsl_fft_complex_radix2_forward (data, 1, 128);
35 for (i = 0; i < 128; i++)
37 printf ("%d %e %e\n", i,
38 REAL(data,i)/sqrt(128),
39 IMAG(data,i)/sqrt(128));