app-kersten/Maxwell2D_H.hh

Go to the documentation of this file.

#ifndef Maxwell2D_hh
#define Maxwell2D_hh

#include <set>
#include "basics/typedefs.hh"
#include "geometry/mesh.hh"
#include "space/formula.hh"
#include "formula/formula.hh"
#include "operator/sparseMatrix.hh"
#include "hp2D/bf_identity.hh"
#include "hp2D/bf_laplace.hh"
#include "hp2D/linearForm.hh"
#include "maxwellConstants.hh"
#include "Eddy2D_geometries.hh"
#include "Models.hh"

namespace hp2D {

  using concepts::Real;
  using concepts::Cmplx;

  class InputMaxwell2D_H;

  // *********************************************************** Maxwell2D_H **
  
  class Maxwell2D_H : public Models<Cmplx> {
  public:
    enum solverType { SUPERLU = 0, SUPERLU2 = 1, BICGSTAB = 2, BICGSTAB2 = 3,
          BICGSTABSUPERLU = 4};

    Maxwell2D_H(concepts::EddyGeometry2D& geom,
    enum solverType type = SUPERLU, bool diagPrecond = true,
    bool afterIter = false,
    const Real eps = EPS0, const Real omega = OMEGA50,
    const Real mu = MU0);
    Maxwell2D_H(concepts::EddyGeometry2D& geom, InputMaxwell2D_H& input);
    virtual ~Maxwell2D_H();
    Real dissipation();
    Real magnEnergy();
    virtual void graphicsOut(const uint numPoints,
           const std::string path = "");
  protected:
    virtual std::ostream& info(std::ostream& os) const;
    virtual concepts::Mesh2& mesh_() { return geom_.mesh(); }
    virtual const std::string mshAbbr_() { return geom_.meshAbbreviation(); }
  private:
    void constructor_();
    concepts::EddyGeometry2D& geom_;
    enum solverType type_;
    bool diagPrecond_;
    bool afterIter_;
    bool statusAfterIter_;
    std::auto_ptr<concepts::SparseMatrix<Cmplx> > A_;
    std::auto_ptr<concepts::SparseMatrix<Real> > M_;
    std::auto_ptr<concepts::SparseMatrix<Cmplx> > S_;
    // load vector
    std::auto_ptr<concepts::Vector<Cmplx> > rhs_;
    concepts::Vector<Cmplx>* rhsScaled_;
    std::auto_ptr<concepts::Vector<Cmplx> > residual_;
    std::auto_ptr<Real> residualNorm_;
    concepts::PiecewiseFormulaFun<Cmplx, Real> iOmegaEps_plus_Sigma_Inv_;
    const Real eps_;
    const Real omega_;
    const Real mu_;
    virtual void solve_();
    virtual void matrices_();
    void laplaceMatrix_();
    void identityMatrix_();
    void systemMatrix_();
    void linearform_();
    virtual void dataOut_();
    std::auto_ptr<Real> dissipation_;
    std::auto_ptr<Real> magnEnergy_;
  };


  // ****************************************************** InputMaxwell2D_H **

  class InputMaxwell2D_H : public concepts::InputParameter {
  public:
    InputMaxwell2D_H(concepts::InOutParameters& input);
    virtual std::ostream& letters(std::ostream& os) const;
    virtual std::ostream& arguments(std::ostream& os) const;
    virtual std::ostream& description(std::ostream& os) const;
    virtual int input(int opt, const char* optarg);
    enum Maxwell2D_H::solverType type() const { return type_; } 
    bool solving() const { return solving_; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    enum Maxwell2D_H::solverType type_;
    bool solving_;
  };

} // namespace hp2D

#endif // Maxwell2D_hh