Maxwell2D_E_hField.hh

Go to the documentation of this file.

 
#ifndef Maxwell2D_E_hField_hh
#define Maxwell2D_E_hField_hh
 
#include "geometry/formula.hh"
#include "hp2Dedge/function.hh"
#include "hp2Dedge/bilinearForm.hh"
#include "basics.hh"
#include "toolbox.hh"
 
namespace hp2D {
 
  using concepts::Real;
  using concepts::Cmplx;
 
  class Maxwell2D_hField : 
    public hp2D::ElementFunction<1, Cmplx, hp2Dedge::Quad<> > {
  public:
    Maxwell2D_hField(const Real& omega, const Real& mu);
    virtual ~Maxwell2D_hField() {}
 
    virtual Maxwell2D_hField* clone() const {
      return new Maxwell2D_hField(1.0 / std::imag(factor_), 1.0); }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    virtual void compute_(const hp2Dedge::Quad<>& elm,
        const concepts::Array<Cmplx>& coeff,
        Cmplx& val, const uint i[2]) const;
    virtual void compute_(const hp2Dedge::Quad<>& elm, const uint j[3],
        Cmplx& val, const uint i[2]) const;
    virtual void compute_(const hp2Dedge::Quad<>& elm,
        const concepts::Array<Cmplx>& coeff,
        Cmplx& val,
        const Real2d& p, const Real t = 0.0) const;
    hp2Dedge::Rot<Cmplx> scurl_;
    Cmplx factor_;
  };
 
} // namespace hp2D
 
#endif // Maxwell2D_E_hField_hh