A = sparse.csc_matrix([[1., 5., 0.], [5., 0., 2.], [0., 2., 0.]], dtype=float) B = linalg.spilu(A) #不完全LU分解 b = np.array([1., 2., 3.], dtype=float) x=B.solve(b) # solve函數求解Ax=b

B的類型是SuprLU

scipy.sparse.linalg.SuperLU — SciPy v1.7.1 Manual

Sparse linear algebra (scipy.sparse.linalg) — SciPy v1.7.1 Manual

循環矩陣的生成和求解

scipy.linalg.circulant — SciPy v1.7.1 Manual

scipy.linalg.solve_circulant — SciPy v1.7.1 Manual

總結

以上是生活随笔為你收集整理的Python 不完全LU分解的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。