17

看看这段代码：

#include <cmath>
#include <iostream>
using namespace std;
class Sphere
{
    double r;
public:
    double V() const { return (4/3) * 3.14 * pow(r,3); }
    bool equal(const Sphere& s) const
    {
        cout  << V() << " == " << s.V() << " : " << ( V() == s.V() );
        return ( V() == s.V() );

    }

    explicit Sphere(double rr = 1): r(rr){}

};
main()
{
    Sphere s(3);
    s.equal(s);
}

输出84.78 == 84.78 : 0意味着相同的方法不会每次都返回相同的值，即使所有参数都是静态的？

但是如果我写3.0而不是3.14在V()方法定义中，像这样：

double V() const { return (4/3) * 3.0 * pow(r,3); }

然后，输出是：84.78 == 84.78 : 1

这里发生了什么？对于我的程序，我需要这种方法，它将比较两个对象的体积，但这是不可能的？我敲了很长时间才弄清楚问题的原因，幸运的是我找到了，但现在我不明白为什么？它与编译器（GCC）有关还是我在这里遗漏了一些重要的东西？

4

2 回答 2

34

使用运算符比较浮点值==很容易出错；应该相等的两个值可能不是由于算术舍入误差。比较这些的常用方法是使用 epsilon：

bool double_equals(double a, double b, double epsilon = 0.001)
{
    return std::abs(a - b) < epsilon;
}
于 2013-09-24T07:05:37.467 回答
4

浮点比较有两个问题：

(1) 浮点运算通常至少涉及难以预测的微小舍入误差。因此，在数学上应该给出相同结果的两个浮点运算（如 4.7 * (1.0 / 3.14) 与 4.7 / 3.14）可能会给出不同的结果。

(2) 有时允许编译器以高于所需的精度进行浮点运算。还允许以其他时间所需的精度执行完全相同的浮点运算。因此，完全相同的操作可能会产生略微不同的结果，这就是您在此处看到的。

为了解决OP的问题，这看起来是由（2）引起的。我会尝试查找是否有任何编译器选项可以阻止编译器使用比需要更高的精度。

于 2014-07-10T00:40:35.540 回答