Eddy2D_H_eField.hh

Go to the documentation of this file.

 
#ifndef Eddy2D_H_eField_hh
#define Eddy2D_H_eField_hh
 
#include "basics.hh"
#include "geometry/formula.hh"
#include "hp2D/function.hh"
#include "models/maxwellConstants.hh"
 
namespace hp2D {
 
  using concepts::Real;
  using concepts::Real2d;
  using concepts::Cmplx;
  using concepts::Cmplx2d;
 
  // ********************************************************* Eddy2D_eField **
 
  class Eddy2D_eField : public ElementFunction<2,Cmplx> {
  public:
    Eddy2D_eField(const concepts::PiecewiseFormulaBase<Real>& Sigma_Inv,
                  const concepts::PiecewiseFormulaBase<Real2d>* curlH0 = 0);
    virtual ~Eddy2D_eField();
 
    virtual Eddy2D_eField* clone() const {
      return new Eddy2D_eField(*Sigma_Inv_, curlH0_.get());
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    virtual void compute_(const Quad<>& elm,
        const concepts::Array<Cmplx>& coeff,
        concepts::Array<Cmplx>& val,
        const uint i[2]) const;
    virtual void compute_(const Quad<>& elm, const uint j[2],
        concepts::Array<Cmplx>& val,
        const uint i[2]) const;
    virtual void compute_(const Quad<>& elm,
        const concepts::Array<Cmplx>& coeff,
        concepts::Array<Cmplx>& val,
        const Real2d& p, const Real t = 0.0) const;
    void convert_(const Quad<>& elm, concepts::Array<Cmplx>& val,
      const Real2d& p, const Real t = 0.0) const;
    std::unique_ptr<concepts::PiecewiseFormulaBase<Real> > Sigma_Inv_;
    std::unique_ptr<const concepts::PiecewiseFormulaBase<Real2d> > curlH0_;
    PlCurl<Cmplx> plcurl_;
  };
 
  // **************************************************** Eddy2D_dissipation **
 
  class Eddy2D_dissipation : public ElementFunction<1,Cmplx> {
  public:
    Eddy2D_dissipation(const concepts::PiecewiseConstFormula<Real>& sigma,
           const concepts::ElementFunction<Cmplx>& eField);
    virtual ~Eddy2D_dissipation() {}
 
    virtual Eddy2D_dissipation* clone() const {
      return new Eddy2D_dissipation(*sigma_, *eField_);
    }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    virtual void compute_(const Quad<>& elm,
        const concepts::Array<Cmplx>& coeff,
        Cmplx& val, const uint i[2]) const;
    virtual void compute_(const Quad<>& elm, const uint j[2],
        Cmplx& val, const uint i[2]) const;
    virtual void compute_(const Quad<>& elm,
        const concepts::Array<Cmplx>& coeff, Cmplx& val,
        const Real2d& p, const Real t = 0.0) const;
    void convert_(const concepts::Array<Cmplx>& eField,
      const Real& sigma, Cmplx& val) const;
    std::unique_ptr<concepts::PiecewiseConstFormula<Real> > sigma_;
    std::unique_ptr<concepts::ElementFunction<Cmplx> > eField_;
  };
 
} // namespace hp2D
 
#endif // Eddy2D_H_eField_hh