这个矩阵肯定是稀疏的,S有M^2-M非零项,我的块矩阵将有N(M^2-M) out (NM)^2或关于1/N%非零条目,但是我将把它添加到另一个NMxNM矩阵中,而我并不期望它是稀疏的。,既然我要把我的块矩阵添加到一个完整的矩阵中,那么用“稀疏”的方式来写我的代码会有什么速度上的提高吗?我不停地重复,
我想要为x解决以下线性系统其中a是稀疏的,b是正则列矩阵。但是,当我插入通常的np.linalg.solve(A,b)例程时,它给我一个错误。然而,当我使用np.linalg.solve(A.todense(),b)时,它工作得很好。如何在保持A的sparseness的情况下使用这种线性求解呢?原因是关于150 x 150的A相当大,大约有50个这样的