/afs/in2p3.fr/throng/virgo/virgoDev/GWOLLUM/v1r3/src/Utils/FFT.h

Go to the documentation of this file.

//  Author : florent robinet (LAL - Orsay): robinet@lal.in2p3.fr
#ifndef __fft__
#define __fft__

#include "CUtils.h"
//#include "TMath.h"

// __float128 is not supported in CINT
#ifdef __CINT__
typedef struct {char a[7];} __float128; /* 7 chars have the same size as one __float128 */
#endif
#include "fftw3.h"

using namespace std;


class fft {
  
 public:
  
  fft(const int asize, const string aplan = "FFTW_ESTIMATE");
  
  virtual ~fft(void);// destructor
  bool Forward(double *aReDataIn=NULL, double *aImDataIn=NULL);

  bool Backward(double *aReDataIn=NULL, double *aImDataIn=NULL);

  double* GetNorm2(void);

  double* GetNorm(void);

  double* GetRe(void);

  double* GetIm(void);
 
 private:
  bool status_OK;         
  int fSize;              
  string fPlan;           

  fftw_complex* FFTin;    
  fftw_complex* FFTout;   
  fftw_plan FFT_Forward;  
  fftw_plan FFT_Backward; 

  int last_transform;     

};

#endif

---