R稀疏矩阵功率

Question

我使用Matrix包创建一个包含大量零的大型（~14000x14000）稀疏矩阵。有谁知道计算这个矩阵功效的最佳方法？

我试过A_pow2 = A％^％2但是我得到了错误：A％^％2中的错误：不是矩阵。这是一个返回相同错误的简单示例：

A = matrix(3,2,2)
A = Matrix(A,sparse=TRUE)
Apow2 = A%^%2

Answer 1

（感谢@ Roland的评论）

自定义功能可能能够解决您的问题。根据{{​​1}}

的文档

  

计算矩阵的k次幂。而x ^ k计算元素   功率，x％^％k对应于k-1矩阵乘法，x％*％x   ％*％...％*％x。

我们可以写一个新的中缀运算符来执行乘法k-1次。不确定它的扩展程度，但它适用于较小的例子。

?expm::`%^%`