R语言大尺寸矩阵乘向量报non-conformable arrays错误求解

报错核心原因 #

• R中*为逐元素乘法，运算时会将矩阵按列优先规则展开为长向量，再基于R的循环补齐规则对两个向量逐元素运算，要求短向量长度可被长向量长度整除。
• 小样本测试时100×100矩阵展开后长度为10000，刚好是长度100的test.scores的整数倍，符合循环补齐要求，因此无报错。
• 37000×37000矩阵展开后总元素数达13.69亿，两种情况会触发报错：一是32位R版本或内存不足时，运算前的预校验逻辑会误判为维度不兼容；二是scores被隐形设置了矩阵维度属性，不符合逐元素运算的长度要求。

最优实现方案 #

无需先生成单位矩阵再做乘法，直接调用diag()函数传入scores向量即可一步生成目标对角矩阵，运算效率是原有写法的数十倍，也不会触发维度兼容问题：

# 直接生成对角线为scores取值的对角矩阵
mat <- diag(scores)

低内存场景优化方案 #

37000×37000的稠密矩阵需要约11GB内存存储，如果设备内存不足，可以使用稀疏对角矩阵结构，仅存储非零的对角线元素，内存占用不到1MB：

library(Matrix)
mat <- Diagonal(x = scores)

备选修改方案 #

如果确实需要基于已生成的单位矩阵调整，直接对对角线赋值即可，性能远高于逐元素乘法：

mat <- diag(length(scores))
diag(mat) <- scores

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