CGAL中CDT插入约束触发is_finite(d)断言异常的规避方法

问题根源 #

触发is_finite(d)断言，本质是构造MP_Float时传入了非有限浮点数（无穷大Inf或非数值NaN），常见诱因：

• 输入点/约束线段的坐标数值过大，计算时溢出
• 存在近乎重合的顶点，导致计算中出现除以极小值的情况
• 大量约束插入后，拓扑计算出现数值不稳定

具体解决办法 #

1. 预处理输入数据 #

• 给所有点坐标做归一化：把点集缩放到[-1e6, 1e6]这类合理区间，避免数值过大溢出。
• 清理重复点：用CGAL的compare函数或1e-8左右的阈值，合并距离极近的顶点，消除病态计算的源头。

2. 调整CGAL内核配置 #

• 替换默认内核为精确构造内核，比如Exact_predicates_exact_constructions_kernel_with_sqrt，用精确有理数运算替代浮点数，从根本上避免精度问题。
• 若必须用近似内核，可把MP_Float换成Gmpzf（需链接GMP库），提升数值处理范围。

3. 优化约束插入逻辑 #

• 不要一次性插完所有约束，分批插入，每插一批就调用triangulation->normalize()或triangulation->optimize()稳定拓扑结构。
• 插约束前先校验：确保线段两端点不重合，且不与已有约束线段出现近乎共线的病态相交。

4. 主动拦截异常数值 #

• 在构造MP_Float前手动检查数值合法性：

  if (!std::isfinite(d)) {
      // 跳过该约束或修正坐标，根据业务逻辑处理
      continue;
  }
  

• 临时调试可禁用CGAL断言（生产环境不推荐）：编译时加宏定义CGAL_NO_ASSERTIONS，但这会掩盖潜在问题，需谨慎。

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