#include <stdio.h>
#include <stdlib.h>
#include "Fourier.h"
#include "SphericalGrid.h"
#include "cmdLineParser.h"
#include <math.h>

int main(int argc,char* argv[]){
	// Read the parameters in from the command line
	cmdLineString In,Out;
	char* paramNames[]={"in","out"};
	cmdLineReadable* params[]={&In,&Out};
	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;
	}
	if(!Out.set){
		fprintf(stderr,"You must specify an output file\n");
		return EXIT_FAILURE;
	}

	SphericalGrid<> in,out,temp;
	FourierKeyS2<> key,tempKey;
	HarmonicTransform<> hForm;

	// Read in the input spherical function
	if(!in.read(In.value)){
		fprintf(stderr,"Failed to read in: %s\n",In.value);
		return EXIT_FAILURE;
	}

	// Kill off the 0th order frequency component and normalize so that
	// the function has unit norm.
	hForm.ForwardFourier(in,key);
	key(0,0)=0;
	hForm.InverseFourier(key,in);
	float l=float(sqrt(in.squareNorm()));
	for(int i=0;i<in.resolution();i++){for(int j=0;j<in.resolution();j++){in(i,j)/=l;}}

	// Allocate memory for the axial symmetry descriptor
	out.resize(in.resolution());

	// Compute the spherical harmonic transform
	hForm.ForwardFourier(in,key);

	// Allocate memory for the temporary key
	tempKey.resize(key.resolution());

	// Do your axial symmetry detection here, setting the axial symmetry
	// descriptor values into the SphericalGrid "out".
	// ...

	// Write out the axial symmetry descriptor here
	out.write(Out.value);
	return EXIT_SUCCESS;
}