标签: python sparse-matrix linear-algebra
此刻,我正试图提高程序效率。我必须解决一个线性系统Ax = b。但是A有很多零条目,因此我使用scipy.sparse.lil_matrix来消耗更少的内存,因为A是80k x 80k矩阵或更多。
scipy.sparse.lil_matrix
现在,由于我只需要解决方案的第一要素,所以我试图更有效地解决系统问题。到目前为止,我已经使用了以下内容:np.linalg.solve(A, b)[0]或sp.sparse.linalg.spsolve(A, b)[0]
np.linalg.solve(A, b)[0]
sp.sparse.linalg.spsolve(A, b)[0]
还有其他解决方案吗?我曾考虑过应用克莱默法则,但这需要行列式。