example1r.cc

#include "Array.h"
#include "fftw++.h"

// Compile with:
// g++ -I .. -fopenmp example1r.cc ../fftw++.cc -lfftw3 -lfftw3_omp

using namespace std;
using namespace Array;
using namespace fftwpp;

int main()
{
  cout << "1D real to complex out-of-place FFT" << endl;
  
  fftw::maxthreads=get_max_threads();
  
  unsigned int n=4;
  unsigned int np=n/2+1;
  size_t align=sizeof(Complex);
  
  array1<double> f(n,align);
  array1<Complex> g(np,align);
  
  rcfft1d Forward(n,f,g);
  crfft1d Backward(n,g,f);
  
  for(unsigned int i=0; i < n; i++) f[i]=i;

  cout << endl << "input:" << endl << f << endl;
	
  Forward.fft(f,g);
  
  cout << endl << "output:" << endl << g << endl;
  
  Backward.fftNormalized(g,f);
  
  cout << endl << "back to input:" << endl << f << endl;
}