Function split_by_binarytree

src/structure_dump.cpp:467–488

}

void split_by_binarytree(grid_volume gvol, std::vector<grid_volume> &result_gvs,
                         std::vector<int> &result_ids, const binary_partition *bp) {
  // reached a leaf
  if (bp->is_leaf()) {
    result_gvs.push_back(gvol);
    result_ids.push_back(bp->get_proc_id());
    return;
  }

  const auto &plane = bp->get_plane();
  int split_point = static_cast<int>((plane.pos - gvol.surroundings().in_direction_min(plane.dir)) /
                                         gvol.surroundings().in_direction(plane.dir) *
                                         gvol.num_direction(plane.dir) +
                                     0.5);
  // traverse left branch
  grid_volume left_gvol = gvol.split_at_fraction(false, split_point, plane.dir);
  split_by_binarytree(left_gvol, result_gvs, result_ids, bp->left_tree());

  // traverse right branch
  grid_volume right_gvol = gvol.split_at_fraction(true, split_point, plane.dir);
  split_by_binarytree(right_gvol, result_gvs, result_ids, bp->right_tree());
}

void structure::load_chunk_layout(const char *filename, boundary_region &br) {
  // Load chunk grid_volumes from a file