diff options
-rw-r--r-- | Docs/source/running_cpp/parameters.rst | 2 | ||||
-rw-r--r-- | Source/Parallelization/WarpXRegrid.cpp | 1 | ||||
-rw-r--r-- | Source/WarpX.cpp | 5 |
3 files changed, 3 insertions, 5 deletions
diff --git a/Docs/source/running_cpp/parameters.rst b/Docs/source/running_cpp/parameters.rst index b0791cc75..59bfbb14e 100644 --- a/Docs/source/running_cpp/parameters.rst +++ b/Docs/source/running_cpp/parameters.rst @@ -608,8 +608,6 @@ Numerics and algorithms - ``direct``: simpler current deposition algorithm, described in the section :doc:`../theory/picsar_theory`. Note that this algorithm is not strictly charge-conserving. - v - * ``algo.charge_deposition`` (`string`, optional) The algorithm for the charge density deposition. Available options are: diff --git a/Source/Parallelization/WarpXRegrid.cpp b/Source/Parallelization/WarpXRegrid.cpp index 649ea80de..2ae167283 100644 --- a/Source/Parallelization/WarpXRegrid.cpp +++ b/Source/Parallelization/WarpXRegrid.cpp @@ -1,6 +1,5 @@ #include <WarpX.H> -#include <WarpXAlgorithmSelection.H> #include <AMReX_BLProfiler.H> using namespace amrex; diff --git a/Source/WarpX.cpp b/Source/WarpX.cpp index a0362dd57..9ba0741bb 100644 --- a/Source/WarpX.cpp +++ b/Source/WarpX.cpp @@ -876,7 +876,8 @@ WarpX::AllocLevelMFs (int lev, const BoxArray& ba, const DistributionMapping& dm // if (aux_is_nodal and !do_nodal) { - BoxArray const nba = amrex::convert(ba,IntVect::TheUnitVector()); + // Create aux multifabs on Nodal Box Array + BoxArray const nba = amrex::convert(ba,IntVect::TheNodeVector()); Bfield_aux[lev][0].reset( new MultiFab(nba,dm,ncomps,ngE)); Bfield_aux[lev][1].reset( new MultiFab(nba,dm,ncomps,ngE)); Bfield_aux[lev][2].reset( new MultiFab(nba,dm,ncomps,ngE)); @@ -967,7 +968,7 @@ WarpX::AllocLevelMFs (int lev, const BoxArray& ba, const DistributionMapping& dm if (n_field_gather_buffer > 0 || mypc->nSpeciesGatherFromMainGrid() > 0) { if (aux_is_nodal) { - BoxArray const& cnba = amrex::convert(cba,IntVect::TheUnitVector()); + BoxArray const& cnba = amrex::convert(cba,IntVect::TheNodeVector()); Bfield_cax[lev][0].reset( new MultiFab(cnba,dm,ncomps,ngE)); Bfield_cax[lev][1].reset( new MultiFab(cnba,dm,ncomps,ngE)); Bfield_cax[lev][2].reset( new MultiFab(cnba,dm,ncomps,ngE)); |