elementPair.hh

Go to the documentation of this file.

 
#ifndef dgElementPair_hh
#define dgElementPair_hh
 
#include "basics/typedefs.hh"
#include "space/elementPairs.hh"
#include "element.hh"
 
#include <vector>
 
#define DEBUG_FILLDATA 0
#define DEBUG_FILLDATA_CONSISTENCY 0
 
namespace linDG3D {
  using concepts::Real;
 
  // ******************************************************* FvdgElementPair **
  
  class FvdgElementPair : public concepts::ElementPair<Real> {
  public:
    FvdgElementPair(const FvdgElement& elm1,
        const FvdgElement& elm2,
        uint faceIdx1, uint faceIdx2);
    FvdgElementPair(const FvdgElement& elm, uint faceIdx);
    virtual ~FvdgElementPair();
    
    virtual const FvdgElement& elm1() const { return elm1_; }
    virtual const FvdgElement& elm2() const { return elm2_; }
    
    Real area() const { return area_; }
    Real diam() const { return diam_; }
    concepts::Real3d normal() const { return normal_; }
    uint numPairs() const { return edgePairs_.size(); }
    uint edgeIndex( uint elem, uint edge ) const;
    void quadrature(std::vector<Real>& weights,
        std::vector<concepts::Real3d>& xi0,
        std::vector<concepts::Real3d>& xi1) const;
    void quadrature(std::vector<Real>& weights,
        std::vector<concepts::Real3d>& xi0,
        std::vector<concepts::Real3d>& xi1,
        std::vector<concepts::Real3d>& x) const;
    uint faceIdx1() const { return faceIdx1_; };
    uint faceIdx2() const { return faceIdx2_; };
  private:
    const FvdgElement& elm1_;
    const FvdgElement& elm2_;
    Real area_;
    Real diam_;
    concepts::Real3d normal_;
    uint faceIdx1_;
    uint faceIdx2_;
    static std::vector<Real> weights_;
    static std::vector<concepts::Real3d> zetaAbsc_;
    static concepts::Real3d xiEdge_[6];
    std::vector<std::pair<uint, uint> > edgePairs_;
    void fillData_();
    void fillQuadrature_();
    concepts::Real3d
    barycentric_(concepts::Real3d zeta, concepts::Real3d* xVtx) const;
    void edgVtx_(concepts::Real3d* xEdge, concepts::Real3d* xVtx) const;
  };
  
} // namespace linDG3D
 
#endif // dgElementPair_hh