naginterfaces.library.lapacklin.dspsv

naginterfaces.library.lapacklin.dspsv(uplo, ap, b)[source]

dspsv computes the solution to a real system of linear equations

where is an symmetric matrix stored in packed format and and are matrices.

For full information please refer to the NAG Library document for f07pa

https://support.nag.com/numeric/nl/nagdoc_30.2/flhtml/f07/f07paf.html

Parameters
uplostr, length 1

If , the upper triangle of is stored.

If , the lower triangle of is stored.

apfloat, array-like, shape

The symmetric matrix , packed by columns.

bfloat, array-like, shape

The right-hand side matrix .

Returns
apfloat, ndarray, shape

The block diagonal matrix and the multipliers used to obtain the factor or from the factorization or as computed by dsptrf(), 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.

bfloat, 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

dspsv 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