das2/dft.h

Go to the documentation of this file.

#ifndef _das2_dft_h_
#define _das2_dft_h_

typedef struct das2_dft_t{

   /* FFTW variables */
   void* vpPlan;
   void* vpIn;
   void* vpOut;
   
   /* Input vector length */
   size_t uLen;
   
   /* Input vector is real only*/
   bool bRealOnly;
   
   /* Holder for the window function and name*/
   char* sWindow;
   double* pWnd;
   
   /* Holder for the magnitude result */
   bool bNewMag;
   double* pMag;
   size_t uMagLen;
   
   /* Holder for continuous real and imaginary results */
   bool bNewCmp[2];   /* fftw convention: 0 = reals, 1 = img */
   double* pCmpOut[2];
   size_t uCmpLen[2];
   
} Das2Dft;

Das2Dft* new_Dft(size_t uLen, const char* sWindow);

void del_Dft(Das2Dft* pThis);

ErrorCode Dft_calculate(
   Das2Dft* pThis, const double* pReal, const double* pImg, size_t uLen
);

const double* Dft_getReal(Das2Dft* pThis, size_t* pLen);

const double* Dft_getImg(Das2Dft* pThis, size_t* pLen);

const double* Dft_getMagnitude(Das2Dft* pThis, size_t* pLen);

typedef struct das2_psd_t{
   
   /* FFTW variables */
   void* vpPlan;
   void* vpIn;
   void* vpOut;
   
   /* Input vector information */
   size_t uLen;
   bool bRealOnly;
   
   /* Center data about average first */
   bool bCenter;
   
   /* Holder for up conversion arrays, helps Psd_calculate_f*/
   size_t uUpConvLen;
   double* pUpConvReal;
   double* pUpConvImg;
   
   /* Holder for the window function and name */
   char* sWindow;
   double* pWnd;
   double rWndSqSum;
   
   /* Holder for the PSD result */
   double* pMag;
   size_t uMagLen;
   
   /* Total Energy calculations */
   double rPwrIn;
   double rPwrOut;
   
} Das2Psd;

Das2Psd* new_Psd(size_t uLen,  bool bCenter, const char* sWindow);

void del_Das2Psd(Das2Psd* pThis);


ErrorCode Psd_calculate(
   Das2Psd* pThis, const double* pReal, const double* pImg, size_t uLen
);

ErrorCode Psd_calculate_f(
   Das2Psd* pThis, const float* pReal, const float* pImg, size_t uLen
);

double Psd_powerRatio(const Das2Psd* pThis, double* pInput, double* pOutput);


const double* Psd_get(const Das2Psd* pThis, size_t* pLen);

#endif

---