diff options
Diffstat (limited to 'Source/Particles/PhysicalParticleContainer.cpp')
-rw-r--r-- | Source/Particles/PhysicalParticleContainer.cpp | 20 |
1 files changed, 10 insertions, 10 deletions
diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp index f4cdbcb92..0daa9fb5f 100644 --- a/Source/Particles/PhysicalParticleContainer.cpp +++ b/Source/Particles/PhysicalParticleContainer.cpp @@ -1000,7 +1000,7 @@ PhysicalParticleContainer::AssignExternalFieldOnParticles(WarpXParIter& pti, Real time = warpx.gett_new(lev); amrex::ParallelFor(pti.numParticles(), [=] AMREX_GPU_DEVICE (long i) { - Real x, y, z; + ParticleReal x, y, z; get_position(i, x, y, z); Exp_data[i] = xfield_partparser->getField(x, y, z, time); Eyp_data[i] = yfield_partparser->getField(x, y, z, time); @@ -1022,7 +1022,7 @@ PhysicalParticleContainer::AssignExternalFieldOnParticles(WarpXParIter& pti, Real time = warpx.gett_new(lev); amrex::ParallelFor(pti.numParticles(), [=] AMREX_GPU_DEVICE (long i) { - Real x, y, z; + ParticleReal x, y, z; get_position(i, x, y, z); Bxp_data[i] = xfield_partparser->getField(x, y, z, time); Byp_data[i] = yfield_partparser->getField(x, y, z, time); @@ -1489,7 +1489,7 @@ PhysicalParticleContainer::SplitParticles(int lev) auto& uzp = attribs[PIdx::uz]; const long np = pti.numParticles(); for(int i=0; i<np; i++){ - Real xp, yp, zp; + ParticleReal xp, yp, zp; get_position(i, xp, yp, zp); auto& p = particles[i]; if (p.id() == DoSplitParticleID){ @@ -1667,7 +1667,7 @@ PhysicalParticleContainer::PushPX (WarpXParIter& pti, Real dt, DtType a_dt_type) Ex[i], Ey[i], Ez[i], Bx[i], By[i], Bz[i], q, m, dt); } - Real x, y, z; + ParticleReal x, y, z; get_position(i, x, y, z); UpdatePosition(x, y, z, ux[i], uy[i], uz[i], dt ); set_position(i, x, y, z); @@ -1680,7 +1680,7 @@ PhysicalParticleContainer::PushPX (WarpXParIter& pti, Real dt, DtType a_dt_type) UpdateMomentumBorisWithRadiationReaction( ux[i], uy[i], uz[i], Ex[i], Ey[i], Ez[i], Bx[i], By[i], Bz[i], q, m, dt); - Real x, y, z; + ParticleReal x, y, z; get_position(i, x, y, z); UpdatePosition(x, y, z, ux[i], uy[i], uz[i], dt ); set_position(i, x, y, z); @@ -1697,7 +1697,7 @@ PhysicalParticleContainer::PushPX (WarpXParIter& pti, Real dt, DtType a_dt_type) UpdateMomentumBorisWithRadiationReaction( ux[i], uy[i], uz[i], Ex[i], Ey[i], Ez[i], Bx[i], By[i], Bz[i], qp, m, dt); - Real x, y, z; + ParticleReal x, y, z; get_position(i, x, y, z); UpdatePosition(x, y, z, ux[i], uy[i], uz[i], dt ); set_position(i, x, y, z); @@ -1713,7 +1713,7 @@ PhysicalParticleContainer::PushPX (WarpXParIter& pti, Real dt, DtType a_dt_type) UpdateMomentumBoris( ux[i], uy[i], uz[i], Ex[i], Ey[i], Ez[i], Bx[i], By[i], Bz[i], qp, m, dt); - Real x, y, z; + ParticleReal x, y, z; get_position(i, x, y, z); UpdatePosition(x, y, z, ux[i], uy[i], uz[i], dt ); set_position(i, x, y, z); @@ -1728,7 +1728,7 @@ PhysicalParticleContainer::PushPX (WarpXParIter& pti, Real dt, DtType a_dt_type) UpdateMomentumVay( ux[i], uy[i], uz[i], Ex[i], Ey[i], Ez[i], Bx[i], By[i], Bz[i], qp, m, dt); - Real x, y, z; + ParticleReal x, y, z; get_position(i, x, y, z); UpdatePosition(x, y, z, ux[i], uy[i], uz[i], dt ); set_position(i, x, y, z); @@ -1743,7 +1743,7 @@ PhysicalParticleContainer::PushPX (WarpXParIter& pti, Real dt, DtType a_dt_type) UpdateMomentumHigueraCary( ux[i], uy[i], uz[i], Ex[i], Ey[i], Ez[i], Bx[i], By[i], Bz[i], qp, m, dt); - Real x, y, z; + ParticleReal x, y, z; get_position(i, x, y, z); UpdatePosition(x, y, z, ux[i], uy[i], uz[i], dt ); set_position(i, x, y, z); @@ -1937,7 +1937,7 @@ void PhysicalParticleContainer::copy_attribs (WarpXParIter& pti) ParallelFor( np, [=] AMREX_GPU_DEVICE (long i) { - Real x, y, z; + ParticleReal x, y, z; get_position(i, x, y, z); xpold[i]=x; ypold[i]=y; |