WarpX
Loading...
Searching...
No Matches
ParticleUtils Namespace Reference

Typedefs

using ParticleType = typename WarpXParticleContainer::ParticleType
 
using ParticleBins = DenseBins<ParticleTileDataType>
 
using ParticleTileType = typename WarpXParticleContainer::ParticleTileType
 
using ParticleTileDataType = typename ParticleTileType::ParticleTileDataType
 

Functions

amrex::DenseBins< ParticleTileDataTypefindParticlesInEachCell (amrex::Geometry const &geom_lev, amrex::MFIter const &mfi, WarpXParticleContainer::ParticleTileType &ptile)
 Find the particles and count the particles that are in each cell. More specifically this function returns an amrex::DenseBins object containing an offset array and a permutation array which can be used to loop over all the cells in a tile and apply an algorithm to particles of a given species present in each cell. Note that this does not rearrange particle arrays.
 
AMREX_GPU_HOST_DEVICE AMREX_INLINE void getCollisionEnergy (amrex::ParticleReal const u2, double const m, double const M, double &gamma, double &energy)
 Return (relativistic) collision energy assuming the target (with mass M) is stationary and the projectile is approaching with the the given speed and mass m. Note the use of double since this calculation is prone to error with single precision.
 
AMREX_GPU_HOST_DEVICE AMREX_INLINE void doLorentzTransform (amrex::ParticleReal &ux, amrex::ParticleReal &uy, amrex::ParticleReal &uz, amrex::ParticleReal const Vx, amrex::ParticleReal const Vy, amrex::ParticleReal const Vz)
 Perform a Lorentz transformation of the given velocity to a frame moving with velocity (Vx, Vy, Vz) relative to the present one.
 
AMREX_GPU_HOST_DEVICE AMREX_INLINE void doLorentzTransformWithU (amrex::ParticleReal &ux, amrex::ParticleReal &uy, amrex::ParticleReal &uz, amrex::ParticleReal const Ux, amrex::ParticleReal const Uy, amrex::ParticleReal const Uz)
 Perform a Lorentz transformation of the given velocity to a frame moving with gamma*velocity (Ux, Uy, Uz) relative to the present one.
 
AMREX_GPU_HOST_DEVICE AMREX_INLINE void doLorentzTransformWithP (amrex::ParticleReal &px, amrex::ParticleReal &py, amrex::ParticleReal &pz, amrex::ParticleReal mass, amrex::ParticleReal const Ux, amrex::ParticleReal const Uy, amrex::ParticleReal const Uz)
 Perform a Lorentz transformation of the given momentum to a frame moving with gamma*velocity (Ux, Uy, Uz) relative to the present one.
 
AMREX_GPU_HOST_DEVICE AMREX_INLINE void getRandomVector (amrex::ParticleReal &x, amrex::ParticleReal &y, amrex::ParticleReal &z, amrex::RandomEngine const &engine)
 Generate random unit vector in 3 dimensions https://mathworld.wolfram.com/SpherePointPicking.html.
 
AMREX_GPU_HOST_DEVICE AMREX_INLINE void RandomizeVelocity (amrex::ParticleReal &ux, amrex::ParticleReal &uy, amrex::ParticleReal &uz, const amrex::ParticleReal vp, amrex::RandomEngine const &engine)
 Function to perform scattering of a particle that results in a random velocity vector with given magnitude. This is used in isotropic collision events.
 
AMREX_GPU_HOST_DEVICE AMREX_INLINE bool containsInclusive (amrex::RealBox const &tilebox, amrex::XDim3 const point)
 
template<typename T_index>
AMREX_GPU_HOST_DEVICE AMREX_INLINE bool ModifyEnergyPairwise (amrex::ParticleReal *const AMREX_RESTRICT uxp, amrex::ParticleReal *const AMREX_RESTRICT uyp, amrex::ParticleReal *const AMREX_RESTRICT uzp, amrex::ParticleReal *const AMREX_RESTRICT w, T_index const *const AMREX_RESTRICT indices, T_index const cell_start, T_index const cell_stop, amrex::ParticleReal const mass, amrex::ParticleReal const energy_fraction, amrex::ParticleReal const energy_fraction_max, amrex::ParticleReal &deltaE)
 

Typedef Documentation

◆ ParticleBins

◆ ParticleTileDataType

using ParticleUtils::ParticleTileDataType = typename ParticleTileType::ParticleTileDataType

◆ ParticleTileType

using ParticleUtils::ParticleTileType = typename WarpXParticleContainer::ParticleTileType

◆ ParticleType

Function Documentation

◆ containsInclusive()

AMREX_GPU_HOST_DEVICE AMREX_INLINE bool ParticleUtils::containsInclusive ( amrex::RealBox const & tilebox,
amrex::XDim3 const point )

◆ doLorentzTransform()

AMREX_GPU_HOST_DEVICE AMREX_INLINE void ParticleUtils::doLorentzTransform ( amrex::ParticleReal & ux,
amrex::ParticleReal & uy,
amrex::ParticleReal & uz,
amrex::ParticleReal const Vx,
amrex::ParticleReal const Vy,
amrex::ParticleReal const Vz )

Perform a Lorentz transformation of the given velocity to a frame moving with velocity (Vx, Vy, Vz) relative to the present one.

Parameters
[in,out]ux,uy,uzcomponents of velocity vector in the current frame - importantly these quantities are gamma * velocity
[in]Vx,Vy,Vzvelocity of the new frame relative to the current one, NOT gamma*velocity!

◆ doLorentzTransformWithP()

AMREX_GPU_HOST_DEVICE AMREX_INLINE void ParticleUtils::doLorentzTransformWithP ( amrex::ParticleReal & px,
amrex::ParticleReal & py,
amrex::ParticleReal & pz,
amrex::ParticleReal mass,
amrex::ParticleReal const Ux,
amrex::ParticleReal const Uy,
amrex::ParticleReal const Uz )

Perform a Lorentz transformation of the given momentum to a frame moving with gamma*velocity (Ux, Uy, Uz) relative to the present one.

Parameters
[in,out]px,py,pzcomponents of momentum vector in the current frame
[in]massthe particle mass
[in]Ux,Uy,Uzvelocity of the new frame relative to the current one, importantly these quantities are gamma * velocity

◆ doLorentzTransformWithU()

AMREX_GPU_HOST_DEVICE AMREX_INLINE void ParticleUtils::doLorentzTransformWithU ( amrex::ParticleReal & ux,
amrex::ParticleReal & uy,
amrex::ParticleReal & uz,
amrex::ParticleReal const Ux,
amrex::ParticleReal const Uy,
amrex::ParticleReal const Uz )

Perform a Lorentz transformation of the given velocity to a frame moving with gamma*velocity (Ux, Uy, Uz) relative to the present one.

Parameters
[in,out]ux,uy,uzcomponents of velocity vector in the current frame - importantly these quantities are gamma * velocity
[in]Ux,Uy,Uzvelocity of the new frame relative to the current one, importantly these quantities are gamma * velocity

◆ findParticlesInEachCell()

amrex::DenseBins< ParticleTileDataType > ParticleUtils::findParticlesInEachCell ( amrex::Geometry const & geom_lev,
amrex::MFIter const & mfi,
WarpXParticleContainer::ParticleTileType & ptile )

Find the particles and count the particles that are in each cell. More specifically this function returns an amrex::DenseBins object containing an offset array and a permutation array which can be used to loop over all the cells in a tile and apply an algorithm to particles of a given species present in each cell. Note that this does not rearrange particle arrays.

Parameters
[in]geom_levthe geometry of the current refinement level.
[in]mfithe MultiFAB iterator.
[in]ptilethe particle tile.

◆ getCollisionEnergy()

AMREX_GPU_HOST_DEVICE AMREX_INLINE void ParticleUtils::getCollisionEnergy ( amrex::ParticleReal const u2,
double const m,
double const M,
double & gamma,
double & energy )

Return (relativistic) collision energy assuming the target (with mass M) is stationary and the projectile is approaching with the the given speed and mass m. Note the use of double since this calculation is prone to error with single precision.

Parameters
[in]u2square of particle speed (i.e. u dot u where u = gamma*v)
[in]m,Mmass of projectile and target, respectively
[out]gammarelativistic factor of the particle
[out]energyparticle energy in eV

◆ getRandomVector()

AMREX_GPU_HOST_DEVICE AMREX_INLINE void ParticleUtils::getRandomVector ( amrex::ParticleReal & x,
amrex::ParticleReal & y,
amrex::ParticleReal & z,
amrex::RandomEngine const & engine )

Generate random unit vector in 3 dimensions https://mathworld.wolfram.com/SpherePointPicking.html.

Parameters
[out]xx-component of resulting random vector
[out]yy-component of resulting random vector
[out]zz-component of resulting random vector
[in]enginethe random-engine

◆ ModifyEnergyPairwise()

template<typename T_index>
AMREX_GPU_HOST_DEVICE AMREX_INLINE bool ParticleUtils::ModifyEnergyPairwise ( amrex::ParticleReal *const AMREX_RESTRICT uxp,
amrex::ParticleReal *const AMREX_RESTRICT uyp,
amrex::ParticleReal *const AMREX_RESTRICT uzp,
amrex::ParticleReal *const AMREX_RESTRICT w,
T_index const *const AMREX_RESTRICT indices,
T_index const cell_start,
T_index const cell_stop,
amrex::ParticleReal const mass,
amrex::ParticleReal const energy_fraction,
amrex::ParticleReal const energy_fraction_max,
amrex::ParticleReal & deltaE )

◆ RandomizeVelocity()

AMREX_GPU_HOST_DEVICE AMREX_INLINE void ParticleUtils::RandomizeVelocity ( amrex::ParticleReal & ux,
amrex::ParticleReal & uy,
amrex::ParticleReal & uz,
const amrex::ParticleReal vp,
amrex::RandomEngine const & engine )

Function to perform scattering of a particle that results in a random velocity vector with given magnitude. This is used in isotropic collision events.

Parameters
[in,out]ux,uy,uzcolliding particle's velocity
[in]vpvelocity magnitude of the colliding particle after collision.
[in]enginethe random-engine