为何转换一致Delaunay/Gabriel三角剖分时触发CGAL断言异常？

这个断言触发于CGAL的一致Delaunay/Gabriel三角剖分细化逻辑中，核心含义是：当当前三角面f与待插入点p无冲突时，f的邻面n必须是定位p后得到的面zone.fh。不满足这个条件，说明三角剖分的内部状态或输入数据存在问题，常见原因如下：

1. 初始约束三角剖分存在非法结构 #

CGAL的约束Delaunay三角剖分（CDT）要求：

• 约束边不能交叉（除非在端点处相交）
• 不存在重复点或长度为0的退化约束边
• 所有约束边都被正确嵌入到三角剖分中

如果初始CDT存在上述问题，在执行一致化操作时，内部的点定位（locate）和冲突检测（test_conflict）会出现逻辑不一致，导致断言失败。

2. 浮点精度误差影响 #

如果使用的是浮点计算核（如Simple_cartesian<double>），浮点舍入误差可能导致：

• 点p的定位结果（zone.fh）与理论上的邻面n不匹配
• 冲突检测的结果出现误判（明明f与p无冲突，但定位结果却不在邻面n内）

这种精度问题会直接触发断言中的等式检查失败。

3. 三角剖分状态被非法修改 #

在调用make_conforming_Delaunay_2或make_conforming_Gabriel_2之前，如果手动修改了CDT的内部结构（比如直接删除面/边、插入点但未维护三角剖分完整性），会导致三角剖分的邻接关系混乱，进而让定位和邻面查询的结果不一致。

4. CGAL版本的已知缺陷 #

CGAL 5.3.1可能存在特定场景下的逻辑漏洞，比如处理某些特殊约束配置时，细化算法的状态判断出现错误。升级到更高版本的CGAL（如5.4及以上）可能修复该问题。

排查与解决步骤 #

• 验证初始CDT合法性：调用cdt.is_valid(true)检查三角剖分状态，同时确保所有约束边无交叉、无重复点/退化边。
• 切换精确计算核：替换浮点核为精确核，比如Exact_predicates_inexact_constructions_kernel，排除精度误差影响。
• 检查前置操作：确保在调用一致化函数前，CDT未被非法修改，所有操作都通过CGAL提供的API完成。
• 升级CGAL版本：尝试更新到最新稳定版CGAL，验证是否修复该断言问题。

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