如何使用SuperLU对象得到矩阵的逆?

时间:2017-03-21 03:40:08

标签: python numpy matrix scipy linear

我试图从SuperLU对象获得矩阵逆。 我尝试过以下方法:

lu = sp.linalg.splu(A)
I = np.eye(A.shape[0])
invA = lu.solve(I)
x = invA.dot(b)

我得到的x与执行x = sp.linalg.solve(b)不同。我也试过做invA.dot(A),结果矩阵单位矩阵。

0 个答案:

没有答案