aglowav2/trafo.hh

Go to the documentation of this file.

/* Transformation operators for agglomerated wavelets.
 */

#ifndef aglowav2trafo_hh
#define aglowav2trafo_hh

#ifdef __GUNG__
#pragma interface
#endif

#include "aglowav2/space.hh"
#include "operator/compositions.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::C2W(" << 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::W2C(" << wavspc_ << " --> ";
    return os << wavspc_.space() << ')';
  }

} // namespace aglowav

#endif // aglowav2trafo_hh