Perl中正向与反向累加浮点数组结果不一致问题排查

你碰到的这个情况完全是浮点数运算的正常精度特性导致的，和代码逻辑无关，放心，你的代码没有bug。

为什么会出现差异？ #

计算机中的双精度浮点数（Perl默认使用的浮点类型）是用二进制来近似表示十进制小数的，很多日常用到的十进制小数（比如0.1）根本无法被二进制精确存储，每次进行加减运算时都会产生极其微小的舍入误差。

当你改变累加顺序时，这些微小误差的累积路径就变了：正向累加是从数组第一个元素开始逐步加，反向累加是从最后一个元素倒着加，每次相加的数值组合不同，舍入的方向和幅度也会有细微差别，最终就导致了结果出现像-2.27374e-13这种极小的差异——这个量级的误差远低于双精度浮点数15-17位有效数字的精度范围，属于完全可接受的正常现象。

代码没问题，是比较方式错了 #

你的两种累加逻辑都是正确的：$sum += $_ for @x;和反向循环累加的代码没有逻辑漏洞。真正的问题出在你用!=直接比较两个浮点数是否相等——因为浮点数的精度特性，几乎所有情况下，理论上相等的浮点计算结果都会存在这种极小的误差，直接用!=判断肯定会触发“不相等”的提示。

正确的做法是比较两个值的绝对差是否小于一个极小的阈值，比如：

my $epsilon = 1e-9; # 可以根据你的业务场景调整这个阈值
if (abs($total - $sum) > $epsilon) {
    printf("SMA($what, $periods) total ($total) != sum ($sum) %g\n", ($total - $sum));
}

这样就能过滤掉这种正常的精度差异，只在真正出现计算错误时才触发提示。

要不要追求更高精度？ #

如果你的场景对精度要求极高（比如某些高精度科学计算），可以考虑使用Perl的Math::BigFloat模块来进行高精度运算，但从你的变量名closing、SMA来看是金融场景，双精度浮点数的精度已经完全够用，没必要额外引入模块增加复杂度。

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