trafo.hh

Go to the documentation of this file.

 
#ifndef aglowav2trafo_hh
#define aglowav2trafo_hh
 
#ifdef __GUNG__
#pragma interface
#endif
 
#include "aglowav2/space.hh"
#include "operator/compositions.hh"
#include "basics/output.hh"
 
namespace aglowav2 {
 
  // ******************************************************************* C2W **
 
  template<class F = concepts::Real>
  class C2W : public concepts::Operator<F> {
  public:
    C2W(const Haar3d0<F>& wavspc);
    ~C2W() {delete[] leven_;  delete[] lodd_;}
 
    void operator()(const concepts::Function<F>& fncC,
        concepts::Function<F>& fncW);
    void operator()(const concepts::Vector<F>& fncC,
                    concepts::Vector<F>& fncW);
    const Haar3d0<F>& spaceX() const {return wavspc_;}
    const concepts::Space<F>&  spaceY() const {return wavspc_.space();}
 
  protected:
    std::ostream& info(std::ostream& os) const;
 
  private:
    const Haar3d0<F>& wavspc_;
 
    F* leven_;
    F* lodd_;
  };
 
  template<class F>
  C2W<F>::C2W(const Haar3d0<F>& wavspc) 
    : concepts::Operator<F>(wavspc.dim(), wavspc.space().dim())
    , wavspc_(wavspc)
  {
    leven_ = new F[wavspc_.dim()];
    lodd_  = new F[wavspc_.dim()];
  }
 
  template<class F>
  std::ostream& C2W<F>::info(std::ostream& os) const {
    os << "aglowav2::" << concepts::typeOf(*this) << "(" << wavspc_.space() << " --> ";
    return os << wavspc_ << ')';
  }
 
  // ******************************************************************* W2C **
 
  template<class F = concepts::Real>
  class W2C : public concepts::Operator<F> {
  public:
    W2C(const Haar3d0<F>& wavspc);
    ~W2C() {delete[] ptrs_;  delete[] data_;  delete[] locdata_;}
 
    void operator()(const concepts::Function<F>& fncW,
        concepts::Function<F>& fncC);
    void operator()(const concepts::Vector<F>& fncW,
                    concepts::Vector<F>& fncC);
    const concepts::Space<F>& spaceX() const {return wavspc_.space();}
    const Haar3d0<F>& spaceY() const {return wavspc_;}
 
  protected:
    std::ostream& info(std::ostream& os) const;
 
  private:
    const Haar3d0<F>& wavspc_;
 
    uint maxdim_;
    uint nlvl_;
 
    F** ptrs_;
    F*  data_;
    F*  locdata_;
  };
 
  template<class F>
  W2C<F>::W2C(const Haar3d0<F>& wavspc)
    : concepts::Operator<F>(wavspc.space().dim(), wavspc.dim())
    , wavspc_(wavspc), maxdim_(wavspc_.maxLocDim()), nlvl_(wavspc_.nlvl() + 1)
  {
    ptrs_    = new F*[nlvl_];
    data_    = new F[maxdim_ * nlvl_];
    locdata_ = new F[maxdim_];
  }
 
  template<class F>
  std::ostream& W2C<F>::info(std::ostream& os) const {
    os << "aglowav2::" <<concepts::typeOf(*this)<<"(" << wavspc_ << " --> ";
    return os << wavspc_.space() << ')';
  }
 
} // namespace aglowav
 
#endif // aglowav2trafo_hh