wavelet.hh

#ifndef __WT_WAVELET_HH__
#define __WT_WAVELET_HH__

#include "types.hh"

namespace wt {

class WaveletObj: public Object
{
protected:
  WaveletObj();

public:
  virtual ~WaveletObj();
  virtual CScalar evalAnalysis(const Scalar &t) const = 0;
  virtual CScalar evalSynthesis(const Scalar &t) const = 0;
  virtual double cutOffTime() const = 0;
  virtual double cutOffFreq() const = 0;
};


class Wavelet : public Container
{
public:
  typedef WaveletObj ObjectType;

public:
  Wavelet();
  Wavelet(WaveletObj *obj);
  Wavelet(const Wavelet &other);
  virtual ~Wavelet();

  Wavelet &operator=(const Wavelet &other);

  inline CScalar evalAnalysis(const Scalar &t) const {
    return _wavelet->evalAnalysis(t);
  }

  inline CScalar evalSynthesis(const Scalar &t) const {
    return _wavelet->evalSynthesis(t);
  }

  inline double cutOffTime() const {
    return _wavelet->cutOffTime();
  }

  inline double cutOffFreq() const {
    return _wavelet->cutOffFreq();
  }

protected:
  WaveletObj *_wavelet;
};


class MorletObj: public WaveletObj
{
public:
  explicit MorletObj(double dff=2);
  virtual ~MorletObj();
  virtual CScalar evalAnalysis(const Scalar &t) const;
  virtual CScalar evalSynthesis(const Scalar &t) const;
  virtual double cutOffTime() const;
  virtual double cutOffFreq() const;

protected:
  double _dff;
};


class Morlet: public Wavelet
{
public:
  Morlet(double dff=2.0);
  explicit Morlet(MorletObj *obj);
  Morlet(const Morlet &other);
  virtual ~Morlet();
  Morlet &operator=(const Morlet &other);

protected:
  MorletObj *_morlet;
};


class CauchyObj: public WaveletObj
{
public:
  explicit CauchyObj(double alpha=1);
  virtual ~CauchyObj();
  virtual CScalar evalAnalysis(const Scalar &t) const;
  virtual CScalar evalSynthesis(const Scalar &t) const;
  virtual double cutOffTime() const;
  virtual double cutOffFreq() const;

protected:
  double _alpha;
  double _norm;
};

class Cauchy: public Wavelet
{
public:
  explicit Cauchy(double alpha=1.0);
  Cauchy(CauchyObj *obj);
  Cauchy(const Cauchy &other);
  virtual ~Cauchy();
  Cauchy &operator=(const Cauchy &other);

protected:
  CauchyObj *_cauchy;
};

}

#endif // __WT_WAVELET_HH__