hpExtension.hh

Go to the documentation of this file.

 
#ifndef hpExtension3D_hh
#define hpExtension3D_hh
 
#include "basics/typedefs.hh"
#include "basics/exceptions.hh"
#include "space/postProcess.hh"
 
namespace concepts {
  // forward declaration
  template<typename F>
  class Vector;
 
  template<typename F>
  class Array;
}
 
namespace hp3D {
  using concepts::Real;
 
  // ******************************************************** RefineOrRaise  **
 
  template<typename F, typename G = typename concepts::Realtype<F>::type>
  class RefineOrRaise : public concepts::CellPostprocess<F> {
  public:
    RefineOrRaise(const concepts::Vector<F>& solution, const Real theta,
                  concepts::AdaptiveSpace<G,
                  concepts::AdaptiveAdjustP<3> >& spc);
    virtual void operator() (const concepts::Element<F>& elm);
    virtual void operator() (const concepts::Cell& cell) {}
  private:
    concepts::Array<F> coeff_;
    concepts::Array<F> legendre_;
    concepts::Array<F> regr_;
    const concepts::Vector<F>& solution_;
    concepts::AdaptiveSpace<G, concepts::AdaptiveAdjustP<3> >& spc_;
    const Real theta_;
    Real transCoeff_(const uint i, const uint j) const;
    bool regression_(concepts::Array<F>& coeff, const uint p) const;
  };
 
} // namespace hp3D
 
#endif // hpExtension3D_hh