numpy and PyTorch may compute matrix multiplications using different BLAS implementations with different accumulation orders. Floating-point addition is not associative, so $(a + b) + c \neq a + (b + c)$ in general. (2) **Fused operations** — PyTorch may use fused multiply-add (FMA) instructions tha