Perl PDL是否有等效Math::Round::nearest()的向量化取整方法？

PDL核心库没有直接提供与Math::Round::nearest()功能完全等价的内置方法，但完全可以通过原生向量化操作实现同等功能，不需要做unpdl/re-pdl格式转换，没有额外运行开销。

原有实现的问题 #

你之前使用rint($pdl/$n)*$n的写法无法得到正确结果，核心原因有两个：

• rint()采用「四舍六入五成双」的银行家舍入规则，遇到恰好处于两个整数中点的值时，会取相邻的偶数结果，和Math::Round::nearest()「中点值向远离0方向进位」的规则不一致
• 二进制浮点数本身存在存储精度误差，部分十进制小数无法用二进制浮点数精确表示，会导致除法计算结果出现微小偏差，进一步影响舍入结果

原生向量化实现代码 #

Math::Round::nearest()的核心逻辑是：返回离目标值最近的指定精度的整数倍，当目标值恰好落在两个倍数的中点时，取绝对值更大的结果。对应的PDL向量化实现如下：

use PDL;
use PDL::Math; # 仅依赖内置的floor函数

sub pdl_nearest {
    my ($val, $precision) = @_;
    my $quotient = $val / $precision;
    # 加入极小容差抵消浮点数存储误差，可根据业务精度需求调整
    $quotient += 1e-10 * sign($quotient);
    my $rounded = floor($quotient + 0.5);
    return $rounded * $precision;
}

# 测试用例
my $pdl = pdl [4.45, 5.55, 45];
my $n = pdl [0.1, 0.2, 10]; # 注：原示例中第二个精度值若为0.3，5.55的中点舍入结果应为5.7，若预期为5.6则对应精度是0.2
print pdl_nearest($pdl, $n);
# 输出 [4.5 5.6 50]，和预期结果一致

实现说明 #

• 该实现完全支持PDL的向量化广播规则，既支持标量精度（所有元素按同一精度舍入），也支持和输入数组同维度的精度向量（每个元素按对应精度舍入），性能和PDL内置函数无差异
• 代码中加入的1e-10量级容差是为了抵消二进制浮点数的存储误差，不会对正常舍入逻辑产生影响，可根据实际业务的精度要求调整大小
• 如果你的使用场景只涉及非负数，可以进一步简化代码为floor($pdl/$n + 0.5 + 1e-10)*$n，不需要额外处理符号逻辑

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