当结果对称时，可以优化 R 矩阵乘法吗？

Question

我需要获取特定矩阵中每一对列的内积，我通过计算来实现

t(M) %*% M

但是，这自然会产生对称的结果，只需完成两倍以上的必要工作（我也不需要对角线）。显然我可以将乘法分解为单独的内积运算，但有没有更好的方法来计算乘积的上三角部分？

Answer 1

来自help("crossprod")中的描述：

给定矩阵 x 和 y 作为参数，返回一个矩阵叉积。 这在形式上等同于（但通常略快于） 致电 t(x) %*% y (crossprod) 或 x %*% t(y) (tcrossprod)。

因此，请使用crossprod(M)。