Functions
AMREX_GPU_DEVICE AMREX_FORCE_INLINE bool	is_boundary_PEC (amrex::GpuArray< int, 3 > const &fboundary, int dir)
	Determines if the field boundary condition stored in fboundary is PEC in direction, dir, is PEC. More...

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	SetEfieldOnPEC (const int icomp, const amrex::IntVect &dom_lo, const amrex::IntVect &dom_hi, const amrex::IntVect &ijk_vec, const int n, amrex::Array4< amrex::Real > const &Efield, const amrex::IntVect &is_nodal, amrex::GpuArray< int, 3 > const &fbndry_lo, amrex::GpuArray< int, 3 > const &fbndry_hi)
	Sets the electric field value tangential to the PEC boundary to zero. The tangential Efield components in the guard cells outside the domain boundary are set equal and opposite to the field in the valid cells at their mirrored locations. The normal Efield components in the guard cells are set equal to the field in the valid cells at their mirrored locations. The number or depth of guard cells updated is equal to the shape factor of particles in each dimension. For corner cells with mixed boundaries, the mirror location could be outside valid region, while still ensuring PEC condition is maintained across the PEC boundary, and the necessary sign change is accounted for depending on if the component, icomp, is tangential or normal to the PEC boundary. More...

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	SetBfieldOnPEC (const int icomp, const amrex::IntVect &dom_lo, const amrex::IntVect &dom_hi, const amrex::IntVect &ijk_vec, const int n, amrex::Array4< amrex::Real > const &Bfield, const amrex::IntVect &is_nodal, amrex::GpuArray< int, 3 > const &fbndry_lo, amrex::GpuArray< int, 3 > const &fbndry_hi)
	Sets the magnetic field value normal to the PEC boundary to zero. The tangential (and normal) field value of the guard cells outside the domain boundary are set equal (and opposite) to the respective field components in the valid cells at their mirrored locations. The number or depth of guard cells updated is equal to the shape factor of particles in each dimension. More...

bool	isAnyBoundaryPEC ()

void	ApplyPECtoEfield (std::array< amrex::MultiFab *, 3 > Efield, const int lev, PatchType patch_type, const bool split_pml_field=false)
	Sets the tangential electric field at the PEC boundary to zero. The guard cell values are set equal and opposite to the valid cell field value at the respective mirror locations. More...

void	ApplyPECtoBfield (std::array< amrex::MultiFab *, 3 > Bfield, const int lev, PatchType patch_type)
	Sets the normal component of the magnetic field at the PEC boundary to zero. The guard cell values are set equal and opposite to the valid cell field value at the respective mirror locations. More...

Function Documentation

◆ ApplyPECtoBfield()

void PEC::ApplyPECtoBfield	(	std::array< amrex::MultiFab *, 3 >	Bfield,
		const int	lev,
		PatchType	patch_type
	)

Sets the normal component of the magnetic field at the PEC boundary to zero. The guard cell values are set equal and opposite to the valid cell field value at the respective mirror locations.

Parameters

[in,out]	Bfield	Boundary values of normal Bfield are set to zero.
[in]	lev	level of the Multifab
[in]	patch_type	coarse or fine

◆ ApplyPECtoEfield()

void PEC::ApplyPECtoEfield	(	std::array< amrex::MultiFab *, 3 >	Efield,
		const int	lev,
		PatchType	patch_type,
		const bool	split_pml_field = `false`
	)

Sets the tangential electric field at the PEC boundary to zero. The guard cell values are set equal and opposite to the valid cell field value at the respective mirror locations.

Parameters

[in,out]	Efield	Boundary values of tangential Efield are set to zero
[in]	lev	level of the Multifab
[in]	patch_type	coarse or fine
[in]	split_pml_field	whether pml the multifab is the regular Efield or split pml field

◆ is_boundary_PEC()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE bool PEC::is_boundary_PEC	(	amrex::GpuArray< int, 3 > const &	fboundary,
		int	dir
	)

Determines if the field boundary condition stored in fboundary is PEC in direction, dir, is PEC.

Parameters

[in]	fboundary	Value containing boundary type
[in]	dir	direction

Returns: 1 if the boundary type is PEC else 0

◆ isAnyBoundaryPEC()

bool PEC::isAnyBoundaryPEC ( )

Returns 1 if any domain boundary is set to PEC, else returns 0.

◆ SetBfieldOnPEC()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void PEC::SetBfieldOnPEC	(	const int	icomp,
		const amrex::IntVect &	dom_lo,
		const amrex::IntVect &	dom_hi,
		const amrex::IntVect &	ijk_vec,
		const int	n,
		amrex::Array4< amrex::Real > const &	Bfield,
		const amrex::IntVect &	is_nodal,
		amrex::GpuArray< int, 3 > const &	fbndry_lo,
		amrex::GpuArray< int, 3 > const &	fbndry_hi
	)

Sets the magnetic field value normal to the PEC boundary to zero. The tangential (and normal) field value of the guard cells outside the domain boundary are set equal (and opposite) to the respective field components in the valid cells at their mirrored locations. The number or depth of guard cells updated is equal to the shape factor of particles in each dimension.

For 3D : x component is tangential to the y-boundary and z-boundary y component is tangential to the x-boundary and z-boundary z component is tangential to the x-boundary and y-boundary x component is normal to the x-boundary y component is normal to the y-boundary z component is normal to the z-boundary where, x-boundary is the yz-plane at x=xmin and x=xmax y-boundary is the xz-plane at y=ymin and y=ymax z-boundary is the xy-plane at z=zmin and z=zmax

For 2D : WarpX uses X-Z as the two dimensions x component is tangential to the z-boundary y component is tangential to the x-boundary and z-boundary z component is tangential to the x-boundary x component is normal to the x-boundary y component is not normal to any boundary (Only xz dimensions in 2D) z component is normal to the z-boundary where, x-boundary is along the line z at x=xmin and x=xmax z-boundary is along the line x at z=zmin and z=zmax

For 1D : WarpX uses Z as the only dimension x component is tangential to the z-boundary y component is tangential to the z-boundary z component is not tangential to the z-boundary x component is not normal to any boundary (Only z dimension in 1D) y component is not normal to any boundary (Only z dimension in 1D) z component is normal to the z-boundary where, z-boundary is a point at z=zmin and z=zmax

For RZ : WarpX uses R-Z as the two dimensions r component is tangential to the z-boundary theta_component is tangential to the r-boundary and z-boundary z component is tangential to the r-boundary r component is normal to the r-boundary theta_component is not normal to any boundary (on RZ dimensions are modeled) z component is normal to the z-boundary where, r-boundary is along the line z at r=rmin and r=rmax z-boundary is along the line r at z=zmin and z=zmax

Parameters

[in]	icomp	component of the Bfield being updated (0=x, 1=y, 2=z in Cartesian) (0=r, 1=theta, 2=z in RZ)
[in]	dom_lo	index value of the lower domain boundary (cell-centered)
[in]	dom_hi	index value of the higher domain boundary (cell-centered)
[in]	ijk_vec	indices along the x(i), y(j), z(k) of Efield Array4
[in]	n	index of the MultiFab component being updated
[in]	Bfield	field data to be updated if (ijk) is at the boundary or a guard cell
[in]	is_nodal	staggering of the field data being updated.
[in]	fbndry_lo	Field boundary type at the lower boundaries
[in]	fbndry_hi	Field boundary type at the upper boundaries

◆ SetEfieldOnPEC()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void PEC::SetEfieldOnPEC	(	const int	icomp,
		const amrex::IntVect &	dom_lo,
		const amrex::IntVect &	dom_hi,
		const amrex::IntVect &	ijk_vec,
		const int	n,
		amrex::Array4< amrex::Real > const &	Efield,
		const amrex::IntVect &	is_nodal,
		amrex::GpuArray< int, 3 > const &	fbndry_lo,
		amrex::GpuArray< int, 3 > const &	fbndry_hi
	)

Sets the electric field value tangential to the PEC boundary to zero. The tangential Efield components in the guard cells outside the domain boundary are set equal and opposite to the field in the valid cells at their mirrored locations. The normal Efield components in the guard cells are set equal to the field in the valid cells at their mirrored locations. The number or depth of guard cells updated is equal to the shape factor of particles in each dimension. For corner cells with mixed boundaries, the mirror location could be outside valid region, while still ensuring PEC condition is maintained across the PEC boundary, and the necessary sign change is accounted for depending on if the component, icomp, is tangential or normal to the PEC boundary.

For 3D : x component is tangential to the y-boundary and z-boundary y component is tangential to the x-boundary and z-boundary z component is tangential to the x-boundary and y-boundary x component is normal to the x-boundary y component is normal to the y-boundary z component is normal to the z-boundary where, x-boundary is the yz-plane at x=xmin and x=xmax y-boundary is the xz-plane at y=ymin and y=ymax z-boundary is the xy-plane at z=zmin and z=zmax

For 2D : WarpX uses X-Z as the two dimensions x component is tangential to the z-boundary y component is tangential to the x-boundary and z-boundary z component is tangential to the x-boundary x component is normal to the x-boundary y component is not normal to any boundary (Only xz dimensions in 2D) z component is normal to the z-boundary where, x-boundary is along the line z at x=xmin and x=xmax z-boundary is along the line x at z=zmin and z=zmax

For 1D : WarpX uses Z as the only dimension x component is tangential to the z-boundary y component is tangential to the z-boundary z component is not tangential to the z-boundary x component is not normal to any boundary (Only z dimension in 1D) y component is not normal to any boundary (Only z dimension in 1D) z component is normal to the z-boundary where, z-boundary is a point at z=zmin and z=zmax

For RZ : WarpX uses R-Z as the two dimensions r component is tangential to the z-boundary theta_component is tangential to the r-boundary and z-boundary z component is tangential to the r-boundary r component is normal to the r-boundary theta_component is not normal to any boundary (on RZ dimensions are modeled) z component is normal to the z-boundary where, r-boundary is along the line z at r=rmin and r=rmax z-boundary is along the line r at z=zmin and z=zmax

Parameters

[in]	icomp	component of the Efield being updated (0=x, 1=y, 2=z in Cartesian) (0=r, 1=theta, 2=z in RZ)
[in]	dom_lo	index value of the lower domain boundary (cell-centered)
[in]	dom_hi	index value of the higher domain boundary (cell-centered)
[in]	ijk_vec	indices along the x(i), y(j), z(k) of Efield Array4
[in]	n	index of the MultiFab component being updated
[in]	Efield	field data to be updated if (ijk) is at the boundary or a guard cell
[in]	is_nodal	staggering of the field data being updated.
[in]	fbndry_lo	Field boundary type at the lower boundaries
[in]	fbndry_hi	Field boundary type at the upper boundaries

Functions

Function Documentation

◆ ApplyPECtoBfield()

◆ ApplyPECtoEfield()

◆ is_boundary_PEC()

◆ isAnyBoundaryPEC()

◆ SetBfieldOnPEC()

◆ SetEfieldOnPEC()