diff options
| author |  Edoardo Zoni <59625522+EZoni@users.noreply.github.com>
| 2021-03-04 10:23:06 -0800 |
|---|---|---|
| committer |  GitHub <noreply@github.com>
| 2021-03-04 10:23:06 -0800 |
| commit | f51c8bdf36077b3f53cc60ea9633f42a6346f103 (patch) | |
| tree | 5b6f92bf03d9256346fce5a9403b4400b8a4ab09 /Source/WarpX.cpp | |
| parent | 279e5266ab50977ba8ea7a3ca20e848eb9ba16d2 (diff) | |
| download | WarpX-f51c8bdf36077b3f53cc60ea9633f42a6346f103.tar.gz WarpX-f51c8bdf36077b3f53cc60ea9633f42a6346f103.tar.zst WarpX-f51c8bdf36077b3f53cc60ea9633f42a6346f103.zip | |
Remove extra ghost cell used with momentum-conserving gathering (#1758)
Diffstat (limited to '')
| -rw-r--r-- | Source/WarpX.cpp | 18 |
1 files changed, 8 insertions, 10 deletions
diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp index d24afce67..5d2c1e9e9 100644 --- a/Source/WarpX.cpp +++ b/Source/WarpX.cpp @@ -1046,7 +1046,6 @@ WarpX::AllocLevelData (int lev, const BoxArray& ba, const DistributionMapping& d WarpX::use_fdtd_nci_corr, do_nodal, do_moving_window, - aux_is_nodal, moving_window_dir, WarpX::nox, nox_fft, noy_fft, noz_fft, @@ -1075,14 +1074,13 @@ WarpX::AllocLevelData (int lev, const BoxArray& ba, const DistributionMapping& d } AllocLevelMFs(lev, ba, dm, guard_cells.ng_alloc_EB, guard_cells.ng_alloc_J, - guard_cells.ng_alloc_Rho, guard_cells.ng_alloc_F, - guard_cells.ng_Extra, aux_is_nodal); + guard_cells.ng_alloc_Rho, guard_cells.ng_alloc_F, aux_is_nodal); } void WarpX::AllocLevelMFs (int lev, const BoxArray& ba, const DistributionMapping& dm, const IntVect& ngE, const IntVect& ngJ, const IntVect& ngRho, - const IntVect& ngF, const IntVect& ngextra, const bool aux_is_nodal) + const IntVect& ngF, const bool aux_is_nodal) { // Declare nodal flags IntVect Ex_nodal_flag, Ey_nodal_flag, Ez_nodal_flag; @@ -1164,13 +1162,13 @@ WarpX::AllocLevelMFs (int lev, const BoxArray& ba, const DistributionMapping& dm // std::array<Real,3> dx = CellSize(lev); - Bfield_fp[lev][0] = std::make_unique<MultiFab>(amrex::convert(ba,Bx_nodal_flag),dm,ncomps,ngE+ngextra,tag("Bfield_fp[x]")); - Bfield_fp[lev][1] = std::make_unique<MultiFab>(amrex::convert(ba,By_nodal_flag),dm,ncomps,ngE+ngextra,tag("Bfield_fp[y]")); - Bfield_fp[lev][2] = std::make_unique<MultiFab>(amrex::convert(ba,Bz_nodal_flag),dm,ncomps,ngE+ngextra,tag("Bfield_fp[z]")); + Bfield_fp[lev][0] = std::make_unique<MultiFab>(amrex::convert(ba,Bx_nodal_flag),dm,ncomps,ngE,tag("Bfield_fp[x]")); + Bfield_fp[lev][1] = std::make_unique<MultiFab>(amrex::convert(ba,By_nodal_flag),dm,ncomps,ngE,tag("Bfield_fp[y]")); + Bfield_fp[lev][2] = std::make_unique<MultiFab>(amrex::convert(ba,Bz_nodal_flag),dm,ncomps,ngE,tag("Bfield_fp[z]")); - Efield_fp[lev][0] = std::make_unique<MultiFab>(amrex::convert(ba,Ex_nodal_flag),dm,ncomps,ngE+ngextra,tag("Efield_fp[x]")); - Efield_fp[lev][1] = std::make_unique<MultiFab>(amrex::convert(ba,Ey_nodal_flag),dm,ncomps,ngE+ngextra,tag("Efield_fp[y]")); - Efield_fp[lev][2] = std::make_unique<MultiFab>(amrex::convert(ba,Ez_nodal_flag),dm,ncomps,ngE+ngextra,tag("Efield_fp[z]")); + Efield_fp[lev][0] = std::make_unique<MultiFab>(amrex::convert(ba,Ex_nodal_flag),dm,ncomps,ngE,tag("Efield_fp[x]")); + Efield_fp[lev][1] = std::make_unique<MultiFab>(amrex::convert(ba,Ey_nodal_flag),dm,ncomps,ngE,tag("Efield_fp[y]")); + Efield_fp[lev][2] = std::make_unique<MultiFab>(amrex::convert(ba,Ez_nodal_flag),dm,ncomps,ngE,tag("Efield_fp[z]")); current_fp[lev][0] = std::make_unique<MultiFab>(amrex::convert(ba,jx_nodal_flag),dm,ncomps,ngJ,tag("current_fp[x]")); current_fp[lev][1] = std::make_unique<MultiFab>(amrex::convert(ba,jy_nodal_flag),dm,ncomps,ngJ,tag("current_fp[y]")); |