Add initial side length to map when building octree

CGAL · Aug 6, 2020 · 25c1aaa · 25c1aaa
1 parent 6cca8b2
commit 25c1aaa
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/Octree/include/CGAL/Octree.h b/Octree/include/CGAL/Octree.h
@@ -191,6 +191,7 @@ class Octree {
     m_bbox_min = bbox.min();
     m_bbox_side = bbox.max()[0] - m_bbox_min[0];
     m_root.points() = {point_range.begin(), point_range.end()};
+    m_side_per_depth.push_back(m_bbox_side);
 
   }
 
@@ -217,6 +218,9 @@ class Octree {
     // Make sure the max depth is reset as well
     m_max_depth_reached = 0;
 
+    // Reset the side length map, too
+    m_side_per_depth.resize(0);
+
     // create a side length map
     for (int i = 0; i <= (int) 32; i++)
       m_side_per_depth.push_back(m_bbox_side / (FT) (1 << i));

diff --git a/Octree/include/CGAL/Octree/IO.h b/Octree/include/CGAL/Octree/IO.h
@@ -22,7 +22,6 @@
 
 using std::ostream;
 
-
 template<typename Value>
 ostream &operator<<(ostream &os, const CGAL::Octree::Node<Value> &node) {