标签: 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),结果矩阵不单位矩阵。
x
x = sp.linalg.solve(b)
invA.dot(A)