Handle boundary vertices in angle_sum

2025-01-10 13:10:33 +01:00 · 2025-01-10 13:10:33 +01:00 · db297d7eff
parent 26005e59f3
commit db297d7eff
1 changed files with 3 additions and 0 deletions
--- a/Polygon_mesh_processing/include/CGAL/Polygon_mesh_processing/curvature.h
+++ b/Polygon_mesh_processing/include/CGAL/Polygon_mesh_processing/curvature.h
@ -92,6 +92,9 @@ angle_sum(typename boost::graph_traits<PolygonMesh>::vertex_descriptor v,
  FT angle_sum = 0;
  for(auto h : halfedges_around_source(v, pmesh))
  {
+    if(is_border(h, pmesh))
+      continue;
+
    angle_sum += approx_angle(get(vpm, target(h, pmesh)),
                              get(vpm, source(h, pmesh)),
                              get(vpm, source(prev(h,pmesh), pmesh)));