linDG1D/elementPair.hh

Go to the documentation of this file.

// element pair for linear DG FEM in 1D

#ifndef elementPairDG1D_hh
#define elementPairDG1D_hh

#include "basics/typedefs.hh"
#include "space/elementPairs.hh"
#include "linearFEM/element1D.hh"

// debugging
#include "basics/debug.hh"

#define DGElementPairConstr_D 0

namespace linDG1D {

  using concepts::Real;

  // ********************************************************* DGElementPair **

  class DGElementPair : public concepts::ElementPair<Real> {
  public:
    DGElementPair(const linearFEM::Line& elm1, const uint idx1,
      const linearFEM::Line& elm2, const uint idx2) : 
      elm1_(elm1), elm2_(elm2), idx1_(idx1), idx2_(idx2),
      size1_(0), size2_(0) {
      DEBUGL(DGElementPairConstr_D, elm1 << ", " << elm2);
    }
    DGElementPair(const linearFEM::Line& elm, const uint idx,
      const Real size1, const Real size2) : 
      elm1_(elm), elm2_(elm), idx1_(idx), idx2_(idx),
      size1_(size1), size2_(size2) {
      DEBUGL(DGElementPairConstr_D, "Element 1:" << size1 << ", Element2:"
       << size2);
    }
    virtual const linearFEM::Line& elm1() const { return elm1_; }
    virtual const linearFEM::Line& elm2() const { return elm2_; }
    uint idx1() const { return idx1_; }
    uint idx2() const { return idx2_; }
    Real size1() const { return size1_; }
    Real size2() const { return size2_; }
  private:
    const linearFEM::Line& elm1_;
    const linearFEM::Line& elm2_;
    const uint idx1_;
    const uint idx2_;
    const Real size1_, size2_;
  };

} // namespace linDG1D

#endif // elementPairDG1D_hh