models/testsuite/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