如何对80位ieee浮点进行相等测试?

既然d内置了80位real(据我所知),为什么不直接使用与epsilon值比较的标准方法呢?如果您事先知道大致范围(如美元),则这可能是一个固定值:

if (abs (a - b) < 1e-6) // effectively equal

或可接受的相对误差,如百万分之一的平均值:

if (abs (a - b) < ((a + b) / 2) / 1e6) // effectively equal

记住,我不知道D,上面的代码只是为了演示。

我现在的解决方案是

bool Near(real a, real b, int count = 5)
{
    // Returns the number of mantissa bits which are equal in x and y.
    int i = std.math.feqrel!(real)(a,b);
    return i + count >= real.mant_dig;
}

它给出输入中未匹配的位数。我不确定这在2次方附近会有多好。