#include "header.h"

int main(void){
  double    x1[size], x2[size], x3[size];
  double 	x4[size][2];
  int 	i, n = 14000;
  long int  diff;
  double *  mat;
  double *	newmat[n];
  FILE   *	fp;
  
  mat = (double *) malloc(sizeof(double)*n*n);
  for(i=0; i<n; i++){
  	newmat[i] = mat + n*i;
  	}  
  
  clock_t   c0, c1; /* clock_t is defined on <time.h> and <sys/types.h> as int */
  // Using the FFTW timer ticks t = getticks();
  ticks 	t1, t2;
  
  /* Zeroing everything */
  for(i=0; i<size; i++){
  	x1[i]=0.0;
  	x2[i]=0.0;
  	x3[i]=0.0;
  	x4[i][1]=0.0; x4[i][2]=0.0;
  	}
  	
  	
  /* Reading in input */
  /*for(i=0; i<size; i++){
  	scanf("%lf", x1+i); // or &x1[i] or x1+i explain
  	}
  	*/
  fp = fopen("Test","r");
  for(i=0; i<size; i++){
  	fscanf(fp, "%lf", x1+i);
  	}
  fclose(fp);
  
  fp = fopen("TestToo.txt","w");
  for(i=0; i<size; i++){
  	fprintf(fp, "%e\n", x1[i]);
  	}
  fclose(fp);
  	
  /* Making copies: some of these calls are in place! */
  for(i=0; i<size; i++){
  	x2[i]=x1[i];
  	x3[i]=x1[i];
  	}
  
  //GSL Version	
  // add timing -0.045513
  gsl_complex_packed_array data = x1;
  
  t1 = getticks();
  i = gsl_fft_real_radix2_transform (data, 1, size); 
  t2 = getticks();
  printf("time %e (arbitrary units) in GSL\n", elapsed(t2,t1));
  /* Prototype is 
  int gsl_fft_real_radix2_transform (double data[], size_t stride, size_t n)
  */
  
  //FFTW version
  //See FFTW docs for prototype etc.
  fftw_plan plan;
  plan = fftw_plan_dft_r2c_1d(size, x2, x4, FFTW_MEASURE);
  
  t1 = getticks();
  fftw_execute(plan);
  t2 = getticks();
  printf("time %e (arbitrary units) in FFTW\n", elapsed(t2,t1));
  
  
  // Checking Ouput
  printf("GSL\n");
  for(i=0; i<size; i++){
  	printf("%d: %f \n", i, x1[i]);
  	}
  	
  printf("FFTW\n");
  for(i=0; i<size; i++){
  	printf("%d: %f \t %f \n", i, x4[i][1], x4[i][2]);
  	}
  	
  printf("%d \t %d \n", c0, c1);
  return 0;
  
  //Demo writing to a file!
}