我需要在一个列表中solve超过上千个矩阵。但是,我得到了错误Lapack routine dgesv: system is exactly singular。这里是我问题的一个基本例子(A是经过一些计算后的矩阵,R是我需要做的下一个计算):R = solve(diag(length(p))-A)的函数
我有一个大小为1393 x 1393的邻接矩阵,它看起来像:
我正在做基于这个矩阵的空间建模,以后会有大量的矩阵逆。我听说这样的稀疏矩阵,我们可以重新排序行和列,以达到最小的带宽,这将大大减少矩阵逆所需的时间。我发现的算法叫做Gibbs-Poole-Stockmeyer,我找到了一个FORTRAN 。我是一个R用户,目前正