#include <iostream.h>
#include <console.h>

#include "transform.h"

int main(int argc, char *argv[]) {

argc = ccommand(&argv);

Array<Complex> v; cout << "\nEnter data vector..."; cin >> v;
double delta; cout << "\nEnter spacing..."; cin >> delta;
int forward; cout << "\nEnter 1/0 for FFT forward/back..."; cin >> forward;

double deltat; Array<Complex> vft;
fftcomp1d(v,delta,forward,vft,deltat);
cout << "vft: " << vft << endl << "deltat = " << deltat << endl;
int n = vft.n, n1 = n / 2;
for (int i=n1-1;i>=0;i--) cout << -(i+1) * deltat << ' ' << vft[n-1-i] << endl;
for (int i=0;i<n1;i++) cout << i * deltat << ' ' << vft[i] << endl;

return 0;
}