如何判断double类型变量是否为零（规避浮点数异常行为）

这可能是个老生常谈的问题，但我没找到满意的解决方案。我需要判断浮点运算的结果是否为零，同时避开所有浮点数的异常行为。我试了下面的代码但完全失败了，求帮忙！

#include <iostream>
#include <cmath>
#include <limits>
#include <iomanip>
#include <iostream>

using std::cin;
using std::cout;
using std::endl;

int main()
{
    double d1 = 0.2;
    double d2 = 1.0;
    double d3 = 0.3;

    d1 = d1 + 0.1;
    d1 = d1 / d3;
    d1 = d1 - d2;

    // print d1 with max precision
    cout << std::setprecision(std::numeric_limits<double>::max_digits10) << d1 << endl;

    // accurately compare d1 with zero with regard to precision
    if (std::abs(d1) < std::numeric_limits<double>::epsilon())
    {
        cout << "d1 is zero" << endl;
    }
    else
    {
        cout << "d1 is not zero" << endl;
    }

    return 0;
}

问题分析 #

你代码里用std::numeric_limits<double>::epsilon()判断是否为零，这是典型的错误用法。epsilon是double类型能表示的最小非零值与1.0的差值，它仅适用于判断两个数是否“近似相等”（基于1.0的尺度），但完全不适合判断一个数是否接近零。

你的运算过程中：

1. 0.2 + 0.1无法得到精确的0.3（二进制浮点数无法精确表示0.1、0.2这类十进制小数）
2. 用这个近似值除以0.3的近似值，结果接近1.0但并非精确值
3. 减去1.0后得到的差值绝对值远大于epsilon，导致判断条件失效

解决方案 #

不要用epsilon，而是根据业务场景定义合理的误差阈值，比如针对小数运算场景，可选用1e-9或更小的数值（具体取决于你需要的精度）。

修改后的代码示例：

#include <iostream>
#include <cmath>
#include <limits>
#include <iomanip>

using std::cin;
using std::cout;
using std::endl;

int main()
{
    double d1 = 0.2;
    double d2 = 1.0;
    double d3 = 0.3;

    d1 = d1 + 0.1;
    d1 = d1 / d3;
    d1 = d1 - d2;

    // 以最大精度打印d1
    cout << std::setprecision(std::numeric_limits<double>::max_digits10) << d1 << endl;

    // 自定义误差阈值判断是否近似为零
    const double tolerance = 1e-9;
    if (std::abs(d1) < tolerance)
    {
        cout << "d1 近似为零" << endl;
    }
    else
    {
        cout << "d1 不为零" << endl;
    }

    return 0;
}

额外说明 #

• epsilon的适用场景是：判断abs(a - b) <= epsilon * abs(a)这类相对误差，而非判断绝对数值是否接近零
• 阈值的选择要匹配业务需求：金融计算可能需要更小的阈值（如1e-12），图形学等场景1e-6或1e-9即可满足要求

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