app-pfrauenf/maxwell-twod-jump.hh

Go to the documentation of this file.

// jump constraints in 2D maxwell transmission problem

#ifndef maxwell_jump_hh
#define maxwell_jump_hh

#include <iostream>
#include <map>
#include "geometry/connector.hh"
#include "toolbox/dynArray.hh"
#include "constraints.hh"

using concepts::Real;

namespace hp2D {
  template<typename F>
  class Element;

  template<typename F>
  class Quad;
}

// ************************************************************* EpsJumpZero **

template<typename F>
class EpsJumpZero : public constraints::AnalyticalConstraint<F> {
public:
  EpsJumpZero(std::map<uint, concepts::Real2d>& coord,
        concepts::Attribute a1, Real eps1,
        concepts::Attribute a2, Real eps2);
  EpsJumpZero(const EpsJumpZero& z);
  virtual void lengthT(const concepts::Element<F>& e, uint& length);
  virtual void assembly(const concepts::Element<F>& e,
      concepts::TColumn<F>*& t, uint& n, uint& j,
      uint length);
  virtual constraints::AnalyticalConstraint<F>* clone() const;
  virtual void localMatrix(concepts::ElementMatrix<F>& em, uint& j,
         const constraints::Element<F>& elm) const;
  virtual void localVector(concepts::ElementMatrix<F>& em, uint& j,
         const constraints::Element<F>& elm) const;
protected:
  virtual std::ostream& info(std::ostream& os) const;
private:
  class CellData {
  public:
    CellData();
    CellData(const CellData& c);
    void addVertex(uint idx);
    uint vertex(uint i) const { return vtxIdx_[i]; }
    void addEdge(uint idx);
    uint edge() const { return edgIdx_; }
  private:
    uint vtxIdx_[2];
    uint edgIdx_;
  };


  concepts::Attribute attr_[2];
  Real eps_[2];

  std::map<uint, concepts::Real2d>& coordinates_;
  std::map<uint, CellData> indices_[2];

  const hp2D::Quad<Real>* getScalarQuad_(const concepts::Element<F>& e) const;
  bool compareCoord_(const concepts::Real2d& s1,
         const concepts::Real2d& s2) const;
  uint findCoord_(const concepts::Real2d& center) const;
};

#endif // maxwell_jump_hh