Merge pull request #329 from JamesMcClung/pr/fix-sparse-particles

Fix particle weight when 0<n<.5
psc-code · Apr 17, 2024 · 23f175a · 23f175a
2 parents 62a4ece + da98c4f
commit 23f175a
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/include/setup_particles.hxx b/src/include/setup_particles.hxx
@@ -1,6 +1,7 @@
 
 #pragma once
 
+#include <algorithm>
 #include <functional>
 #include <type_traits>
 #include "centering.hxx"
@@ -122,7 +123,7 @@ struct SetupParticles
     if (fractional_n_particles_per_cell) {
       return np.n / norm_.cori + dist.get();
     }
-    return np.n / norm_.cori + .5;
+    return std::max(1, int(np.n / norm_.cori + .5));
   }
 
   // ----------------------------------------------------------------------