From 0a6df0877a5c5f9069b8985fbd98a003eae974ab Mon Sep 17 00:00:00 2001 From: Jane Tournois Date: Thu, 13 Nov 2025 15:38:20 +0100 Subject: [PATCH] wip using locate_corner() --- Mesh_3/include/CGAL/Sizing_field_with_aabb_tree.h | 13 +++++++++---- 1 file changed, 9 insertions(+), 4 deletions(-) diff --git a/Mesh_3/include/CGAL/Sizing_field_with_aabb_tree.h b/Mesh_3/include/CGAL/Sizing_field_with_aabb_tree.h index 74346f5ce94..2471f6028b1 100644 --- a/Mesh_3/include/CGAL/Sizing_field_with_aabb_tree.h +++ b/Mesh_3/include/CGAL/Sizing_field_with_aabb_tree.h @@ -228,7 +228,8 @@ public: //fill incidences of corners with curves d_ptr->corners_incident_curves.resize(d_ptr->corners.size()); - for(const typename Corners_indices::value_type& pair : d_ptr->corners_indices) { + for(const typename Corners_indices::value_type& pair : d_ptr->corners_indices) + { d_ptr->dt.insert(pair.first); // Fill `corners_incident_curves[corner_id]` @@ -236,15 +237,19 @@ public: d_ptr->domain.get_corner_incident_curves(pair.second, std::inserter(incident_curves, incident_curves.end())); - // For each incident loops, insert a point on the loop, as far as + // For each incident loop, insert a point on the loop, as far as // possible. - for(Curve_index curve_index : incident_curves) { + for(Curve_index curve_index : incident_curves) + { if(domain.is_loop(curve_index)) { FT curve_length = d_ptr->domain.curve_length(curve_index); + auto loc = + d_ptr->domain.locate_corner(curve_index, pair.first); Point_3 other_point = d_ptr->domain.construct_point_on_curve(pair.first, curve_index, - curve_length / 2); + curve_length / 2, + loc); d_ptr->dt.insert(other_point); } }