face.hh

Go to the documentation of this file.

 
#ifndef HP3D_FACE_HH
#define HP3D_FACE_HH
 
#include "basics/typedefs.hh"
#include "basics/outputOperator.hh"
#include "basics/Zm.hh"
#include "geometry/normalVector.hh"
 
 
namespace hp2D {
  // forward declaration
  template <class F>
  class Quad;
}
 
namespace concepts {
  // forward declarations
  class Quad3d;
  template <class F>
  class TMatrix;
}
 
namespace hp3D {
 
  using concepts::Real;
  using concepts::Z2;
  
  // forward declaration
  class Hexahedron;
  
  // ******************************* HexahedronFaceBase *************************
  
  class HexahedronFaceBase : public concepts::OutputOperator {
  public:
    HexahedronFaceBase(const concepts::FaceNormalVectorRule& 
                       normalVectorRule = concepts::FaceNormalVectorRule())
      : normalVectorRule_(normalVectorRule.clone())
    {} 
    
    hp2D::Quad<Real>* face(const Hexahedron &elm, const ushort k, 
                           bool emptyElm = true,
                           hp2D::Quad<Real>* face = 0);
    
  protected:
    virtual std::ostream& info(std::ostream &os) const;
    
    virtual Real weight_(bool first, Z2 dir) const = 0;
    
    std::unique_ptr<concepts::FaceNormalVectorRule> normalVectorRule_;
    
  private:
    void transferDof_(const Hexahedron &elm, const ushort k,
                      const bool first, const Z2 dir,
                      hp2D::Quad<Real> &face);
    
    bool faceDof_(const Hexahedron &elm, const ushort k, 
                  const uint dofHex, uint *dofFace);
    
  };
    
  // **************************************** HexahedronFaceFirst ***************
  
  class HexahedronFaceFirst : public HexahedronFaceBase {
  public:
    HexahedronFaceFirst(const concepts::FaceNormalVectorRule&
                        normalVectorRule = concepts::FaceNormalVectorRule())
      : HexahedronFaceBase(normalVectorRule)
    {
    }
    
  protected:
    virtual std::ostream& info(std::ostream &os) const;
  private:
    virtual Real weight_(bool first, Z2 dir) const;
  };
    
 
  
}
 
#endif // HP3D_FACE_HH