Forgot to let the original test :)

cpp/particles/p2p/p2p_distr_mpi.hpp

@@ -231,7 +231,9 @@ public:
         const long double field_section_width_z = spatial_box_width[IDXC_Z]/(long double)(field_grid_dim[IDXC_Z]);
         const long int limite = static_cast<long int>((field_section_width_z*(long double)(partition_interval_offset_per_proc[dest_proc+1])
                                      - std::numeric_limits<long double>::epsilon())/cutoff_radius);
-        if(limite == nb_cell_levels[IDXC_Z]){
+        if(static_cast<long double>(limite)*cutoff_radius
+                == field_section_width_z*(long double)(partition_interval_offset_per_proc[dest_proc+1])
+                || limite == nb_cell_levels[IDXC_Z]){
             return limite-1;
         }
         return limite;