naginterfaces.library.lapacklin.dsptrs

naginterfaces.library.lapacklin.dsptrs(uplo, n, ap, ipiv, b)[source]

dsptrs solves a real symmetric indefinite system of linear equations with multiple right-hand sides,

where has been factorized by dsptrf(), using packed storage.

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

https://www.nag.com/numeric/nl/nagdoc_27.3/flhtml/f07/f07pef.html

Parameters
uplostr, length 1

Specifies how has been factorized.

, where is upper triangular.

, where is lower triangular.

nint

, the order of the matrix .

apfloat, array-like, shape

The factorization of stored in packed form, as returned by dsptrf().

ipivint, array-like, shape

Details of the interchanges and the block structure of , as returned by dsptrf().

bfloat, array-like, shape

The right-hand side matrix .

Returns
bfloat, ndarray, shape

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: .

Notes

dsptrs is used to solve a real symmetric indefinite system of linear equations , the function must be preceded by a call to dsptrf() which computes the Bunch–Kaufman factorization of , using packed storage.

If , , where is a permutation matrix, is an upper triangular matrix and is a symmetric block diagonal matrix with and blocks; the solution is computed by solving and then .

If , , where is a lower triangular matrix; the solution is computed by solving and then .

References

Golub, G H and Van Loan, C F, 1996, Matrix Computations, (3rd Edition), Johns Hopkins University Press, Baltimore