/* This reads one of our sample transform data files, consisting of a table of
     65536 complex numbers assumed to have conjugate symmetry,
     and produces the inverse transform, which will be
     all real.
   Note that the input is read from stdin. */

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "complexdefs.h"
#define numcells 4096
void initialize_exps(int N);
void fft(Complex *sequence,int N,int mode);

main()
{
  Complex *sequence;
  float fdummy,min,max;
  int j,dummy;
 
  initialize_exps(numcells);
  sequence=malloc(numcells*sizeof(Complex));
  for(j=0;j<numcells;j++)
    {
      sequence[j].Im=0.0;
      scanf("%d%g%g",&dummy,&(sequence[j].Re),&(sequence[j].Im));
    }
  fft(sequence,numcells,-1);

 for(j=0;j<numcells;j++)
    printf("%d\t%g\n",j,sequence[j].Re);
  
}