From 3c3bd8f98deb864c0b1f1f7cd6693d92e46fa044 Mon Sep 17 00:00:00 2001 From: Menno Fraters Date: Fri, 24 Jan 2020 18:39:20 -0800 Subject: [PATCH] Fix particle probability density mpi problem --- source/particles/generators.cc | 11 +++++++---- 1 file changed, 7 insertions(+), 4 deletions(-) diff --git a/source/particles/generators.cc b/source/particles/generators.cc index 1d3713b255..60d1b36f41 100644 --- a/source/particles/generators.cc +++ b/source/particles/generators.cc @@ -304,13 +304,16 @@ namespace Particles } #endif - // Calculate start id and number of local particles + // Calculate start id start_particle_id = std::llround(static_cast(n_particles_to_create) * local_start_weight / global_weight_integral); - n_local_particles = - std::llround(static_cast(n_particles_to_create) * - local_weight_integral / global_weight_integral); + + // Calcualate number of local particles + const types::particle_index end_particle_id = + llround(static_cast(n_particles_to_create) * + local_end_weight / global_weight_integral); + n_local_particles = end_particle_id - start_particle_id; if (random_cell_selection) { -- 2.39.5