Eddy2D_H_Ring.hh

Go to the documentation of this file.

 
#ifndef MeshImp2DTest_hh
#define MeshImp2DTest_hh
 
#include <iostream>
#include <cstdlib>
#include "space/space.hh"
#include "basics/testcase.hh"
#include "basics/exceptions.hh"
#include "models/maxwell.hh"
 
namespace test {
 
  // ******************************************************** Eddy2DRingTest **
 
  class Eddy2DRingTest : public TestCase {
  public:
    Eddy2DRingTest(const uint pmax = 3, const Real sigma = 1,
       const Real epsilon = 1e-8, const bool graphics = false);
    virtual ~Eddy2DRingTest() {}
    virtual void run();
 
 
 
    void testCircle();
    void testMaxwell_PMC_H_Circle(concepts::EddyGeometry2D& geom);
    void testMaxwell_PEC_H_Circle(concepts::EddyGeometry2D& geom);
    void testMaxwell_PMC_E_Circle(concepts::EddyGeometry2D& geom);
    void testMaxwell_PEC_E_Circle(concepts::EddyGeometry2D& geom);
    void testMaxwell_PMC_H_Circle_DD(concepts::EddyGeometry2D& geom);
 
    void testRing();
    void testMaxwell_PMC_H_Ring(concepts::EddyGeometry2D& geom);
    void testMaxwell_PMC_H_Ring_DD(concepts::EddyGeometry2D& geom);
 
    template<class F>
    void testEnergies(F& model, const std::string& modelStr,
          const concepts::Array<Real>& P,
          const concepts::Array<Real>& W,
          const concepts::Array<uint>& Dim,
          const uint nelm);
    void testMaxwellRingSolvers();
  private:
    bool originalConstants_() const;
 
    concepts::InOutParameters output_;
 
    const uint pmax_;
    const bool graphics_;
 
    Real omega_;   // angular frequency in 1/s
    Real sigma_;   // conductivity constant in 1/(Ohm*m)
    Real eps_;     // dielectricity constant in s/(Ohm*m)
    Real mu_;      // permeability constant in Ohm*s/m
    Real j0_;      // current density in A/m^2
    concepts::MaxwellBoundary bc_;
 
    void testMaxwellSolvers(concepts::EddyGeometry2D& geom,
          std::string probStr);
    void testMaxwellSolversDD(concepts::EddyGeometry2D& geom,
            std::string probStr);
 
    void testMaxwellSolver
    (concepts::ModelControl<hp2D::Maxwell2D_H_Base>& control,
     enum hp2D::Maxwell2D_H_Base::solverType type, bool diagPrecond,
     Real& residual, Real& solvetime);
  };
 
} // namespace test
 
# endif // MeshImp2DTest_hh