为何9元素NumPy数组除法触发RuntimeWarning溢出，其他长度则无？

问题现象 #

用NumPy做数组除法时碰到个诡异情况：

• 创建全0的float64数组，执行数组 / np.float64(1e-309)操作，当数组长度为1、4、7、8、12、16、100、200时，无警告正常运行；
• 但数组长度为9时，会触发RuntimeWarning: overflow encountered in divide警告；
• 总结规律：只要数组元素数超过8且无法被4整除，就会出现该警告。

原因解析 #

这是NumPy底层SIMD（单指令多数据）硬件优化导致的逻辑差异：
NumPy处理float64类型数组运算时，默认会用128位SIMD指令，一次批量处理4个元素：

• 若数组长度是4的倍数，所有元素都走批量运算逻辑，此时0除以极小值的操作不会触发溢出警告；
• 若数组长度不是4的倍数，剩余的元素会走单元素运算分支，这个分支的溢出检测逻辑更严格——1e-309已经接近float64能表示的最小正数（约2.2e-308），0除以它会被判定为溢出，从而弹出警告。

代码复现 #

import numpy as np

# 无警告正常运行
arr_normal = np.zeros(8, dtype=np.float64)
arr_normal / np.float64(1e-309)

# 触发溢出警告
arr_warning = np.zeros(9, dtype=np.float64)
arr_warning / np.float64(1e-309)

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