hp2D/edge.hh

Go to the documentation of this file.

#ifndef hp2dedge_hh
#define hp2dedge_hh

#include <iostream>
#include <memory>
#include "basics/typedefs.hh"
#include "basics/outputOperator.hh"
#include "basics/vectorsMatrices.hh"
#include "geometry/cell1D.hh"
#include "geometry/elementMaps.hh"
#include "geometry/integral.hh"
#include "integration/quadRule.hh"
#include "space/tmatrix.hh"
#include "space/hpMethod.hh"
#include "space/integral.hh"
#include "hp1D/element.hh"
#include "integration/karniadakis.hh"

namespace hp2D {

  using concepts::Real;

  // forward declaration
  template<typename F>
  class Quad;


  // ************************************************************ QuadEdgeBase **

  class QuadEdgeBase : public concepts::OutputOperator {
  public:
    template<class F>
    std::pair<hp1D::Element<F>*,concepts::Edge2d*> 
    edge(const Quad<F>& elm, const ushort k, bool emptyElm = true,
         hp1D::Element<F>* edge = 0);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
    
    virtual Real weight_(bool first, bool rho) const = 0;
  private:
    template<class F>
    void transferDof_(const ushort* p, const ushort k, const bool rho,
                      const bool first, const concepts::TMatrix<F>& T,
                      hp1D::Element<F>& edge);
    
    bool edgeDof_(const ushort* p, const ushort k, const bool rho,
                  const uint dofQuad, uint& dofEdge);
  };  

  class QuadEdgeFirst : public QuadEdgeBase {
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    virtual Real weight_(bool first, bool rho) const;
  };

  class QuadEdgeMean : public QuadEdgeBase {
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    virtual Real weight_(bool first, bool rho) const;
  };

  class QuadEdgeJump : public QuadEdgeBase {
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    virtual Real weight_(bool first, bool rho) const;
  };


} // namespace hp2D

#endif // hp2dedge_hh