相同64位双精度浮点数为何出现不同内存位模式？

你遇到的核心问题是：两个显示为1020.0的双精度浮点数，实际内存位模式存在微小差异，导致结构体句柄判定不一致。

为什么会出现这种情况？ #

从给出的十六进制位模式来看：

• 第一个数的位模式是0x408FE00000000000，这是1020.0的精确双精度表示。
• 第二个数的位模式是0x408FDFFFFFFF，这个值实际是略小于1020.0的浮点数，只是调试器在显示时做了四舍五入，才呈现为1020.00000000000。

导致这种差异的常见原因有：

• 计算过程的精度损失：其中一个数可能不是直接赋值1020.0，而是通过浮点运算得到（比如看似简单的1000.0 + 20.0，如果中间涉及其他精度受限的运算，可能产生微小误差）。
• 编译器优化与扩展精度：Intel编译器默认可能使用x87的80位扩展精度寄存器进行计算，当把寄存器中的值存储回64位双精度内存时，会截断低位，可能产生和直接赋值不同的位模式。
• 初始化方式差异：一个结构体成员是直接初始化1020.0，另一个是通过其他带微小误差的变量赋值。

解决办法 #

1. 放弃位模式比较浮点数：浮点数的本质是近似表示（即使是整数，也可能因为运算或存储方式产生微小差异），直接比较位模式判定结构体相等是不可靠的。
2. 使用容差比较：对浮点数成员，判断两个值的绝对差或相对差是否小于一个合理阈值，比如：

   ABS(a%wd - b%wd) < 1e-9
   

   阈值可以根据你的精度需求调整。
3. 强制严格浮点语义：使用Intel Fortran的编译选项-fp-model strict（Linux/macOS）或/fp:strict（Windows），强制编译器严格遵循浮点数标准，避免因扩展精度或优化导致的存储差异。
4. 统一初始化方式：确保所有结构体的浮点数成员都通过相同的方式初始化（比如直接赋值常量，而非运算结果），减少差异来源。

内容的提问来源于stack exchange，提问作者Arthur