实现 FMA 却揪出 C 与 Rust 标准库的隐藏 bug

概述

作者为 fearless_simd 实现软件版 FMA（fused multiply-add，一次舍入算出 a*b+c）时，循着一篇 2008 年用 Coq 形式化验证过的论文，用 f64 算完再舍入回 f32。正当他以为大功告成，一位陌生网友贴出一组让他的实现与硬件结果分道扬镳的输入——问题出在他漏掉了对亚正常数（subnormal）的特殊处理。深入排查后他发现：Rust 标准库内置的 f32::mul_add 与 std::simd 存在完全相同的 bug，而它又源自 musl libc 的 fmaf.c，后者版权标注又追溯到 FreeBSD。这段带着 bug 的代码在标准库与工具链里被复制粘贴了二十年。作者逐一修复并提交 PR，还顺带发现 32 位 ARM 上 nightly 的 f128 类型因 rustc ABI bug 彻底错乱。他提醒：在没有硬件 FMA 的旧 Intel 与国产 Hygon 芯片上，这个错误舍入虽只差一个最低有效位，却可能被放大，让确定性模拟在不同机器上分道扬镳。