python - 测试两个 float 是否在 1 ULP 之内

Question

如何判断两个 float 之差是否最多为1 ULP？

类似于下面的代码，它(忽略符号和指数，并且)提取尾数的二进制，但稍微不那么糟糕:

a = int(<float>.hex().lstrip('-').lstrip('0x')[1:].lstrip('.')[:-3].rstrip('+').rstrip('p'), 16)
b = int(<float>.hex().lstrip('-').lstrip('0x')[1:].lstrip('.')[:-3].rstrip('+').rstrip('p'), 16)
print abs(a - b) <= 1

Answer 1

使用 NumPy，您可以获得给定数字之后的下一个 float :

# Test whether taking the smallest possible step from b in the direction of a gets you a.
# Special case: If b == a, then np.nextafter(b, a) == a
a == np.nextafter(b, a)

(如果两个数字具有不同的 ULP，或者其中一个数字是无穷大，这可能不是您想要的。)