AVX512中vrndscalepd无法实现保留两位小数的正确舍入问题

vrndscalepd无法直接满足你的需求——它的设计目标是将数值舍入到2的负整数次幂的倍数（如0.5、0.25、0.125等），而0.01（1/100）并非2的幂次，因此所有模式都会得到0.25倍数的结果，无法精准保留两位小数。

正确实现步骤 #

要实现保留两位小数的舍入，需要通过「缩放→舍入整数→还原」的流程，利用AVX512的通用浮点运算和舍入指令完成：

1. 放大数值：将所有数乘以100，把需要保留的两位小数转化为整数部分；
2. 舍入到整数：使用vroundpd指令选择合适的舍入模式（如最近偶数、四舍五入等）；
3. 还原数值：将结果除以100（或乘以0.01，更高效），恢复到两位小数格式。

代码示例 #

假设原始数据存储在zmm8，掩码k1用于部分元素操作：

方案1：最近偶数舍入（IEEE标准舍入，与常规四舍五入几乎一致）

; 提前定义常量（可放在数据段）
one_hundred dq 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0, 100.0
point_zero_one dq 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01

; 执行舍入流程
vmulpd zmm10, zmm8, zmmword ptr [one_hundred]   ; 乘以100，放大数值
vroundpd zmm9{k1}{z}, zmm10, 0b0001             ; 舍入到最近偶数（低4位为0001）
vmulpd zmm9{k1}{z}, zmm9, zmmword ptr [point_zero_one] ; 乘以0.01还原

方案2：传统四舍五入（0.5始终向上舍入）

如果需要严格遵循传统四舍五入规则（小数部分≥0.5时向上舍入），可以先给放大后的数值加0.5再截断：

point_five dq 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5

vmulpd zmm10, zmm8, zmmword ptr [one_hundred]   ; 乘以100
vaddpd zmm10, zmm10, zmmword ptr [point_five]  ; 加0.5，为截断做准备
vroundpd zmm9{k1}{z}, zmm10, 0b0000             ; 截断到整数（低4位为0000）
vmulpd zmm9{k1}{z}, zmm9, zmmword ptr [point_zero_one] ; 还原为两位小数

验证你的测试数据 #

以第一个数值3.9304247359239284为例：

• 乘以100得到393.04247359239284
• 方案1：最近偶数舍入后为393，还原后是3.93
• 方案2：加0.5后为393.54247359239284，截断后为393，还原后同样是3.93

其他数值也会得到符合预期的两位小数结果，比如0.85804618270500566处理后为0.86，9.1150061014624413处理后为9.12。

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