tracing

src/core/cell_system/CellStructure.hpp

@@ -544,6 +544,10 @@ struct CellStructure : public System::Leaf<CellStructure> {
     for (auto &p : Cells::particles(decomposition->local_cells())) {
       add_particle(std::move(p));
     }
+    for (auto &p : local_particles()) {
+      std::cout << p.id()<< " ";
+    }
+    std::cout << std::endl;
   }
 
 public:

src/core/cell_system/RegularDecomposition.cpp

@@ -443,7 +443,7 @@ void RegularDecomposition::init_cell_interactions() {
   if (fully_connected_boundary()) {
     if ((*fully_connected_boundary()).first ==
         (*fully_connected_boundary()).second) {
-      throw std::runtime_error("fully_connectd_boudnary normal and connection "
+      throw std::runtime_error("fully_connectd_boundary normal and connection "
                                "coordinates need to differ.");
     }
     if (node_grid[(*fully_connected_boundary()).second] != 1) {
@@ -674,6 +674,11 @@ RegularDecomposition::RegularDecomposition(
     boost::optional<std::pair<int, int>> fully_connected)
     : m_comm(std::move(comm)), m_box(box_geo), m_local_box(local_geo),
       m_fully_connected_boundary(fully_connected) {
+
+  std::cout << "Regular: Range: "<<range<<std::endl;
+  if (m_fully_connected_boundary) {
+     std::cout << "fcb "<<(*m_fully_connected_boundary).first<<" "<<(*m_fully_connected_boundary).second<<std::endl;
+  }
   /* set up new regular decomposition cell structure */
   create_cell_grid(range);
 

src/core/forces.cpp

@@ -221,6 +221,7 @@ void System::System::calculate_forces() {
     lb_couple_particles();
   }
 
+  particles = cell_structure->local_particles();
 #ifdef CUDA
 #ifdef CALIPER
   CALI_MARK_BEGIN("copy_forces_from_GPU");

src/script_interface/cell_system/CellSystem.cpp

@@ -275,7 +275,7 @@ void CellSystem::initialize(CellStructureType const &cs_type,
     m_cell_structure->set_hybrid_decomposition(cutoff_regular, n_square_types);
   } else if (cs_type == CellStructureType::REGULAR) {
     boost::optional<std::pair<int, int>> fcb_pair = {};
-    if (params.count("fully_conected_boundary") && is_type<None>(params.at("fully_connected_boudnary"))) {
+    if (params.count("fully_connected_boundary") && is_type<None>(params.at("fully_connected_boundary"))) {
         // pass
     } else {
       auto const fully_connected_boundary =   get_value_or<std::vector<int>>(params, "fully_connected_boundary", {});