如何正确修改CGAL Surface Mesh顶点坐标？修改后计算出现未定义结果

首先，你写的修改顶点坐标的代码本身语法上是正确的：

for (auto vertex_iterator : mesh.vertices()) { 
  Point_3 p(1,0,1); 
  mesh.point(vertex_iterator) = p; 
}

但出现未定义行为，几乎都是因为修改坐标后没有同步更新依赖这些位置信息的辅助结构或属性，以下是几个最常见的场景和解决方法：

1. 同步依赖顶点位置的空间查询结构 #

如果你在修改坐标前构建了诸如AABB_tree（用于射线查询、最近点查找等）这类依赖顶点位置的结构，修改坐标后这些结构的缓存数据会失效，必须手动更新：

// 假设之前构建了AABB树
CGAL::AABB_tree<CGAL::AABB_traits<K, CGAL::Surface_mesh_AABB_primitive<K, Surface_mesh>>> tree(faces(mesh).first, faces(mesh).second, mesh);

// 修改顶点坐标后
tree.refit(faces(mesh).first, faces(mesh).second, mesh);
// 或者如果结构变化很大，直接重新构建更稳妥
tree = CGAL::AABB_tree<CGAL::AABB_traits<K, CGAL::Surface_mesh_AABB_primitive<K, Surface_mesh>>>(faces(mesh).first, faces(mesh).second, mesh);

2. 更新自定义属性映射 #

如果你给mesh附加了依赖顶点坐标的自定义属性（比如法线、曲率、纹理坐标关联的位置信息），CGAL不会自动同步这些属性，修改坐标后需要手动重新计算：

// 示例：更新顶点法线属性
auto normals = mesh.add_property_map<Vertex_index, Vector_3>("v:normals").first;
CGAL::Polygon_mesh_processing::compute_vertex_normals(mesh, normals);

3. 确保点类型匹配 #

检查你使用的Point_3类型是否和Surface Mesh的点类型完全一致。比如如果你的mesh定义是：

using K = CGAL::Simple_cartesian<double>;
using Surface_mesh = CGAL::Surface_mesh<K::Point_3>;

那你修改时用的Point_3必须是K::Point_3，避免隐式转换带来的精度或内存布局问题。

4. 避免退化的拓扑结构 #

修改坐标后如果出现了退化的三角形（三点共线、重合），后续的CGAL算法（比如网格简化、布尔运算）可能会触发未定义行为。可以在修改后添加检查：

bool has_degenerate = false;
for (auto face : mesh.faces()) {
  auto vertices = CGAL::vertices_around_face(mesh.halfedge(face), mesh);
  std::vector<K::Point_3> pts;
  for (auto v : vertices) pts.push_back(mesh.point(v));
  if (CGAL::collinear(pts[0], pts[1], pts[2])) {
    has_degenerate = true;
    break;
  }
}
if (has_degenerate) {
  // 处理退化三角形，比如调整坐标或标记删除
}

5. 检查迭代器有效性（遍历修改时的边界情况） #

如果你的遍历逻辑不是单纯修改坐标，而是结合了顶点的添加/删除，要确保迭代器在修改过程中保持有效。不过你的示例代码只是遍历修改坐标，这一点通常没问题，但如果是复杂遍历，建议使用mesh.vertices()的范围遍历（就像你现在写的），而不是手动管理迭代器。

按照以上几点排查后，应该就能解决修改坐标后出现的未定义结果问题了。

内容的提问来源于stack exchange，提问作者carak