diff options
author | 2019-09-30 11:41:23 -0700 | |
---|---|---|
committer | 2019-09-30 11:41:23 -0700 | |
commit | 9d567430aeb46181f81e3e1f1dbc48b6291b1e95 (patch) | |
tree | aa06411d8b165294fce671ae804899a2578ef44f /Source/Particles/WarpXParticleContainer.cpp | |
parent | 66885c39fd0f58d7e4ad7398ea6a595391ca60c5 (diff) | |
download | WarpX-9d567430aeb46181f81e3e1f1dbc48b6291b1e95.tar.gz WarpX-9d567430aeb46181f81e3e1f1dbc48b6291b1e95.tar.zst WarpX-9d567430aeb46181f81e3e1f1dbc48b6291b1e95.zip |
Generalized direct current deposition allowing different centerings for Jx, y, and z
Diffstat (limited to 'Source/Particles/WarpXParticleContainer.cpp')
-rw-r--r-- | Source/Particles/WarpXParticleContainer.cpp | 11 |
1 files changed, 3 insertions, 8 deletions
diff --git a/Source/Particles/WarpXParticleContainer.cpp b/Source/Particles/WarpXParticleContainer.cpp index 9fd92ceb4..c7dbd09ba 100644 --- a/Source/Particles/WarpXParticleContainer.cpp +++ b/Source/Particles/WarpXParticleContainer.cpp @@ -395,29 +395,24 @@ WarpXParticleContainer::DepositCurrent(WarpXParIter& pti, WarpX::n_rz_azimuthal_modes); } } else { - // Lower corner of tile box physical domain - // Note that this includes guard cells since it is after tilebox.ngrow - std::array<Real, 3> xyzminx = WarpX::LowerCornerWithCentering(tbx, depos_lev); - std::array<Real, 3> xyzminy = WarpX::LowerCornerWithCentering(tby, depos_lev); - std::array<Real, 3> xyzminz = WarpX::LowerCornerWithCentering(tbz, depos_lev); if (WarpX::nox == 1){ doDepositionShapeN<1>( xp, yp, zp, wp.dataPtr() + offset, uxp.dataPtr() + offset, uyp.dataPtr() + offset, uzp.dataPtr() + offset, ion_lev, jx_arr, jy_arr, jz_arr, np_to_depose, dt, dx, - xyzminx, xyzminy, xyzminz, lo, q); + tbx, tby, tbz, depos_lev, lo, q); } else if (WarpX::nox == 2){ doDepositionShapeN<2>( xp, yp, zp, wp.dataPtr() + offset, uxp.dataPtr() + offset, uyp.dataPtr() + offset, uzp.dataPtr() + offset, ion_lev, jx_arr, jy_arr, jz_arr, np_to_depose, dt, dx, - xyzminx, xyzminy, xyzminz, lo, q); + tbx, tby, tbz, depos_lev, lo, q); } else if (WarpX::nox == 3){ doDepositionShapeN<3>( xp, yp, zp, wp.dataPtr() + offset, uxp.dataPtr() + offset, uyp.dataPtr() + offset, uzp.dataPtr() + offset, ion_lev, jx_arr, jy_arr, jz_arr, np_to_depose, dt, dx, - xyzminx, xyzminy, xyzminz, lo, q); + tbx, tby, tbz, depos_lev, lo, q); } } BL_PROFILE_VAR_STOP(blp_deposit); |