7 #ifndef WARPX_PARTICLE_UTILS_H_
8 #define WARPX_PARTICLE_UTILS_H_
46 void getEnergy ( amrex::ParticleReal
const u2,
double const mass,
50 using namespace amrex::literals;
69 double const M,
double&
gamma,
double& energy )
72 using namespace amrex::literals;
78 2.0_rt * m * M * u2 / (
gamma + 1.0_rt)
79 / (M + m +
sqrt(m*m + M*M + 2.0_rt * m * M *
gamma))
94 amrex::ParticleReal& uz,
95 amrex::ParticleReal
const Vx, amrex::ParticleReal
const Vy,
96 amrex::ParticleReal
const Vz )
98 using namespace amrex::literals;
102 const auto V2 = (Vx*Vx + Vy*Vy + Vz*Vz);
103 const auto gamma_V = 1.0_prt / std::sqrt(1.0_prt - V2 / c2);
104 const auto gamma_u = std::sqrt(1.0_prt + (ux*ux + uy*uy + uz*uz) / c2);
111 ux = vx * (1.0_prt + (gamma_V - 1.0_prt) * Vx*Vx/V2)
112 + vy * (gamma_V - 1.0_prt) * Vx*Vy/V2
113 + vz * (gamma_V - 1.0_prt) * Vx*Vz/V2
114 - gamma_V * Vx * gamma_u;
116 uy = vy * (1.0_prt + (gamma_V - 1.0_prt) * Vy*Vy/V2)
117 + vx * (gamma_V - 1.0_prt) * Vx*Vy/V2
118 + vz * (gamma_V - 1.0_prt) * Vy*Vz/V2
119 - gamma_V * Vy * gamma_u;
121 uz = vz * (1.0_prt + (gamma_V - 1.0_prt) * Vz*Vz/V2)
122 + vx * (gamma_V - 1.0_prt) * Vx*Vz/V2
123 + vy * (gamma_V - 1.0_prt) * Vy*Vz/V2
124 - gamma_V * Vz * gamma_u;
143 using namespace amrex::literals;
145 auto const theta =
amrex::Random(engine) * 2.0_prt * MathConst::pi;
147 auto const xy =
sqrt(1_prt - z*z);
163 amrex::ParticleReal& uz,
164 const amrex::ParticleReal vp,
167 amrex::ParticleReal
x, y, z;
185 const auto *
const xlo = tilebox.
lo();
186 const auto *
const xhi = tilebox.
hi();
188 && (xlo[1] <= point.
y) && (point.
y <= xhi[1]),
189 && (xlo[2] <= point.
z) && (point.
z <= xhi[2]));
#define AMREX_GPU_HOST_DEVICE
#define AMREX_D_TERM(a, b, c)
AMREX_GPU_HOST_DEVICE const Real * hi() const &noexcept
AMREX_GPU_HOST_DEVICE const Real * lo() const &noexcept
Definition: ParticleUtils.cpp:26
ParticleBins findParticlesInEachCell(int lev, MFIter const &mfi, ParticleTileType &ptile)
Find the particles and count the particles that are in each cell. More specifically this function ret...
Definition: ParticleUtils.cpp:40
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 projec...
Definition: ParticleUtils.H:68
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,...
Definition: ParticleUtils.H:93
AMREX_GPU_HOST_DEVICE AMREX_INLINE void getEnergy(amrex::ParticleReal const u2, double const mass, double &energy)
Return (relativistic) particle energy given velocity and mass. Note the use of double since this calc...
Definition: ParticleUtils.H:46
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 magn...
Definition: ParticleUtils.H:162
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....
Definition: ParticleUtils.H:137
AMREX_GPU_HOST_DEVICE AMREX_INLINE bool containsInclusive(amrex::RealBox const &tilebox, amrex::XDim3 const point)
Definition: ParticleUtils.H:184
static constexpr auto c
vacuum speed of light [m/s]
Definition: constant.H:44
static constexpr auto q_e
elementary charge [C]
Definition: constant.H:50
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE GpuComplex< T > sqrt(const GpuComplex< T > &a_z) noexcept
int gamma
boosted frame
Definition: stencil.py:431