PEC_Insulator.H

Go to the documentation of this file.

#ifndef PEC_INSULATOR_H_
#define PEC_INSULATOR_H_
 
#include "Utils/WarpXAlgorithmSelection.H"
 
#include <AMReX_Array.H>
#include <AMReX_Geometry.H>
#include <AMReX_Vector.H>
 
#include <AMReX_BaseFwd.H>
 
#include <array>
#include <memory>
 
class PEC_Insulator
{
public:
 
    PEC_Insulator();

    void ApplyPEC_InsulatortoEfield (std::array<amrex::MultiFab*, 3> Efield,
                                     amrex::Array<FieldBoundaryType,AMREX_SPACEDIM> const & field_boundary_lo,
                                     amrex::Array<FieldBoundaryType,AMREX_SPACEDIM> const & field_boundary_hi,
                                     amrex::IntVect const & ng_fieldgather, amrex::Geometry const & geom,
                                     int lev, PatchType patch_type, amrex::Vector<amrex::IntVect> const & ref_ratios,
                                     amrex::Real time,
                                     bool split_pml_field = false);
    void ApplyPEC_InsulatortoBfield (std::array<amrex::MultiFab*, 3> Bfield,
                                     amrex::Array<FieldBoundaryType,AMREX_SPACEDIM> const & field_boundary_lo,
                                     amrex::Array<FieldBoundaryType,AMREX_SPACEDIM> const & field_boundary_hi,
                                     amrex::IntVect const & ng_fieldgather, amrex::Geometry const & geom,
                                     int lev, PatchType patch_type, amrex::Vector<amrex::IntVect> const & ref_ratios,
                                     amrex::Real time);

    void
    ApplyPEC_InsulatortoField (std::array<amrex::MultiFab*, 3> field,
                               amrex::Array<FieldBoundaryType,AMREX_SPACEDIM> const & field_boundary_lo,
                               amrex::Array<FieldBoundaryType,AMREX_SPACEDIM> const & field_boundary_hi,
                               amrex::IntVect const & ng_fieldgather, amrex::Geometry const & geom,
                               int lev, PatchType patch_type, amrex::Vector<amrex::IntVect> const & ref_ratios,
                               amrex::Real time,
                               bool split_pml_field,
                               bool E_like
#ifndef WARPX_DIM_1D_Z
                               , bool set_F_x_lo, bool set_F_x_hi,
                               std::unique_ptr<amrex::Parser> const & a_Fy_x_lo, std::unique_ptr<amrex::Parser> const & a_Fz_x_lo,
                               std::unique_ptr<amrex::Parser> const & a_Fy_x_hi, std::unique_ptr<amrex::Parser> const & a_Fz_x_hi
#endif
#if defined(WARPX_DIM_3D)
                               , bool set_F_y_lo, bool set_F_y_hi,
                               std::unique_ptr<amrex::Parser> const & a_Fx_y_lo, std::unique_ptr<amrex::Parser> const & a_Fz_y_lo,
                               std::unique_ptr<amrex::Parser> const & a_Fx_y_hi, std::unique_ptr<amrex::Parser> const & a_Fz_y_hi
#endif
#if defined(WARPX_ZINDEX)
                               , bool set_F_z_lo, bool set_F_z_hi,
                               std::unique_ptr<amrex::Parser> const & a_Fx_z_lo, std::unique_ptr<amrex::Parser> const & a_Fy_z_lo,
                               std::unique_ptr<amrex::Parser> const & a_Fx_z_hi, std::unique_ptr<amrex::Parser> const & a_Fy_z_hi
#endif
                               );
 
private:
 
    /* \brief Reads in the parsers for the tangential fields, returning whether
     *        the input parameter was specified.
     * \param[in]  pp_insulator ParmParse instance
     * \param[out] parser the parser generated from the input
     * \param[in]  input_name the name of the input parameter
     * \param[in]  coord1 the first coordinate in the plane
     * \param[in]  coord2 the second coordinate in the plane
     */
    bool ReadTangentialFieldParser (amrex::ParmParse const & pp_insulator,
                                    std::unique_ptr<amrex::Parser> & parser,
                                    std::string const & input_name,
                                    std::string const & coord1,
                                    std::string const & coord2);
 
    std::vector<std::unique_ptr<amrex::Parser>> m_insulator_area_lo;
    std::vector<std::unique_ptr<amrex::Parser>> m_insulator_area_hi;
 
#ifndef WARPX_DIM_1D_Z
    bool m_set_B_x_lo = false, m_set_B_x_hi = false;
    std::unique_ptr<amrex::Parser> m_By_x_lo, m_Bz_x_lo;
    std::unique_ptr<amrex::Parser> m_By_x_hi, m_Bz_x_hi;
#endif
#if defined(WARPX_DIM_3D)
    bool m_set_B_y_lo = false, m_set_B_y_hi = false;
    std::unique_ptr<amrex::Parser> m_Bx_y_lo, m_Bz_y_lo;
    std::unique_ptr<amrex::Parser> m_Bx_y_hi, m_Bz_y_hi;
#endif
    bool m_set_B_z_lo = false, m_set_B_z_hi = false;
    std::unique_ptr<amrex::Parser> m_Bx_z_lo, m_By_z_lo;
    std::unique_ptr<amrex::Parser> m_Bx_z_hi, m_By_z_hi;
 
 
#ifndef WARPX_DIM_1D_Z
    bool m_set_E_x_lo = false, m_set_E_x_hi = false;
    std::unique_ptr<amrex::Parser> m_Ey_x_lo, m_Ez_x_lo;
    std::unique_ptr<amrex::Parser> m_Ey_x_hi, m_Ez_x_hi;
#endif
#if defined(WARPX_DIM_3D)
    bool m_set_E_y_lo = false, m_set_E_y_hi = false;
    std::unique_ptr<amrex::Parser> m_Ex_y_lo, m_Ez_y_lo;
    std::unique_ptr<amrex::Parser> m_Ex_y_hi, m_Ez_y_hi;
#endif
    bool m_set_E_z_lo = false, m_set_E_z_hi = false;
    std::unique_ptr<amrex::Parser> m_Ex_z_lo, m_Ey_z_lo;
    std::unique_ptr<amrex::Parser> m_Ex_z_hi, m_Ey_z_hi;
 
 
};
#endif // PEC_INSULATOR_H_