#include <stdio.h>
#include <complex>
#include <fftw3.h>

int main(){
	fftw_complex *in, *out;
	fftw_plan p;
	int N = 100;
	
	in = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * N);
	out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * N);
	for(int i=0;i<N;i++){
		(*(in+i))[0] = i;
		(*(in+i))[1] = i;
	}
	p = fftw_plan_dft_1d(N, in, out, FFTW_FORWARD, FFTW_ESTIMATE);
	for(int i=0;i<N;i++)
		printf("(%lf,%lf)\n",(*(out+i))[0], (*(out+i))[1] );
	fftw_execute(p);
	for(int i=0;i<N;i++)
		printf("(%lf,%lf)\n",(*(out+i))[0], (*(out+i))[1] );
}