Eddy2D_geometries.hh

Go to the documentation of this file.

 
#ifndef Eddy2D_geometries_hh
#define Eddy2D_geometries_hh
 
#include <set>
#include "toolbox/inputOutput.hh"
#include "toolbox/inputParam.hh"
#include "toolbox/set.hh"
#include "formula/piecewiseFormula.hh"
#include "geometry/mesh.hh"
#include "geometry/formula.hh"
#include "geometry/circle.hh"
#include "maxwellConstants.hh"
 
namespace concepts {
 
  // ******************************************************** EddyGeometry2D **
 
  class EddyGeometry2D : public OutputOperator {
  public:
    EddyGeometry2D(const std::string meshAbbr, const Real sigma0 = 0.0) 
      : sigma_(sigma0), meshAbbr_(meshAbbr) {}
    virtual ~EddyGeometry2D() {}
    virtual Mesh2& mesh() = 0;
    Set<uint> outerBoundary() const { return outerBoundary_; }
    Set<uint> innerBoundary() const { return innerBoundary_; }
    virtual const PiecewiseFormulaBase<Real>* J0x() const = 0;
    virtual const PiecewiseFormulaBase<Real>* J0y() const = 0;
    virtual const PiecewiseFormulaBase<Real2d>* J0() const = 0;
    PiecewiseConstFormula<Real>& sigma() { return sigma_;}
    const PiecewiseConstFormula<Real>& sigma() const { return sigma_;}
    void meshAbbreviation(const char* meshAbbr) { meshAbbr_ = meshAbbr; }
    const std::string meshAbbreviation() const { return meshAbbr_; }
    virtual void dataout(concepts::InOutParameters& output);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
    Set<uint> outerBoundary_, innerBoundary_;
    PiecewiseConstFormula<Real> sigma_;
    std::string meshAbbr_;
  };
 
  // ********************************************** EddyGeometry2DRotational **
 
  class EddyGeometry2DRotational : public EddyGeometry2D {
  public:
    EddyGeometry2DRotational(const Array<Real> r,
           const std::set<uint> conductive, 
           const std::set<uint> coil,
           const Real j0 = 1.0,
           const Real sigma = SIGMA_CU,
           const Real sigma0 = 0.0);
    EddyGeometry2DRotational(const Array<Real> r, const Real ri,
           const std::set<uint> conductive, 
           const std::set<uint> coil,
           const Real j0 = 1.0,
           const Real sigma = SIGMA_CU,
           const Real sigma0 = 0.0);
    virtual ~EddyGeometry2DRotational();
    virtual Mesh2& mesh() { return *msh_; }
    virtual const PiecewiseFormulaBase<Real>* J0x() const { return &J0x_; }
    virtual const PiecewiseFormulaBase<Real>* J0y() const { return &J0y_; }
    virtual const PiecewiseFormulaBase<Real2d>* J0() const { return 0; }
    virtual void dataout(concepts::InOutParameters& output);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    std::unique_ptr<Mesh2> msh_;
    concepts::PiecewiseFormula<Real> J0x_, J0y_;
    const Real j0_;
    Array<Real> r_;
    const std::set<uint> conductive_;
    const std::set<uint> coil_;
  };
  
  // ********************************************** EddyGeometry2DRectImport **
 
  class EddyGeometry2DRectImport : public EddyGeometry2D {
  public:
    EddyGeometry2DRectImport(const std::string coord, const std::string elms,
           const std::string attributes,
           const uint idxStart = 1,
           const std::string meshAbbr = "",
           const Real j0 = 1.0,
           const Real sigma = SIGMA_CU,
           const Real sigma0 = 0.0);
    EddyGeometry2DRectImport(const std::string path,
           const std::string meshAbbr,
           const Real j0 = 1.0,
           const Real sigma = SIGMA_CU,
           const Real sigma0 = 0.0);
    virtual Mesh2& mesh() { return *msh_; }
    virtual const PiecewiseFormulaBase<Real>* J0x() const { return 0; }
    virtual const PiecewiseFormulaBase<Real>* J0y() const { return 0; }
    virtual const PiecewiseFormulaBase<Real2d>* J0() const { return &J0_; }
    virtual void dataout(concepts::InOutParameters& output);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    std::unique_ptr<Mesh2> msh_;
    concepts::PiecewiseConstFormula<Real2d> J0_;
    const Real j0_;
    const std::string coord_, elms_, attr_;
  };
 
  // ******************************************** EddyGeometry2DRotateImport **
 
  class EddyGeometry2DRotateImport : public EddyGeometry2D {
  public:
    EddyGeometry2DRotateImport(const std::string coord, const std::string elms,
             const std::string attributes,
             const std::string edgRadia,
             const uint idxStart = 1,
             const std::string meshAbbr = "",
             const Real j0 = 1.0,
             const Real sigma = SIGMA_CU,
             const Real sigma0 = 0.0);
    virtual Mesh2& mesh() { return *msh_; }
    virtual const PiecewiseFormulaBase<Real>* J0x() const { return &J0x_; }
    virtual const PiecewiseFormulaBase<Real>* J0y() const { return &J0y_; }
    virtual const PiecewiseFormulaBase<Real2d>* J0() const { return &J0_; }
    virtual void dataout(concepts::InOutParameters& output);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    std::unique_ptr<Mesh2> msh_;
    concepts::PiecewiseFormula<Real> J0x_, J0y_;
    concepts::PiecewiseFormula<Real2d> J0_;
    const Real j0_;
    const std::string coord_, elms_, attr_, edgRadia_;
  };
 
  // ************************************************* InputEddy2DGeometries **
 
  class InputEddy2DGeometries : public concepts::InputParameter {
  public:
    InputEddy2DGeometries(concepts::InOutParameters& input,
        const Real sigma = SIGMA_CU, const Real sigma0 = 0);
    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);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  };
 
} // namespace concepts
 
#endif // Eddy2D_geometries_hh