浮动和双重比较最有效的方法是什么？

Question

比较两个double 或两个float 值的最有效方法是什么？

仅仅这样做是不正确的：

bool CompareDoubles1 (double A, double B)
{
   return A == B;
}

但类似：

bool CompareDoubles2 (double A, double B) 
{
   diff = A - B;
   return (diff < EPSILON) && (-diff < EPSILON);
}

似乎浪费了处理。

有人知道更智能的浮点比较器吗？

Answer 1

如需更深入的方法，请阅读Comparing floating point numbers。这是该链接中的代码 sn-p：

// Usable AlmostEqual function    
bool AlmostEqual2sComplement(float A, float B, int maxUlps)    
{    
    // Make sure maxUlps is non-negative and small enough that the    
    // default NAN won't compare as equal to anything.    
    assert(maxUlps > 0 && maxUlps < 4 * 1024 * 1024);    
    int aInt = *(int*)&A;    
    // Make aInt lexicographically ordered as a twos-complement int    
    if (aInt < 0)    
        aInt = 0x80000000 - aInt;    
    // Make bInt lexicographically ordered as a twos-complement int    
    int bInt = *(int*)&B;    
    if (bInt < 0)    
        bInt = 0x80000000 - bInt;    
    int intDiff = abs(aInt - bInt);    
    if (intDiff <= maxUlps)    
        return true;    
    return false;    
}

Answer 2

你写的代码有问题：

return (diff < EPSILON) && (-diff > EPSILON);

正确的代码是：

return (diff < EPSILON) && (diff > -EPSILON);

（......是的，这是不同的）

我想知道在某些情况下晶圆厂是否不会让你失去惰性评估。我会说这取决于编译器。您可能想尝试两者。如果它们平均相等，请使用 fabs 实施。

如果您有一些关于两个浮点数中哪一个更可能大于另一个的信息，您可以按照比较的顺序进行操作，以更好地利用惰性求值。

最后你可能会通过内联这个函数得到更好的结果。不过不太可能改善很多...

编辑：OJ，感谢您更正您的代码。我相应地删除了我的评论