Starting with the pusher, consistently use ParticleReal (#3259)

* Starting with the pusher, consistently use ParticleReal

* Update benchmarks for background_mcc_dp_psp

1 files changed, 14 insertions, 14 deletions

diff --git a/Source/Particles/Pusher/UpdateMomentumBoris.H b/Source/Particles/Pusher/UpdateMomentumBoris.H
index ecb379956..c26a97aa4 100644
--- a/Source/Particles/Pusher/UpdateMomentumBoris.H
+++ b/Source/Particles/Pusher/UpdateMomentumBoris.H
@@ -17,31 +17,31 @@ void UpdateMomentumBoris(
     amrex::ParticleReal& ux, amrex::ParticleReal& uy, amrex::ParticleReal& uz,
     const amrex::ParticleReal Ex, const amrex::ParticleReal Ey, const amrex::ParticleReal Ez,
     const amrex::ParticleReal Bx, const amrex::ParticleReal By, const amrex::ParticleReal Bz,
-    const amrex::Real q, const amrex::Real m, const amrex::Real dt )
+    const amrex::ParticleReal q, const amrex::ParticleReal m, const amrex::Real dt )
 {
     using namespace amrex::literals;
 
-    const amrex::Real econst = 0.5_rt*q*dt/m;
+    const amrex::ParticleReal econst = 0.5_prt*q*dt/m;
 
     // First half-push for E
     ux += econst*Ex;
     uy += econst*Ey;
     uz += econst*Ez;
     // Compute temporary gamma factor
-    constexpr amrex::Real inv_c2 = 1._rt/(PhysConst::c*PhysConst::c);
-    const amrex::Real inv_gamma = 1._rt/std::sqrt(1._rt + (ux*ux + uy*uy + uz*uz)*inv_c2);
+    constexpr amrex::ParticleReal inv_c2 = 1._prt/(PhysConst::c*PhysConst::c);
+    const amrex::ParticleReal inv_gamma = 1._prt/std::sqrt(1._prt + (ux*ux + uy*uy + uz*uz)*inv_c2);
     // Magnetic rotation
     // - Compute temporary variables
-    const amrex::Real tx = econst*inv_gamma*Bx;
-    const amrex::Real ty = econst*inv_gamma*By;
-    const amrex::Real tz = econst*inv_gamma*Bz;
-    const amrex::Real tsqi = 2._rt/(1._rt + tx*tx + ty*ty + tz*tz);
-    const amrex::Real sx = tx*tsqi;
-    const amrex::Real sy = ty*tsqi;
-    const amrex::Real sz = tz*tsqi;
-    const amrex::Real ux_p = ux + uy*tz - uz*ty;
-    const amrex::Real uy_p = uy + uz*tx - ux*tz;
-    const amrex::Real uz_p = uz + ux*ty - uy*tx;
+    const amrex::ParticleReal tx = econst*inv_gamma*Bx;
+    const amrex::ParticleReal ty = econst*inv_gamma*By;
+    const amrex::ParticleReal tz = econst*inv_gamma*Bz;
+    const amrex::ParticleReal tsqi = 2._prt/(1._prt + tx*tx + ty*ty + tz*tz);
+    const amrex::ParticleReal sx = tx*tsqi;
+    const amrex::ParticleReal sy = ty*tsqi;
+    const amrex::ParticleReal sz = tz*tsqi;
+    const amrex::ParticleReal ux_p = ux + uy*tz - uz*ty;
+    const amrex::ParticleReal uy_p = uy + uz*tx - ux*tz;
+    const amrex::ParticleReal uz_p = uz + ux*ty - uy*tx;
     // - Update momentum
     ux += uy_p*sz - uz_p*sy;
     uy += uz_p*sx - ux_p*sz;