双倍精度算法详解：在不离开FPU的情况下实现31位精度的算术技巧

本文深入探讨了一种名为“双倍精度”的算术技术，旨在填补标准双精度浮点数与任意精度数学库之间的性能与精度空白。该方法通过将两个标准的64位双精度浮点数组合并在逻辑上视为单一数值，从而在不依赖堆分配、外部库或原生二进制文件的情况下，提供约31位十进制精度（约106位有效二进制位）。文章详细解析了实现该技术的核心算法，即基于Knuth的“two-sum”和利用硬件FMA（融合乘加）指令实现的“two-product”，这些误差自由变换确保了运算的精确性。性能基准测试表明，在实际应用中，双倍精度的计算成本约为普通双精度浮点数的9倍，但在同等精度水平下，其性能却比使用MPFR的任意精度库快了近9倍（整体比普通double慢约81倍）。这种技术特别适用于需要深度缩放或高精度中间计算的场景，如分形渲染，且由于其基于值类型的实现，支持向量化并避免了垃圾回收的开销。

事件分析

核心观点：通过算法巧妙利用FMA指令，软件构造的高精度数值类型能以极低成本填补硬件浮点数与重型数学库之间的巨大性能鸿沟。