diff options
| author |  Remi Lehe <remi.lehe@normalesup.org>
| 2023-08-11 06:20:14 -0700 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2023-08-11 15:20:14 +0200 |
| commit | b10502843780f16b02f247c05af8d52b153c2475 (patch) | |
| tree | 7647b0057d8d8e9a0be785fd0de3cec04d69bd7f /Source/Particles/PhysicalParticleContainer.cpp | |
| parent | 7293f20939ca6e45e29d995408825ecd86dd6421 (diff) | |
| download | WarpX-b10502843780f16b02f247c05af8d52b153c2475.tar.gz WarpX-b10502843780f16b02f247c05af8d52b153c2475.tar.zst WarpX-b10502843780f16b02f247c05af8d52b153c2475.zip | |
Use more consistent particle type in pusher (#4175)
Diffstat (limited to 'Source/Particles/PhysicalParticleContainer.cpp')
| -rw-r--r-- | Source/Particles/PhysicalParticleContainer.cpp | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp index ded222ce8..c5c44f9b6 100644 --- a/Source/Particles/PhysicalParticleContainer.cpp +++ b/Source/Particles/PhysicalParticleContainer.cpp @@ -2497,8 +2497,8 @@ PhysicalParticleContainer::PushP (int lev, Real dt, } // Loop over the particles and update their momentum - const amrex::Real q = this->charge; - const amrex::Real m = this-> mass; + const amrex::ParticleReal q = this->charge; + const amrex::ParticleReal m = this-> mass; const auto pusher_algo = WarpX::particle_pusher_algo; const auto do_crr = do_classical_radiation_reaction; @@ -2535,25 +2535,25 @@ PhysicalParticleContainer::PushP (int lev, Real dt, } if (do_crr) { - amrex::Real qp = q; + amrex::ParticleReal qp = q; if (ion_lev) { qp *= ion_lev[ip]; } UpdateMomentumBorisWithRadiationReaction(ux[ip], uy[ip], uz[ip], Exp, Eyp, Ezp, Bxp, Byp, Bzp, qp, m, dt); } else if (pusher_algo == ParticlePusherAlgo::Boris) { - amrex::Real qp = q; + amrex::ParticleReal qp = q; if (ion_lev) { qp *= ion_lev[ip]; } UpdateMomentumBoris( ux[ip], uy[ip], uz[ip], Exp, Eyp, Ezp, Bxp, Byp, Bzp, qp, m, dt); } else if (pusher_algo == ParticlePusherAlgo::Vay) { - amrex::Real qp = q; + amrex::ParticleReal qp = q; if (ion_lev){ qp *= ion_lev[ip]; } UpdateMomentumVay( ux[ip], uy[ip], uz[ip], Exp, Eyp, Ezp, Bxp, Byp, Bzp, qp, m, dt); } else if (pusher_algo == ParticlePusherAlgo::HigueraCary) { - amrex::Real qp = q; + amrex::ParticleReal qp = q; if (ion_lev){ qp *= ion_lev[ip]; } UpdateMomentumHigueraCary( ux[ip], uy[ip], uz[ip], Exp, Eyp, Ezp, Bxp, |