Functions
template<typename T_PC >
static void	deposit_charge (typename T_PC::ParIterType &pti, typename T_PC::RealVector const &wp, amrex::Real const charge, int const const ion_lev, amrex::MultiFab rho, amrex::FArrayBox &local_rho, int const particle_shape, std::array< amrex::Real, 3 > const &dx, std::array< amrex::Real, 3 > const &xyzmin, int const n_rz_azimuthal_modes=0, std::optional< amrex::IntVect > num_rho_deposition_guards=std::nullopt, std::optional< int > depos_lev=std::nullopt, std::optional< amrex::IntVect > rel_ref_ratio=std::nullopt, long const offset=0, std::optional< long > np_to_depose=std::nullopt, int const icomp=0, int const nc=1, amrex::Real *const AMREX_RESTRICT cost=nullptr, long const load_balance_costs_update_algo=0, bool const do_device_synchronize=true)

constexpr uint64_t	localIDtoGlobal (int const id, int const cpu)

AMREX_GPU_HOST_DEVICE AMREX_INLINE void	compute_weights_nodal (const amrex::ParticleReal xp, const amrex::ParticleReal yp, const amrex::ParticleReal zp, amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &plo, amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &dxi, int &i, int &j, int &k, amrex::Real W[AMREX_SPACEDIM][2]) noexcept
	Compute weight of each surrounding node in interpolating a nodal field to the given coordinates. More...

AMREX_GPU_HOST_DEVICE AMREX_INLINE amrex::Real	interp_field_nodal (int i, int j, int k, const amrex::Real W[AMREX_SPACEDIM][2], amrex::Array4< const amrex::Real > const &scalar_field) noexcept
	Interpolate nodal field value based on surrounding indices and weights. More...

AMREX_GPU_HOST_DEVICE AMREX_INLINE amrex::Real	doGatherScalarFieldNodal (const amrex::ParticleReal xp, const amrex::ParticleReal yp, const amrex::ParticleReal zp, amrex::Array4< const amrex::Real > const &scalar_field, amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &dxi, amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &lo) noexcept
	Scalar field gather for a single particle. The field has to be defined at the cell nodes (see https://amrex-codes.github.io/amrex/docs_html/Basics.html#id2) More...

AMREX_GPU_HOST_DEVICE AMREX_INLINE amrex::GpuArray< amrex::Real, 3 >	doGatherVectorFieldNodal (const amrex::ParticleReal xp, const amrex::ParticleReal yp, const amrex::ParticleReal zp, amrex::Array4< const amrex::Real > const &vector_field_x, amrex::Array4< const amrex::Real > const &vector_field_y, amrex::Array4< const amrex::Real > const &vector_field_z, amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &dxi, amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &lo) noexcept
	Vector field gather for a single particle. The field has to be defined at the cell nodes (see https://amrex-codes.github.io/amrex/docs_html/Basics.html#id2) More...

template<typename T_PC >
static std::tuple< amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal >	MinAndMaxPositions (T_PC const &pc)

template<typename T_PC , int T_RealSoAWeight>
static std::tuple< amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal >	MeanAndStdPositions (T_PC const &pc)

Function Documentation

◆ compute_weights_nodal()

AMREX_GPU_HOST_DEVICE AMREX_INLINE void ablastr::particles::compute_weights_nodal	(	const amrex::ParticleReal	xp,
		const amrex::ParticleReal	yp,
		const amrex::ParticleReal	zp,
		amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &	plo,
		amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &	dxi,
		int &	i,
		int &	j,
		int &	k,
		amrex::Real	W[AMREX_SPACEDIM][2]
	)

noexcept

Compute weight of each surrounding node in interpolating a nodal field to the given coordinates.

This currently only does linear order.

Parameters

xp,yp,zp	Particle position coordinates
plo	Index lower bounds of domain.
dxi	inverse 3D cell spacing
i,j,k	Variables to store indices of position on grid
W	2D array of weights to store each neighbouring node

◆ deposit_charge()

template<typename T_PC >

static void ablastr::particles::deposit_charge	(	typename T_PC::ParIterType &	pti,
		typename T_PC::RealVector const &	wp,
		amrex::Real const	charge,
		int const *const	ion_lev,
		amrex::MultiFab *	rho,
		amrex::FArrayBox &	local_rho,
		int const	particle_shape,
		std::array< amrex::Real, 3 > const &	dx,
		std::array< amrex::Real, 3 > const &	xyzmin,
		int const	n_rz_azimuthal_modes = `0`,
		std::optional< amrex::IntVect >	num_rho_deposition_guards = `std::nullopt`,
		std::optional< int >	depos_lev = `std::nullopt`,
		std::optional< amrex::IntVect >	rel_ref_ratio = `std::nullopt`,
		long const	offset = `0`,
		std::optional< long >	np_to_depose = `std::nullopt`,
		int const	icomp = `0`,
		int const	nc = `1`,
		amrex::Real *const AMREX_RESTRICT	cost = `nullptr`,
		long const	load_balance_costs_update_algo = `0`,
		bool const	do_device_synchronize = `true`
	)

static

Perform charge deposition for the particles on a tile.

Template Parameters

T_PC	a type of amrex::ParticleContainer

Parameters

pti	an amrex::ParIter pointing to the tile to operate on
wp	vector of the particle weights for those particles.
charge	charge of the particle species
ion_lev	pointer to array of particle ionization level. This is required to have the charge of each macroparticle since q is a scalar. For non-ionizable species, ion_lev is a null pointer.
rho	MultiFab of the charge density
local_rho	temporary FArrayBox for deposition with OpenMP
particle_shape	shape factor in each direction
dx	cell spacing at level lev
xyzmin	lo corner of the current tile in physical coordinates.
n_rz_azimuthal_modes	number of azimuthal modes in use, irrelevant outside RZ geometry (default: 0)
num_rho_deposition_guards	number of ghost cells to use for rho (default: rho.nGrowVect())
depos_lev	the level to deposit the particles to (default: lev)
rel_ref_ratio	mesh refinement ratio between lev and depos_lev (default: 1)
offset	index to start at when looping over particles to depose (default: 0)
np_to_depose	number of particles to depose (default: pti.numParticles())
icomp	component in MultiFab to start depositing to
nc	number of components to deposit
cost	pointer to (load balancing) cost corresponding to box where present particles deposit current. If nullptr, costs are not updated. (default: nullptr)
load_balance_costs_update_algo	selected method for updating load balance costs (default: 0)
do_device_synchronize	call amrex::Gpu::synchronize() for tiny profiler regions (default: true)

◆ doGatherScalarFieldNodal()

AMREX_GPU_HOST_DEVICE AMREX_INLINE amrex::Real ablastr::particles::doGatherScalarFieldNodal	(	const amrex::ParticleReal	xp,
		const amrex::ParticleReal	yp,
		const amrex::ParticleReal	zp,
		amrex::Array4< const amrex::Real > const &	scalar_field,
		amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &	dxi,
		amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &	lo
	)

noexcept

Scalar field gather for a single particle. The field has to be defined at the cell nodes (see https://amrex-codes.github.io/amrex/docs_html/Basics.html#id2)

Parameters

xp,yp,zp	Particle position coordinates
scalar_field	Array4 of the nodal scalar field, either full array or tile.
dxi	inverse 3D cell spacing
lo	Index lower bounds of domain.

◆ doGatherVectorFieldNodal()

AMREX_GPU_HOST_DEVICE AMREX_INLINE amrex::GpuArray<amrex::Real, 3> ablastr::particles::doGatherVectorFieldNodal	(	const amrex::ParticleReal	xp,
		const amrex::ParticleReal	yp,
		const amrex::ParticleReal	zp,
		amrex::Array4< const amrex::Real > const &	vector_field_x,
		amrex::Array4< const amrex::Real > const &	vector_field_y,
		amrex::Array4< const amrex::Real > const &	vector_field_z,
		amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &	dxi,
		amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > const &	lo
	)

noexcept

Vector field gather for a single particle. The field has to be defined at the cell nodes (see https://amrex-codes.github.io/amrex/docs_html/Basics.html#id2)

Parameters

xp,yp,zp	Particle position coordinates
vector_field_x,vector_field_y,vector_field_z	Array4 of nodal scalar fields, either full array or tile.
dxi	inverse 3D cell spacing
lo	Index lower bounds of domain.

◆ interp_field_nodal()

AMREX_GPU_HOST_DEVICE AMREX_INLINE amrex::Real ablastr::particles::interp_field_nodal	(	int	i,
		int	j,
		int	k,
		const amrex::Real	W[AMREX_SPACEDIM][2],
		amrex::Array4< const amrex::Real > const &	scalar_field
	)

noexcept

Interpolate nodal field value based on surrounding indices and weights.

Parameters

i,j,k	Indices of position on grid
W	2D array of weights for each neighbouring node
scalar_field	Array4 of the nodal scalar field, either full array or tile.

◆ localIDtoGlobal()

constexpr uint64_t ablastr::particles::localIDtoGlobal	(	int const	id,
		int const	cpu
	)

constexpr

A helper function to derive a globally unique particle ID

Parameters

[in]	id	AMReX particle ID (on local cpu/rank), AoS .id
[in]	cpu	AMReX particle CPU (rank) at creation of the particle, AoS .cpu

Returns: global particle ID that is unique and permanent in the whole simulation

◆ MeanAndStdPositions()

template<typename T_PC , int T_RealSoAWeight>

static std::tuple< amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal> ablastr::particles::MeanAndStdPositions ( T_PC const & pc )

static

Compute the mean and std of the particle position in each dimension

Template Parameters

T_PC	a type of amrex::ParticleContainer
T_RealSoAWeight

Parameters

pc	the particle container to operate on

Returns: x_mean, x_std, y_mean, y_std, z_mean, z_std

◆ MinAndMaxPositions()

template<typename T_PC >

static std::tuple< amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal, amrex::ParticleReal> ablastr::particles::MinAndMaxPositions ( T_PC const & pc )

static

Compute the min and max of the particle position in each dimension

Template Parameters

T_PC	a type of amrex::ParticleContainer

Parameters

pc	the particle container to operate on

Returns: x_min, y_min, z_min, x_max, y_max, z_max

Functions

Function Documentation

◆ compute_weights_nodal()

◆ deposit_charge()

◆ doGatherScalarFieldNodal()

◆ doGatherVectorFieldNodal()

◆ interp_field_nodal()

◆ localIDtoGlobal()

◆ MeanAndStdPositions()

◆ MinAndMaxPositions()