14 #ifdef WARPX_USE_PSATD
20 #include <AMReX_Config.H>
81 const amrex::Box& regular_domain,
const amrex::Real v_sigma_sb)
93 const
amrex::FabInfo& ,
int ) const final
117 const amrex::Box& regular_domain,
const amrex::Real v_sigma_sb);
134 amrex::Real
dt,
int nox_fft,
int noy_fft,
int noz_fft,
short grid_type,
135 int do_moving_window,
int pml_has_particles,
int do_pml_in_domain,
136 const int psatd_solution_type,
const int J_in_time,
const int rho_in_time,
137 const bool do_pml_dive_cleaning,
const bool do_pml_divb_cleaning,
140 int max_guard_EB, amrex::Real v_sigma_sb,
146 std::array<amrex::MultiFab*,3>
GetE_fp ();
147 std::array<amrex::MultiFab*,3>
GetB_fp ();
148 std::array<amrex::MultiFab*,3>
Getj_fp ();
149 std::array<amrex::MultiFab*,3>
GetE_cp ();
150 std::array<amrex::MultiFab*,3>
GetB_cp ();
151 std::array<amrex::MultiFab*,3>
Getj_cp ();
155 std::array<amrex::MultiFab*,3> GetH_fp ();
156 std::array<amrex::MultiFab*,3> GetH_cp ();
185 #ifdef WARPX_USE_PSATD
189 void CopyJtoPMLs (
const std::array<amrex::MultiFab*,3>& j_fp,
190 const std::array<amrex::MultiFab*,3>& j_cp);
192 void Exchange (
const std::array<amrex::MultiFab*,3>& mf_pml,
193 const std::array<amrex::MultiFab*,3>& mf,
195 const int do_pml_in_domain);
198 const std::array<amrex::MultiFab*,3>& jp);
204 void ExchangeH (
const std::array<amrex::MultiFab*,3>& H_fp,
205 const std::array<amrex::MultiFab*,3>& H_cp,
int do_pml_in_domain);
207 const std::array<amrex::MultiFab*,3>& Hp,
int do_pml_in_domain);
213 void FillBoundaryE ();
214 void FillBoundaryB ();
215 void FillBoundaryF ();
216 void FillBoundaryG ();
217 void FillBoundaryE (
PatchType patch_type);
218 void FillBoundaryB (
PatchType patch_type);
219 void FillBoundaryF (
PatchType patch_type);
221 void FillBoundaryH ();
222 void FillBoundaryH (
PatchType patch_type);
224 void FillBoundaryG (
PatchType patch_type);
226 bool ok ()
const {
return m_ok; }
228 void CheckPoint (
const std::string& dir)
const;
229 void Restart (
const std::string& dir);
247 std::array<std::unique_ptr<amrex::MultiFab>,3>
pml_E_fp;
248 std::array<std::unique_ptr<amrex::MultiFab>,3>
pml_B_fp;
249 std::array<std::unique_ptr<amrex::MultiFab>,3>
pml_j_fp;
253 std::array<std::unique_ptr<amrex::MultiFab>,3>
pml_E_cp;
254 std::array<std::unique_ptr<amrex::MultiFab>,3>
pml_B_cp;
255 std::array<std::unique_ptr<amrex::MultiFab>,3>
pml_j_cp;
262 std::array<std::unique_ptr<amrex::MultiFab>,3> pml_H_fp;
263 std::array<std::unique_ptr<amrex::MultiFab>,3> pml_H_cp;
280 #ifdef WARPX_USE_PSATD
289 return *pml_field_factory;
302 int do_pml_in_domain,
315 int do_pml_in_domain,
322 #ifdef WARPX_USE_PSATD
325 std::array<std::unique_ptr<amrex::MultiFab>,3>& pml_E,
326 std::array<std::unique_ptr<amrex::MultiFab>,3>& pml_B,
327 std::unique_ptr<amrex::MultiFab>& pml_F,
328 std::unique_ptr<amrex::MultiFab>& pml_G,
void PushPMLPSATDSinglePatch(const int lev, SpectralSolver &solver, std::array< std::unique_ptr< amrex::MultiFab >, 3 > &pml_E, std::array< std::unique_ptr< amrex::MultiFab >, 3 > &pml_B, std::unique_ptr< amrex::MultiFab > &pml_F, std::unique_ptr< amrex::MultiFab > &pml_G, const amrex::IntVect &fill_guards)
Definition: PML.cpp:1622
PatchType
Definition: WarpX.H:73
amrex::Real dt_E
Definition: PML.H:122
void ComputePMLFactorsE(const amrex::Real *dx, amrex::Real dt)
Definition: PML.cpp:533
void ComputePMLFactorsB(const amrex::Real *dx, amrex::Real dt)
Definition: PML.cpp:517
amrex::Real dt_B
Definition: PML.H:121
MultiSigmaBox(const amrex::BoxArray &ba, const amrex::DistributionMapping &dm, const amrex::BoxArray &grid_ba, const amrex::Real *dx, const amrex::IntVect &ncell, const amrex::IntVect &delta, const amrex::Box ®ular_domain, const amrex::Real v_sigma_sb)
Definition: PML.cpp:508
amrex::MultiFab * Getmu_fp()
Definition: PML.cpp:1181
const MultiSigmaBox & GetMultiSigmaBox_fp() const
Definition: PML.H:179
amrex::MultiFab * Getsigma_fp()
Definition: PML.cpp:1187
const amrex::Geometry * m_geom
Definition: PML.H:244
std::array< std::unique_ptr< amrex::MultiFab >, 3 > pml_edge_lengths
Definition: PML.H:251
std::unique_ptr< amrex::MultiFab > pml_eps_cp
Definition: PML.H:270
void ComputePMLFactors(amrex::Real dt)
Definition: PML.cpp:1084
amrex::MultiFab * GetG_cp()
Definition: PML.cpp:1217
std::unique_ptr< amrex::FabFactory< amrex::FArrayBox > > pml_field_factory
Definition: PML.H:286
std::unique_ptr< amrex::MultiFab > pml_F_fp
Definition: PML.H:258
const amrex::IntVect m_fill_guards_fields
Definition: PML.H:241
std::array< amrex::MultiFab *, 3 > GetE_fp()
Definition: PML.cpp:1097
std::array< amrex::MultiFab *, 3 > Getj_cp()
Definition: PML.cpp:1149
std::array< std::unique_ptr< amrex::MultiFab >, 3 > pml_B_fp
Definition: PML.H:248
amrex::MultiFab * GetF_fp()
Definition: PML.cpp:1162
std::unique_ptr< MultiSigmaBox > sigba_fp
Definition: PML.H:277
const amrex::Geometry * m_cgeom
Definition: PML.H:245
bool ok() const
Definition: PML.H:226
std::array< std::unique_ptr< amrex::MultiFab >, 3 > pml_j_fp
Definition: PML.H:249
amrex::FabFactory< amrex::FArrayBox > const & fieldFactory() const noexcept
Definition: PML.H:288
amrex::MultiFab * Geteps_cp()
Definition: PML.cpp:1193
std::unique_ptr< amrex::MultiFab > pml_mu_fp
Definition: PML.H:268
bool m_divb_cleaning
Definition: PML.H:239
amrex::MultiFab * Geteps_fp()
Definition: PML.cpp:1175
bool m_ok
Definition: PML.H:236
std::array< amrex::MultiFab *, 3 > Getj_fp()
Definition: PML.cpp:1123
std::unique_ptr< amrex::MultiFab > pml_eps_fp
Definition: PML.H:267
PML(const int lev, const amrex::BoxArray &ba, const amrex::DistributionMapping &dm, const amrex::Geometry *geom, const amrex::Geometry *cgeom, int ncell, int delta, amrex::IntVect ref_ratio, amrex::Real dt, int nox_fft, int noy_fft, int noz_fft, short grid_type, int do_moving_window, int pml_has_particles, int do_pml_in_domain, const int psatd_solution_type, const int J_in_time, const int rho_in_time, const bool do_pml_dive_cleaning, const bool do_pml_divb_cleaning, const amrex::IntVect &fill_guards_fields, const amrex::IntVect &fill_guards_current, int max_guard_EB, amrex::Real v_sigma_sb, const amrex::IntVect do_pml_Lo=amrex::IntVect::TheUnitVector(), const amrex::IntVect do_pml_Hi=amrex::IntVect::TheUnitVector())
Definition: PML.cpp:548
std::array< amrex::MultiFab *, 3 > Get_edge_lengths()
Definition: PML.cpp:1155
bool m_dive_cleaning
Definition: PML.H:238
void PushPSATD(const int lev)
Definition: PML.cpp:1612
std::array< std::unique_ptr< amrex::MultiFab >, 3 > pml_B_cp
Definition: PML.H:254
amrex::MultiFab * GetF_cp()
Definition: PML.cpp:1168
amrex::MultiFab * Getsigma_cp()
Definition: PML.cpp:1205
amrex::MultiFab * GetG_fp()
Definition: PML.cpp:1211
const MultiSigmaBox & GetMultiSigmaBox_cp() const
Definition: PML.H:182
std::unique_ptr< amrex::MultiFab > pml_sigma_cp
Definition: PML.H:272
const amrex::IntVect m_fill_guards_current
Definition: PML.H:242
std::array< std::unique_ptr< amrex::MultiFab >, 3 > pml_E_cp
Definition: PML.H:253
std::array< std::unique_ptr< amrex::MultiFab >, 3 > pml_j_cp
Definition: PML.H:255
std::unique_ptr< MultiSigmaBox > sigba_cp
Definition: PML.H:278
std::array< amrex::MultiFab *, 3 > GetE_cp()
Definition: PML.cpp:1129
std::unique_ptr< amrex::MultiFab > pml_mu_cp
Definition: PML.H:271
amrex::MultiFab * Getmu_cp()
Definition: PML.cpp:1199
std::unique_ptr< amrex::MultiFab > pml_F_cp
Definition: PML.H:259
std::array< amrex::MultiFab *, 3 > Get_face_areas()
std::unique_ptr< SpectralSolver > spectral_solver_cp
Definition: PML.H:282
std::array< std::unique_ptr< amrex::MultiFab >, 3 > pml_E_fp
Definition: PML.H:247
void Exchange(const std::array< amrex::MultiFab *, 3 > &mf_pml, const std::array< amrex::MultiFab *, 3 > &mf, const PatchType &patch_type, const int do_pml_in_domain)
Definition: PML.cpp:1222
std::unique_ptr< amrex::MultiFab > pml_sigma_fp
Definition: PML.H:269
std::unique_ptr< SpectralSolver > spectral_solver_fp
Definition: PML.H:281
std::array< amrex::MultiFab *, 3 > GetB_fp()
Definition: PML.cpp:1109
std::array< amrex::MultiFab *, 3 > GetB_cp()
Definition: PML.cpp:1135
void CopyJtoPMLs(const std::array< amrex::MultiFab *, 3 > &j_fp, const std::array< amrex::MultiFab *, 3 > &j_cp)
Definition: PML.cpp:1283
std::unique_ptr< amrex::MultiFab > pml_G_fp
Definition: PML.H:274
std::unique_ptr< amrex::MultiFab > pml_G_cp
Definition: PML.H:275
SigmaBoxFactory(const amrex::BoxArray &grid_ba, const amrex::Real *dx, const amrex::IntVect &ncell, const amrex::IntVect &delta, const amrex::Box ®ular_domain, const amrex::Real v_sigma_sb)
Definition: PML.H:79
virtual SigmaBoxFactory * clone() const final
Definition: PML.H:98
amrex::IntVect m_delta
Definition: PML.H:105
SigmaBoxFactory(SigmaBoxFactory &&) noexcept=default
amrex::IntVect m_ncell
Definition: PML.H:104
amrex::Box m_regdomain
Definition: PML.H:106
SigmaBoxFactory(const SigmaBoxFactory &)=default
virtual void destroy(SigmaBox *fab) const final
Definition: PML.H:95
amrex::Real m_v_sigma_sb
Definition: PML.H:107
virtual ~SigmaBoxFactory()=default
const amrex::BoxArray & m_grids
Definition: PML.H:102
virtual SigmaBox * create(const amrex::Box &box, int, const amrex::FabInfo &, int) const final
Definition: PML.H:92
const amrex::Real * m_dx
Definition: PML.H:103
Top-level class for the electromagnetic spectral solver.
Definition: SpectralSolver.H:33
AMREX_GPU_HOST_DEVICE static constexpr AMREX_FORCE_INLINE IntVect TheUnitVector() noexcept
void FillBoundary(amrex::MultiFab &mf, amrex::IntVect ng, bool do_single_precision_comms, const amrex::Periodicity &period, std::optional< bool > nodal_sync)
Definition: Communication.cpp:65
default
Definition: run_alltests.py:113
int dx
Definition: stencil.py:436
int dt
Definition: stencil.py:440
SigmaVect sigma_star_cumsum_fac
Definition: PML.H:70
SigmaVect sigma_star
Definition: PML.H:65
std::array< Sigma, AMREX_SPACEDIM > SigmaVect
Definition: PML.H:59
SigmaVect sigma_star_fac
Definition: PML.H:69
SigmaVect sigma_cumsum
Definition: PML.H:64
SigmaVect sigma_star_cumsum
Definition: PML.H:66
void ComputePMLFactorsB(const amrex::Real *dx, amrex::Real dt)
Definition: PML.cpp:441
void value_type
Definition: PML.H:61
SigmaVect sigma_cumsum_fac
Definition: PML.H:68
SigmaVect sigma_fac
Definition: PML.H:67
void define_multiple(const amrex::Box &box, const amrex::BoxArray &grids, const amrex::IntVect &ncell, const amrex::Array< amrex::Real, AMREX_SPACEDIM > &fac, const amrex::Real v_sigma)
Definition: PML.cpp:240
void ComputePMLFactorsE(const amrex::Real *dx, amrex::Real dt)
Definition: PML.cpp:475
SigmaVect sigma
Definition: PML.H:63
SigmaBox(const amrex::Box &box, const amrex::BoxArray &grids, const amrex::Real *dx, const amrex::IntVect &ncell, const amrex::IntVect &delta, const amrex::Box ®domain, const amrex::Real v_sigma)
Definition: PML.cpp:147
void define_single(const amrex::Box ®domain, const amrex::IntVect &ncell, const amrex::Array< amrex::Real, AMREX_SPACEDIM > &fac, const amrex::Real v_sigma)
Definition: PML.cpp:197
amrex::Real v_sigma
Definition: PML.H:71
int m_lo
Definition: PML.H:39
int lo() const
Definition: PML.H:37
int hi() const
Definition: PML.H:38
int m_hi
Definition: PML.H:39