From fd3f76294b9edfe4a39ae5594b46c1c20308753a Mon Sep 17 00:00:00 2001 From: Luca Fedeli Date: Fri, 12 Mar 2021 19:54:49 +0100 Subject: AddPlasma: Avoid a calculation resulting in NaN when num_ppc is 0 (#1787) * avoid a NaN when num_ppc is 0 * fix style --- Source/Particles/PhysicalParticleContainer.cpp | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) (limited to 'Source/Particles/PhysicalParticleContainer.cpp') diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp index f0320d00f..412ed53d4 100644 --- a/Source/Particles/PhysicalParticleContainer.cpp +++ b/Source/Particles/PhysicalParticleContainer.cpp @@ -541,12 +541,15 @@ PhysicalParticleContainer::AddPlasma (int lev, RealBox part_realbox) const auto dx = geom.CellSizeArray(); const auto problo = geom.ProbLoArray(); - Real scale_fac; + Real scale_fac = 0.0_rt; + + if(num_ppc != 0){ #if AMREX_SPACEDIM==3 - scale_fac = dx[0]*dx[1]*dx[2]/num_ppc; + scale_fac = dx[0]*dx[1]*dx[2]/num_ppc; #elif AMREX_SPACEDIM==2 - scale_fac = dx[0]*dx[1]/num_ppc; + scale_fac = dx[0]*dx[1]/num_ppc; #endif + } defineAllParticleTiles(); -- cgit v1.2.3