7 #ifndef WARPX_PARTICLES_PUSHER_UPDATEMOMENTUM_HIGUERACARY_H_
8 #define WARPX_PARTICLES_PUSHER_UPDATEMOMENTUM_HIGUERACARY_H_
24 const T Ex,
const T Ey,
const T Ez,
25 const T Bx,
const T By,
const T Bz,
26 const T q,
const T m,
const amrex::Real
dt )
28 using namespace amrex::literals;
31 const T qmt = 0.5_prt*q*
dt/m;
35 const T umx = ux + qmt*Ex;
36 const T umy = uy + qmt*Ey;
37 const T umz = uz + qmt*Ez;
39 T
gamma = 1._prt + (umx*umx + umy*umy + umz*umz)*invclightsq;
41 const T betax = qmt*Bx;
42 const T betay = qmt*By;
43 const T betaz = qmt*Bz;
44 const T betam = betax*betax + betay*betay + betaz*betaz;
46 const T sigma =
gamma - betam;
48 const T ust = (umx*betax + umy*betay + umz*betaz)*invclight;
50 gamma = 1._prt/std::sqrt(0.5_prt*(sigma + std::sqrt(sigma*sigma + 4._prt*(betam + ust*ust)) ));
52 const T tx =
gamma*betax;
53 const T ty =
gamma*betay;
54 const T tz =
gamma*betaz;
56 const T s = 1._prt/(1._prt+(tx*tx + ty*ty + tz*tz));
58 const T umt = umx*tx + umy*ty + umz*tz;
60 const T upx = s*( umx + umt*tx + umy*tz - umz*ty );
61 const T upy = s*( umy + umt*ty + umz*tx - umx*tz );
62 const T upz = s*( umz + umt*tz + umx*ty - umy*tx );
64 ux = upx + qmt*Ex + upy*tz - upz*ty;
65 uy = upy + qmt*Ey + upz*tx - upx*tz;
66 uz = upz + qmt*Ez + upx*ty - upy*tx;
#define AMREX_GPU_HOST_DEVICE
AMREX_GPU_HOST_DEVICE AMREX_INLINE void UpdateMomentumHigueraCary(T &ux, T &uy, T &uz, const T Ex, const T Ey, const T Ez, const T Bx, const T By, const T Bz, const T q, const T m, const amrex::Real dt)
Push the particle's positions over one timestep, given the value of its momenta ux,...
Definition: UpdateMomentumHigueraCary.H:22
static constexpr auto c
vacuum speed of light [m/s]
Definition: constant.H:44
float dt
Definition: stencil.py:442
int gamma
boosted frame
Definition: stencil.py:431