我目前正在开发一个需要将某些稀疏矩阵相乘的项目。
使用Matrix
包在R中工作我通过“对角线”或简称“矩阵”之类的命令创建矩阵,例如
A <- Diagonal(10)
B <- kronecker(Matrix(diag(10)),t((rep(1,10))))
C <- A%*%B
矩阵A
属于dgTMatrix
类,而矩阵B
属于类ddiMatrix
。
有没有办法加快乘法的计算?
我正在考虑RcppArmadillo
,但我不知道如何处理这些类,如果它甚至表现更好。