#include "LArFFTWPlan.h"

Public Member Functions
	LArFFTWPlan (int transformSize, const std::string &option)

	~LArFFTWPlan ()

Public Attributes
void *	fPlan

void *	rPlan

void *	fIn

void *	fOut

void *	rIn

void *	rOut

Private Member Functions
unsigned int	MapFFTWOption ()

Private Attributes
int	fSize

int	fFreqSize

int *	fN

std::string	fOption

Static Private Attributes
static std::mutex	mutex_

Detailed Description

Definition at line 13 of file LArFFTWPlan.h.

Constructor & Destructor Documentation

util::LArFFTWPlan::LArFFTWPlan	(	int	transformSize,
		const std::string &	option
	)

Definition at line 6 of file LArFFTWPlan.cxx.

References fFreqSize, fIn, fN, fOut, fPlan, fSize, MapFFTWOption(), mutex_, rIn, rOut, and rPlan.

   : fSize(transformSize), fOption(option)
 {
 
   std::lock_guard<std::mutex> lock(mutex_);
 
   fFreqSize = fSize / 2 + 1;
   fN = new int[1];
   fN[0] = fSize;
 
   fIn = fftw_malloc(sizeof(double) * fSize);
   fOut = fftw_malloc(sizeof(fftw_complex) * fFreqSize);
   fPlan = (void*)fftw_plan_dft_r2c(1, fN, (double*)fIn, (fftw_complex*)fOut, MapFFTWOption());
 
   rIn = fftw_malloc(sizeof(fftw_complex) * fFreqSize);
   rOut = fftw_malloc(sizeof(double) * fSize);
   rPlan = (void*)fftw_plan_dft_c2r(1, fN, (fftw_complex*)rIn, (double*)rOut, MapFFTWOption());
 }

util::LArFFTWPlan::~LArFFTWPlan ( )

Definition at line 25 of file LArFFTWPlan.cxx.

References fIn, fN, fOut, fPlan, rIn, rOut, and rPlan.

 {
   fftw_destroy_plan((fftw_plan)fPlan);
   fPlan = 0;
   fftw_free(fIn);
   fIn = 0;
   fftw_free((fftw_complex*)fOut);
   fOut = 0;
 
   fftw_destroy_plan((fftw_plan)rPlan);
   rPlan = 0;
   fftw_free((fftw_complex*)rIn);
   rIn = 0;
   fftw_free(rOut);
   rOut = 0;
 
   delete[] fN;
   fN = 0;
 }

Member Function Documentation

unsigned int util::LArFFTWPlan::MapFFTWOption ( )

private

Definition at line 45 of file LArFFTWPlan.cxx.

References fOption.

Referenced by LArFFTWPlan().

 {
   std::transform(fOption.begin(), fOption.end(), fOption.begin(), ::toupper);
   if (fOption.find("ES") != string::npos) return FFTW_ESTIMATE;
   if (fOption.find("M") != string::npos) return FFTW_MEASURE;
   if (fOption.find("P") != string::npos) return FFTW_PATIENT;
   if (fOption.find("EX") != string::npos) return FFTW_EXHAUSTIVE;
   return FFTW_ESTIMATE;
 }