IEEE 754浮点除法的可逆性是什么？我的意思是由标准保证如果double y = 1.0 / x然后x == 1.0 / y，即x可以精确地一点一点恢复？y是infinity或NaN的情况是明显的异常(exception)。
最佳答案
是的，IEEE 754 double (*)值x就是x != 1.0 / (1.0 / x)。
用手可以很容易地用此属性构建一个正常值的示例:在C99's hexadecimal notation for floating-point values中写了0x1.fffffffffffffp0的那个就是1.0 / (1.0 / 0x1.fffffffffffffp0) == 0x1.ffffffffffffep0。很自然地期望0x1.fffffffffffffp0是一个反例，因为1.0 / 0x1.fffffffffffffp0落在binade的开头，而binade的浮点数密度较小，因此最里面的除法必须发生较大的相对误差。更准确地说，1.0 / 0x1.fffffffffffffp0恰好位于0.5及其 double 后继字符之间的中点上方，因此1.0 / 0x1.fffffffffffffp0舍入为0.5的后继字符，并且具有较大的相对误差。
采用十进制%.16e格式，0x1.fffffffffffffp0是1.9999999999999998e+00，0x1.ffffffffffffep0是1.9999999999999996e+00。
(*)对于任何IEEE 754格式，没有理由使反函数具有该属性