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

Question

我需要计算特定矩阵中每对列的内积，这是通过计算实现的

t(M) %*% M

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

Answer 1

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

Given matrices x and y as arguments, return a matrix cross-product. This is formally equivalent to (but usually slightly faster than) the call t(x) %*% y (crossprod) or x %*% t(y) (tcrossprod).

因此，请使用crossprod(M)。