#include <stdio.h>
#include <stdlib.h>
#include <Fourier.h>
#include <cmdLineParser.h>


void ForwardFFT(const CircularArray<>& values,Complex<>* coeffs){
}
void InverseFFT(const Complex<>* coeffs,CircularArray<>& values){
}

int main(int argc,char* argv[]){
	// Read the parameters in from the command line
	cmdLineString In;
	char* paramNames[]={"in"};
	cmdLineReadable* params[]={&In};
	cmdLineParse(argc-1,&argv[1],paramNames,sizeof(paramNames)/sizeof(char*),params);
	if(!In.set){
		fprintf(stderr,"You must specify an input file\n");
		return EXIT_FAILURE;
	}

	int res;
	float l2;
	CircularArray<> cIn,cOut;
	Complex<>* coeffs;

	// Read in the array
	cIn.read(In.value);
	res=cIn.resolution();

	// Allocate space for the Fourier coefficients
	coeffs=new Complex<>[res];

	// Run the forward and inverse Fourier transforms
	ForwardFFT(cIn,coeffs);
	InverseFFT(coeffs,cOut);
	// Correct for the scaling term
	for(int i=0;i<res;i++){cOut(i)/=res;}



	// Test that the Fourier coefficients satisfy the conjugacy relations (difference should be zero)
	l2=0;
	for(int i=0;i<res;i++){l2+=(coeffs[i]-coeffs[(res-i)%res].conjugate()).squareNorm();}
	printf("Conjugate Test: %f\n",l2);



	// Compare the input and output (the difference should be zero)
	printf("Forward-Inverse Test: %f\n",CircularArray<>::SquareDifference(cIn,cOut));



	// Now compare the values of the Fourier coefficients (difference should be zero)
	FourierKey1D<> key;
	FourierTransform<> xForm;

	xForm.ForwardFourier(cIn,key);

	l2=0;
	float n=float(1.0/(2.0*PI))*key.resolution();
	for(int i=0;i<key.size();i++){l2+=(key(i)*n-coeffs[i]).squareNorm();}
	printf("Coefficient test: %f\n",l2);

	return EXIT_SUCCESS;
}