From 23793e198743d086145a6ff378f3f625562988ef Mon Sep 17 00:00:00 2001 From: Edoardo Zoni <59625522+EZoni@users.noreply.github.com> Date: Wed, 25 Aug 2021 16:30:19 -0700 Subject: Multi-J: Optimize Allocations of rho MultiFabs (#2224) * OneStep_multiJ: Always Deposit rho in Component 0 * For Multi-J Algo Allocate only 1 Component Co-authored-by: Remi Lehe --- Source/Particles/MultiParticleContainer.cpp | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) (limited to 'Source/Particles/MultiParticleContainer.cpp') diff --git a/Source/Particles/MultiParticleContainer.cpp b/Source/Particles/MultiParticleContainer.cpp index 5cf9bcb34..549358f30 100644 --- a/Source/Particles/MultiParticleContainer.cpp +++ b/Source/Particles/MultiParticleContainer.cpp @@ -521,13 +521,12 @@ MultiParticleContainer::DepositCurrent ( void MultiParticleContainer::DepositCharge ( amrex::Vector >& rho, - const amrex::Real relative_t, const int icomp) + const amrex::Real relative_t) { // Reset the rho array for (int lev = 0; lev < rho.size(); ++lev) { - int const nc = WarpX::ncomps; - rho[lev]->setVal(0.0, icomp*nc, nc, rho[lev]->nGrowVect()); + rho[lev]->setVal(0.0, 0, WarpX::ncomps, rho[lev]->nGrowVect()); } // Push the particles in time, if needed @@ -541,7 +540,7 @@ MultiParticleContainer::DepositCharge ( bool const do_rz_volume_scaling = false; bool const interpolate_across_levels = false; pc->DepositCharge(rho, local, reset, do_rz_volume_scaling, - interpolate_across_levels, icomp); + interpolate_across_levels); } // Push the particles back in time -- cgit v1.2.3