Kat*_*ine 7 r matrix-inverse
计算对称密集矩阵(2000 * 2000)的逆对角线的最佳方法是什么？目前我首先计算逆solve(x),然后提取对角线(diag(y)).即使它有效,但我想知道是否有更好的方法来做到这一点,所以代码运行得更快.我尝试过,chol2inv()但由于我的矩阵不是正确的,所以它没有用.
2000 * 2000
solve(x)
diag(y)
chol2inv()
更新:对于任何可能感兴趣的人,我都可以通过使用优化的数学库英特尔MKL来加速矩阵求逆.在我的机器上反转2000*2000矩阵需要3秒钟.英特尔MKL适用于Microsoft R Open.
小智 0
如果你的矩阵没有好的属性，比如对称、对角或正定，那么遗憾的是你唯一的选择是sum(diag(solve(x)))
sum(diag(solve(x)))
在你的矩阵上运行需要多长时间？
归档时间：
8 年，4 月 前
查看次数：
545 次
最近记录：