我是 C++ 新手。在阅读 C++ 时，我有一个疑问。 如何判断两个浮点数是否相等？
提前致谢
【问题讨论】：
标签： c++ floating-point
您需要了解一个特殊的常量，称为 DBL_EPSILON（或 FLT_EPSILON）。这是可以添加到 1.0 并更改其值的最小值。 1.0 的值非常重要——添加到 DBL_EPSILON 时较大的数字不会改变。现在，您可以将此值缩放到您要比较的数字，以判断它们是否不同。比较两个双精度的正确表达式是：
if (fabs(a-b) <= DBL_EPSILON * fmax(fabs(a), fabs(b)))
{
// ...
}
【讨论】：
FLT_EPSILON 是不是可以添加到 1.0f 并更改其值的最小值。 blog.frama-c.com/index.php?post/2013/05/09/FLT_EPSILON
如果您的浮点类型使用 IEEE 754 表示（很可能是这种情况），那么您应该使用浮点的 二进制表示 的顺序与 按价值排序。也就是说，如果将浮点数的二进制表示增加一位，则会得到下一个更大的数字。
利用这个事实，我们可以通过计算它们的二进制差异来比较浮点数。这称为“按最后一个单元 (ULP) 进行比较”。有一些涉及符号、零、无穷大和 NaN 的微妙之处，但这就是它的要点。这是comprehensive article 解释这一点。
基本上，我们认为两个浮点数相等，如果它们在最后一个位置的少数几个单位不同。结合编译器关于其数学函数精度的文档和您自己的代码，您可以确定哪个截止值适合您的需求。
在伪代码中：
double x, y;
// this is type punning, should be done differently in reality
uint64_t ux = *reinterpret_cast<const uint64_t*>(&x);
uint64_t uy = *reinterpret_cast<const uint64_t*>(&y);
return abs(ux - uy) < CUT_OFF; // e.g. CUT_OFF = 3;
上面的代码只是一个粗略的例子，它不会工作，在最后的比较之前你必须处理很多特殊情况。有关详细信息，请参阅文章。
【讨论】：
if (sign(a) != sign(b)) return a == b开始，所以我们以后可能不必担心这个。
显然，您不应该使用operator == 来比较它们。
这里的重要概念是，如果您的两个浮点数之差足够小到您要解决的问题的精度要求或小于您的误差范围，我们应该将它们视为相等。
有一些实用的方法建议如
fabs(f1 - f2) < precision-requirement
fabs(f1 - f2) < max(fabs(f1), fabs(f2)) * percentage-precision-requirement
【讨论】：
== 是错误的解决方案，而是根据情况权衡它。当然，这需要更多地了解不精确的浮点表示。
precision_requirement。