estimator.hh

Go to the documentation of this file.

 
#ifndef aglowavEstimator_hh
#define aglowavEstimator_hh
 
#ifdef __GUNG__
  #pragma interface
#endif
 
#include "function/vector.hh"
#include "bem/space.hh"
#include "aglowav/space.hh"
 
namespace aglowav {
 
  // ***************************************************************** Est01 **
 
  template <class F = concepts::Real>
  class Est01 {
  public:
    Est01(bem::AdaptiveSpace<F>& spc, uint maxdim, concepts::Real trshld);
 
    bool refine() const {return ref_;}
 
    void operator()(const Haar3d<F>& wavspc, const concepts::Function<F>& f);
    void operator()(const Haar3d<F>& wavspc, const concepts::Vector<F>& f);
 
  private:
    bool ref_;
    uint maxdim_;
    concepts::Real trshld_;
    bem::AdaptiveSpace<F>& spc_;
  };
 
  template <class F>
  inline Est01<F>::Est01(bem::AdaptiveSpace<F>& spc,
       uint maxdim, concepts::Real trshld)
    : maxdim_(maxdim), trshld_(trshld), spc_(spc) {
    ref_ = (maxdim_ > spc_.dim()) ? 1 : 0;
  }
 
  // ***************************************************************** Est02 **
 
  template <class F = concepts::Real>
  class Est02 {
  public:
    Est02(bem::AdaptiveSpace<F>& spc,
    uint maxdim, concepts::Real trshld, concepts::Real s = 1.0);
 
    bool refine() const {return ref_;}
    void operator()(const Haar3d<F>& wavspc, const concepts::Function<F>& f);
    void operator()(const Haar3d<F>& wavspc, const concepts::Vector<F>& f);
  private:
    bool ref_;
    uint maxdim_;
    concepts::Real trshld2_;
    concepts::Real s_;
    bem::AdaptiveSpace<F>& spc_;
  };
 
  template <class F>
  inline Est02<F>::Est02(bem::AdaptiveSpace<F>& spc,
       uint maxdim, concepts::Real trshld, concepts::Real s)
    : maxdim_(maxdim), trshld2_(trshld*trshld), s_(s), spc_(spc) {
 
    ref_ = (maxdim_ > spc_.dim()) ? 1 : 0;
    if (s_ > 1.0)
      std::cout << "Uuups really s > 1.0 ??" << std::endl;
  }
 
  // ***************************************************************** Est03 **
 
  template <class F = concepts::Real>
  class Est03 {
  public:
    Est03(bem::AdaptiveSpace<F>& spc,
    uint maxdim, concepts::Real trshld);
 
    bool refine() const {return ref_;}
 
    void operator()(const Haar3d<F>& wavspc, const concepts::Function<F>& f);
    void operator()(const Haar3d<F>& wavspc, const concepts::Vector<F>& f);
  private:
    bool ref_;
    uint maxdim_;
    concepts::Real trshld_;
    bem::AdaptiveSpace<F>& spc_;
  };
 
  template <class F>
  inline Est03<F>::Est03(bem::AdaptiveSpace<F>& spc,
       uint maxdim, concepts::Real trshld)
    : maxdim_(maxdim), trshld_(trshld), spc_(spc) {
    ref_ = (maxdim_ > spc_.dim()) ? 1 : 0;
  }
 
  // ***************************************************************** Est04 **
 
  template <class F = concepts::Real>
  class Est04 {
  public:
    Est04(bem::AdaptiveSpace<F>& spc,
    uint maxdim, concepts::Real trshld, concepts::Real s = 1.0);
 
    bool refine() const {return ref_;}
    void operator()(const Haar3d<F>& wavspc, const concepts::Function<F>& f);
    void operator()(const Haar3d<F>& wavspc, const concepts::Vector<F>& f);
  private:
    bool ref_;
    uint maxdim_;
    concepts::Real trshld2_;
    concepts::Real s_;
    bem::AdaptiveSpace<F>& spc_;
  };
 
  template <class F>
  inline Est04<F>::Est04(bem::AdaptiveSpace<F>& spc,
       uint maxdim, concepts::Real trshld, concepts::Real s)
    : maxdim_(maxdim), trshld2_(trshld*trshld), s_(s), spc_(spc) {
    ref_ = (maxdim_ > spc_.dim()) ? 1 : 0;
  }
 
} // namespace aglowav
 
#endif // aglowavEstimator_hh