可平方的最大浮点值：函数默认参数的安全性及替代方案问询

好问题！咱们从你的核心需求——要让默认的平方距离尽可能大，但又不能出问题——一步步分析：

1. 原写法std::sqrt(std::numeric_limits<double>::max())的安全隐患 #

先理清楚几个关键数值：

• std::numeric_limits<double>::max()（简称DBL_MAX）是双精度浮点数能存储的最大有限值，大概是1.79769e+308。
• 它的平方根sqrt(DBL_MAX)大概是1.34078e+154，这个值本身能被double正常存储，不会直接崩溃。

但真正的坑在你算法里的平方运算：当你用这个默认的dist计算dist * dist时，结果很可能会溢出成inf（无穷大）。理论上sqrt(DBL_MAX)的平方等于DBL_MAX，但双精度浮点数只有52位尾数精度，实际计算时没法精确还原DBL_MAX，反而会因为舍入溢出到无穷大。

一旦出现inf，比如你算法里有if (点的平方距离 <= dist²)的判断，任何有限的平方距离和inf比较都会永远为真，这完全违背了你“取尽可能大的有效距离”的初衷。

2. 改进写法std::sqrt(std::numeric_limits<double>::max() - 100)是否靠谱？ #

这个写法确实比原写法安全，但100这个数有点太随意了：

• 减去100后，sqrt(DBL_MAX - 100)的平方就是DBL_MAX - 100，这个值肯定小于DBL_MAX，不会溢出成inf，算法里的平方运算能得到一个有效的最大有限值。
• 不过更严谨的做法是用std::nextafter(std::numeric_limits<double>::max(), 0.0)，也就是比DBL_MAX小一点点的“次大值”，取它的平方根，这样平方后刚好是最接近DBL_MAX的有限值，比硬减100更合理。

3. 更聪明的替代方案 #

既然你的算法实际用的是平方距离，为什么不直接默认平方距离的最大值呢？这样能省掉平方根运算的开销，还彻底消除溢出风险：

Result foo(double squared_dist = std::numeric_limits<double>::max()) const;

然后如果用户习惯传原始距离，你再重载一个函数做转换：

Result foo(double dist) const {
    return foo(dist * dist);
}

这样逻辑更清晰，还避开了默认参数里的浮点运算坑，一举两得。

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