cpp

  // Get cells bounding boxes.
  for (const auto &cell_ite : input_cells) {
    const Tetrahedron tetra(
        cell_ite->vertex(0)->point(), cell_ite->vertex(1)->point(),
        cell_ite->vertex(2)->point(), cell_ite->vertex(3)->point());
    boxes.push_back(Box(std::move(tetra.bbox()), cell_ite));
  }

  // Get the intersections
  CGAL::box_self_intersection_d(boxes.begin(), boxes.end(), called);