1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
|
#include "WarpX.H"
#include "BoundaryConditions/PML.H"
#include "FieldSolver/FiniteDifferenceSolver/FiniteDifferenceSolver.H"
#include "Evolve/WarpXDtType.H"
#include "WarpX_PEC.H"
#include <AMReX.H>
#include <AMReX_Geometry.H>
#include <AMReX_IntVect.H>
#include <AMReX_REAL.H>
#include <AMReX_Vector.H>
#include <AMReX_Print.H>
#include <algorithm>
#include <array>
#include <memory>
using namespace amrex::literals;
using namespace amrex;
void WarpX::ApplyEfieldBoundary(const int lev, PatchType patch_type)
{
if (PEC::isAnyBoundaryPEC()) {
if (patch_type == PatchType::fine) {
PEC::ApplyPECtoEfield( { get_pointer_Efield_fp(lev, 0),
get_pointer_Efield_fp(lev, 1),
get_pointer_Efield_fp(lev, 2) }, lev, patch_type);
if (WarpX::isAnyBoundaryPML()) {
// apply pec on split E-fields in PML region
const bool split_pml_field = true;
PEC::ApplyPECtoEfield( pml[lev]->GetE_fp(), lev, patch_type, split_pml_field);
}
} else {
PEC::ApplyPECtoEfield( { get_pointer_Efield_cp(lev, 0),
get_pointer_Efield_cp(lev, 1),
get_pointer_Efield_cp(lev, 2) }, lev, patch_type);
if (WarpX::isAnyBoundaryPML()) {
// apply pec on split E-fields in PML region
const bool split_pml_field = true;
PEC::ApplyPECtoEfield( pml[lev]->GetE_cp(), lev, patch_type, split_pml_field);
}
}
}
}
void WarpX::ApplyBfieldBoundary (const int lev, PatchType patch_type, DtType a_dt_type)
{
if (PEC::isAnyBoundaryPEC()) {
if (patch_type == PatchType::fine) {
PEC::ApplyPECtoBfield( { get_pointer_Bfield_fp(lev, 0),
get_pointer_Bfield_fp(lev, 1),
get_pointer_Bfield_fp(lev, 2) }, lev, patch_type);
} else {
PEC::ApplyPECtoBfield( { get_pointer_Bfield_cp(lev, 0),
get_pointer_Bfield_cp(lev, 1),
get_pointer_Bfield_cp(lev, 2) }, lev, patch_type);
}
}
// Silver-Mueller boundaries are only applied on the first half-push of B
// This is because the formula used for Silver-Mueller assumes that
// E and B are staggered in time, which is only true after the first half-push
if (lev == 0) {
if (a_dt_type == DtType::FirstHalf) {
bool applySilverMueller = false;
for (int idim = 0; idim < AMREX_SPACEDIM; ++idim) {
if ( (WarpX::field_boundary_lo[idim] == FieldBoundaryType::Absorbing_SilverMueller) ||
(WarpX::field_boundary_hi[idim] == FieldBoundaryType::Absorbing_SilverMueller) ) {
applySilverMueller = true;
}
}
if(applySilverMueller) m_fdtd_solver_fp[0]->ApplySilverMuellerBoundary(
Efield_fp[lev], Bfield_fp[lev],
Geom(lev).Domain(), dt[lev],
WarpX::field_boundary_lo,
WarpX::field_boundary_hi);
}
}
}
|
