QRealComplexFFT.cc

Go to the documentation of this file.

#include "QRealComplexFFT.hh"
 
ClassImp(Diana::QRealComplexFFT)
 
Q_BEGIN_NAMESPACE
 
double QRealComplexFFT::fTimer=0;
size_t QRealComplexFFT::fNCalls=0;
 
QRealComplexFFT::QRealComplexFFT() :
  fFFTW(NULL),
  fGSL(NULL)
{
#ifdef _HAVE_FFTW3_
  fFFTW=new QRealComplexFFTW3();
#else
  fGSL=new QRealComplexFFTGSL();
#endif
}
QRealComplexFFT::QRealComplexFFT(const size_t Size) :
  fFFTW(NULL),
  fGSL(NULL)
{
#ifdef _HAVE_FFTW3_
  fFFTW=new QRealComplexFFTW3(Size);
#else
  fGSL=new QRealComplexFFTGSL(Size);
#endif
}
QRealComplexFFT::~QRealComplexFFT() {
#ifdef _HAVE_FFTW3_
  delete fFFTW;
#else
  delete fGSL;
#endif
}
 
int QRealComplexFFT::TransformToFreq(const QVector& data, QVectorC& FFT,
                                     bool compress){
  int ret=0;
  clock_t t0=::clock();
#ifdef _HAVE_FFTW3_
  ret=fFFTW->TransformToFreq(data,FFT,compress);
#else
  ret=fGSL->TransformToFreq(data,FFT,compress);
#endif
  fTimer+=1.0E6*((UInt_t)(::clock()-t0))/CLOCKS_PER_SEC;
  fNCalls++;
  return ret;
}
 
int QRealComplexFFT::TransformFromFreq(const QVectorC& FT, QVector& spectrum,
                                       bool compress){
  int ret=0;
  clock_t t0=::clock();
#ifdef _HAVE_FFTW3_
  ret=fFFTW->TransformFromFreq(FT,spectrum,compress);
#else
  ret=fGSL->TransformFromFreq(FT,spectrum,compress);
#endif
  fTimer+=1.0E6*((UInt_t)(::clock()-t0))/CLOCKS_PER_SEC;
  fNCalls++;
  return ret;
}
 
void QRealComplexFFT::SetWindowType(WindowType wt, int coherent){
#ifdef _HAVE_FFTW3_
  return fFFTW->SetWindowType(wt,coherent);
#else
  return fGSL->SetWindowType(wt,coherent);
#endif
}
const char *QRealComplexFFT::GetFFTLibrary() {
#ifdef _HAVE_FFTW3_
  return "FFTW3";
#else
  return "GSL";
#endif
}
 
QRealComplexFFT *gRealComplexFFT = new QRealComplexFFT();
 
 
Q_END_NAMESPACE