正如我在标题中提到的，到目前为止，我使用以下代码为R中的矩阵提供动力。此方法使用expm包及其%^%运算符。
Qmp <- expm::`%^%`( as.matrix(Qm), (as.numeric(L)-1) )
最后一天，我正在尝试为尺寸为
的大型矩阵执行此操作> dim(Qm)
[1] 17328 17328
我意识到这需要很长时间才能完成。老实说，我正在运行该代码一小时但仍未完成。
所以我想知道是否还有其他办法来处理这样的任务并且更快。
编辑：
Qm上电之后我也在执行并且两个矩阵相乘，如
Pm <- Qmp %*% as.matrix(Rm)
其中
dim(Rm)
>17328 58