如何设置2D CDT允许约束相交并解决坐标转double编译错误

替换内核后触发该编译错误的核心原因是：

• 原内核CGAL::Exact_predicates_inexact_constructions_kernel（EPICK）仅保证几何谓词计算精确，构造的坐标值底层存储为double类型，支持直接隐式转换为double
• 新使用的CGAL::Exact_predicates_exact_constructions_kernel（EPECK）是全流程精确计算内核，点坐标的底层类型为CGAL封装的惰性精确数CGAL::Lazy_exact_nt<boost::multiprecision::mpq_rational>，该类型没有定义隐式转double的逻辑，直接赋值给double变量就会抛出你看到的类型转换错误。

方案1：显式调用CGAL官方转换函数（推荐） #

引入CGAL数值工具头文件，使用CGAL::to_double()方法做显式类型转换，这是CGAL官方推荐的精确数值转双精度浮点数的标准写法，不会触发未定义行为：

// 先引入必要头文件
#include <CGAL/number_utils.h>

bool operator()(const Point & p1, const Point & p2) const
{
  double x = CGAL::to_double(p1.x());
  // 其余y坐标、其他点的坐标读取都按相同方式处理
  // 原有业务逻辑保持不变
}

方案2：排查类型定义不一致问题 #

批量替换类型时需要确认Point的类型定义没有残留旧内核的声明，必须和当前使用的三角网点类型保持一致，避免类型混用：

// 错误写法：残留旧内核的Point定义，会和CDTP返回的点类型冲突
// typedef K::Point_2 Point;

// 正确写法：统一使用三角网类对应的Point类型
typedef CDTP::Point Point;

注意事项 #

1. 转成double的数值仅可用于业务侧的非几何判断逻辑（比如输出、自定义非几何计算），不要将转换后的double值传回CGAL做几何谓词判断（比如线段相交、点位置判断、三角剖分相关计算），否则会丢失精确计算的特性，导致约束相交场景下出现计算错误。
2. 如果业务侧需要高频读取点的double坐标，可以在点插入三角网时额外存储一份double格式的坐标备份，减少重复转换的性能开销。

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