setup_particles: fix when 0<n<.5

psc-code · Apr 3, 2024 · da98c4f · da98c4f
1 parent 62a4ece
commit da98c4f
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));
   }
 
   // ----------------------------------------------------------------------