Compute particle distribution

experimental/examples/compare_files.cpp

@@ -40,7 +40,7 @@
 //!             second file name to compare (with extension .fma (ascii) or bfma (binary) </li>
 //!
 //!     <li> --sort
-//!      Sort particles according to their morton index. (0 yes, 1 no) </li>
+//!      Sort particles according to their morton index. </li>
 //!
 //!     <li> --index, --index1 value,value,value,...  (required)
 //!       column index for data associated to particle (x,y,z,input data, output data)

experimental/examples/test-build-let.cpp

@@ -78,7 +78,7 @@ auto main([[maybe_unused]] int argc, [[maybe_unused]] char* argv[]) -> int
     ///   2) each processor read N/P particles
     scalfmm::tools::DistFmaGenericLoader<value_type, dimension> loader(input_file, para);
     //
-    const int number_of_particles = loader.getNumberOfParticles();
+    const std::size_t number_of_particles = loader.getNumberOfParticles();
     const int local_number_of_particles = loader.getMyNumberOfParticles();
     value_type width = loader.getBoxWidth();
     auto centre = loader.getBoxCenter();
@@ -122,11 +122,12 @@ auto main([[maybe_unused]] int argc, [[maybe_unused]] char* argv[]) -> int
 
     int rank = para.get_process_id();
 
+    int leaf_level = tree_height - 1;
 
-        scalfmm::tree::let::buildLetTree(para, loader, particles_set, box, tree_height, group_size,
+        scalfmm::tree::let::buildLetTree(para, number_of_particles, particles_set, box, leaf_level, group_size,
                                      group_size, letGroupTree,
                                      mortonCellDistribution, nb_block);
-        scalfmm::tree::print("rank(" + std::to_string(rank) + ") data distrib: ", mortonCellDistribution);
+      //  scalfmm::tree::print("rank(" + std::to_string(rank) + ") data distrib: ", mortonCellDistribution);
 
 
     ///

experimental/include/scalfmm/utils/io_helpers.hpp

@@ -4,6 +4,7 @@
 #include "inria/integer_sequence.hpp"
 #include "inria/ostream_joiner.hpp"
 
+#include <array>
 #include <ostream>
 #include <tuple>