code/include/NewWave/GSLEngine.hh

Go to the documentation of this file.

#ifndef NEWWAVE_GSL_ENGINE_HH
#define NEWWAVE_GSL_ENGINE_HH

#include <vector>

#include <gsl/gsl_wavelet.h>
#include <gsl/gsl_wavelet2d.h>

#include "NewWave/WaveletEngine.hh"
#include "NewWave/PixelDefinition.hh"

namespace NewWave{
  
  using std::size_t;
  using std::vector;
  
  class GSLEngine : public WaveletEngine{

  public:
    
    GSLEngine(const gsl_wavelet_type *type, size_t degree,
              const PixelDefinition &pixelDefn);
    
    ~GSLEngine();
    
    WaveletCoefficients transform(const PixelArray &pixels)const;
    
    PixelArray inverseTransform(const WaveletCoefficients &coefficients)const;

    
  private:
    
    void init(const PixelDefinition &pixelDefn);

    void initYLevels(const PixelDefinition &pixelDefn);

    void initPhiLevels();

    const gsl_wavelet_type *_waveletType;
    
    size_t _degree;
    
    gsl_wavelet *_wavelet;
    gsl_wavelet_workspace *_workspace;
    
    size_t _nSegments;
    size_t _nSegments2;
    size_t _nLevels;
    
    vector<size_t> _yLevels;
    vector<size_t> _phiLevels;
    vector<double> _yMinEdges;
    vector<double> _yMaxEdges;
    vector<double> _phiMinEdges;
    vector<double> _phiMaxEdges;
    
  };

}


#endif