naginterfaces.library.lapacklin.zspsv¶
- naginterfaces.library.lapacklin.zspsv(uplo, ap, b)[source]¶
zspsv
computes the solution to a complex system of linear equationswhere is an symmetric matrix stored in packed format and and are matrices.
For full information please refer to the NAG Library document for f07qn
https://support.nag.com/numeric/nl/nagdoc_30.2/flhtml/f07/f07qnf.html
- Parameters
- uplostr, length 1
If , the upper triangle of is stored.
If , the lower triangle of is stored.
- apcomplex, array-like, shape
The symmetric matrix , packed by columns.
- bcomplex, array-like, shape
The right-hand side matrix .
- Returns
- apcomplex, ndarray, shape
The block diagonal matrix and the multipliers used to obtain the factor or from the factorization or as computed by
zsptrf()
, stored as a packed triangular matrix in the same storage format as .- ipivint, ndarray, shape
Details of the interchanges and the block structure of . More precisely,
if , is a pivot block and the th row and column of were interchanged with the th row and column;
if and , is a pivot block and the th row and column of were interchanged with the th row and column;
if and , is a pivot block and the th row and column of were interchanged with the th row and column.
- bcomplex, ndarray, shape
If no exception or warning is raised, the solution matrix .
- Raises
- NagValueError
- (errno )
On entry, error in parameter .
Constraint: or .
- (errno )
On entry, error in parameter .
Constraint: .
- (errno )
On entry, error in parameter .
Constraint: .
- Warns
- NagAlgorithmicWarning
- (errno )
Element of the diagonal is exactly zero. The factorization has been completed, but the block diagonal matrix is exactly singular, so the solution could not be computed.
- Notes
zspsv
uses the diagonal pivoting method to factor as if or if , where (or ) is a product of permutation and unit upper (lower) triangular matrices, is symmetric and block diagonal with and diagonal blocks. The factored form of is then used to solve the system of equations .
- References
Anderson, E, Bai, Z, Bischof, C, Blackford, S, Demmel, J, Dongarra, J J, Du Croz, J J, Greenbaum, A, Hammarling, S, McKenney, A and Sorensen, D, 1999, LAPACK Users’ Guide, (3rd Edition), SIAM, Philadelphia, https://www.netlib.org/lapack/lug
Golub, G H and Van Loan, C F, 1996, Matrix Computations, (3rd Edition), Johns Hopkins University Press, Baltimore
Higham, N J, 2002, Accuracy and Stability of Numerical Algorithms, (2nd Edition), SIAM, Philadelphia