本文简要介绍 python 语言中 scipy.sparse.linalg.SuperLU
的用法。
用法:
class scipy.sparse.linalg.SuperLU#
稀疏矩阵的 LU 分解。
因式分解表示为:
Pr @ A @ Pc = L @ U
要构造这些
SuperLU
对象,请调用splu
和spilu
函数。注意:
例子:
LU分解可用于求解矩阵方程。考虑:
>>> import numpy as np >>> from scipy.sparse import csc_matrix, linalg as sla >>> A = csc_matrix([[1,2,0,4],[1,0,0,1],[1,0,2,1],[2,2,1,0.]])
对于给定的右侧可以解决这个问题:
>>> lu = sla.splu(A) >>> b = np.array([1, 2, 3, 4]) >>> x = lu.solve(b) >>> A.dot(x) array([ 1., 2., 3., 4.])
lu
对象还包含分解的显式表示。排列表示为索引的映射:>>> lu.perm_r array([0, 2, 1, 3], dtype=int32) >>> lu.perm_c array([2, 0, 1, 3], dtype=int32)
L 和 U 因子是 CSC 格式的稀疏矩阵:
>>> lu.L.A array([[ 1. , 0. , 0. , 0. ], [ 0. , 1. , 0. , 0. ], [ 0. , 0. , 1. , 0. ], [ 1. , 0.5, 0.5, 1. ]]) >>> lu.U.A array([[ 2., 0., 1., 4.], [ 0., 2., 1., 1.], [ 0., 0., 1., 1.], [ 0., 0., 0., -5.]])
可以构造置换矩阵:
>>> Pr = csc_matrix((np.ones(4), (lu.perm_r, np.arange(4)))) >>> Pc = csc_matrix((np.ones(4), (np.arange(4), lu.perm_c)))
我们可以重新组装原始矩阵:
>>> (Pr.T @ (lu.L @ lu.U) @ Pc.T).A array([[ 1., 2., 0., 4.], [ 1., 0., 0., 1.], [ 1., 0., 2., 1.], [ 2., 2., 1., 0.]])
shape
原始矩阵的形状为整数元组。
nnz
矩阵中非零元素的数量。
- scipy.sparse.linalg.SuperLU.perm_c
置换 Pc 表示为索引数组。
- scipy.sparse.linalg.SuperLU.perm_r
置换 Pr 表示为索引数组。
L
单位对角线为
scipy.sparse.csc_matrix
的下三角因子。U
上三角因子为
scipy.sparse.csc_matrix
。
属性 ::
相关用法
- Python SciPy linalg.eigvalsh_tridiagonal用法及代码示例
- Python SciPy linalg.cdf2rdf用法及代码示例
- Python SciPy linalg.LaplacianNd用法及代码示例
- Python SciPy linalg.solve_circulant用法及代码示例
- Python SciPy linalg.polar用法及代码示例
- Python SciPy linalg.clarkson_woodruff_transform用法及代码示例
- Python SciPy linalg.rsf2csf用法及代码示例
- Python SciPy linalg.hessenberg用法及代码示例
- Python SciPy linalg.tril用法及代码示例
- Python SciPy linalg.triu用法及代码示例
- Python SciPy linalg.svd用法及代码示例
- Python SciPy linalg.ishermitian用法及代码示例
- Python SciPy linalg.invhilbert用法及代码示例
- Python SciPy linalg.factorized用法及代码示例
- Python SciPy linalg.lu_factor用法及代码示例
- Python SciPy linalg.lsqr用法及代码示例
- Python SciPy linalg.cho_factor用法及代码示例
- Python SciPy linalg.fractional_matrix_power用法及代码示例
- Python SciPy linalg.eig_banded用法及代码示例
- Python SciPy linalg.tanhm用法及代码示例
- Python SciPy linalg.orthogonal_procrustes用法及代码示例
- Python SciPy linalg.use_solver用法及代码示例
- Python SciPy linalg.qr_multiply用法及代码示例
- Python SciPy linalg.spsolve用法及代码示例
- Python SciPy linalg.LinAlgError用法及代码示例
注:本文由纯净天空筛选整理自scipy.org大神的英文原创作品 scipy.sparse.linalg.SuperLU。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。