aboutsummaryrefslogtreecommitdiff
path: root/Source/Particles/PhysicalParticleContainer.cpp
diff options
context:
space:
mode:
authorGravatar MaxThevenet <mthevenet@lbl.gov> 2019-09-19 13:22:24 -0700
committerGravatar MaxThevenet <mthevenet@lbl.gov> 2019-09-19 13:22:24 -0700
commit60917ab3d92cf5f325468ffd15a9c35b101699a6 (patch)
tree94356d000ffa7e7f8a82e608b91ed8317d47779a /Source/Particles/PhysicalParticleContainer.cpp
parentf43f0264e05d022dbd4d0d5cc7d9d511836f3e33 (diff)
downloadWarpX-60917ab3d92cf5f325468ffd15a9c35b101699a6.tar.gz
WarpX-60917ab3d92cf5f325468ffd15a9c35b101699a6.tar.zst
WarpX-60917ab3d92cf5f325468ffd15a9c35b101699a6.zip
do not store old attribs in second particle push when subcycling
Diffstat (limited to 'Source/Particles/PhysicalParticleContainer.cpp')
-rw-r--r--Source/Particles/PhysicalParticleContainer.cpp10
1 files changed, 6 insertions, 4 deletions
diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp
index 4bc0ee16e..2d0b90902 100644
--- a/Source/Particles/PhysicalParticleContainer.cpp
+++ b/Source/Particles/PhysicalParticleContainer.cpp
@@ -942,7 +942,7 @@ PhysicalParticleContainer::Evolve (int lev,
MultiFab* rho, MultiFab* crho,
const MultiFab* cEx, const MultiFab* cEy, const MultiFab* cEz,
const MultiFab* cBx, const MultiFab* cBy, const MultiFab* cBz,
- Real t, Real dt)
+ Real t, Real dt, DtType a_dt_type)
{
BL_PROFILE("PPC::Evolve()");
BL_PROFILE_VAR_NS("PPC::Evolve::Copy", blp_copy);
@@ -1301,7 +1301,7 @@ PhysicalParticleContainer::Evolve (int lev,
// Particle Push
//
BL_PROFILE_VAR_START(blp_ppc_pp);
- PushPX(pti, m_xp[thread_num], m_yp[thread_num], m_zp[thread_num], dt);
+ PushPX(pti, m_xp[thread_num], m_yp[thread_num], m_zp[thread_num], dt, a_dt_type);
BL_PROFILE_VAR_STOP(blp_ppc_pp);
//
@@ -1526,7 +1526,7 @@ PhysicalParticleContainer::PushPX(WarpXParIter& pti,
Cuda::ManagedDeviceVector<Real>& xp,
Cuda::ManagedDeviceVector<Real>& yp,
Cuda::ManagedDeviceVector<Real>& zp,
- Real dt)
+ Real dt, DtType a_dt_type)
{
// This wraps the momentum and position advance so that inheritors can modify the call.
@@ -1545,8 +1545,10 @@ PhysicalParticleContainer::PushPX(WarpXParIter& pti,
const Real* const AMREX_RESTRICT By = attribs[PIdx::By].dataPtr();
const Real* const AMREX_RESTRICT Bz = attribs[PIdx::Bz].dataPtr();
- if (WarpX::do_boosted_frame_diagnostic && do_boosted_frame_diags)
+ Print()<<"a_dt_type "<<int( a_dt_type )<<'\n';
+ if (WarpX::do_boosted_frame_diagnostic && do_boosted_frame_diags && (a_dt_type!=DtType::SecondHalf))
{
+ Print()<<"copy attribs\n";
copy_attribs(pti, x, y, z);
}