naginterfaces.library.lapacklin.dsptri

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

dsptri computes the inverse of a real symmetric indefinite matrix , where has been factorized by dsptrf(), using packed storage.

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

https://www.nag.com/numeric/nl/nagdoc_27.3/flhtml/f07/f07pjf.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().

Returns
apfloat, ndarray, shape

The factorization is overwritten by the matrix .

Raises
NagValueError
(errno )

On entry, error in parameter .

Constraint: or .

(errno )

On entry, error in parameter .

Constraint: .

Warns
NagAlgorithmicWarning
(errno )

Element of the diagonal is exactly zero. is singular and the inverse of cannot be computed.

Notes

dsptri is used to compute the inverse of a real symmetric indefinite matrix , the function must be preceded by a call to dsptrf(), which computes the Bunch–Kaufman factorization of , using packed storage.

If , and is computed by solving .

If , and is computed by solving .

References

Du Croz, J J and Higham, N J, 1992, Stability of methods for matrix inversion, IMA J. Numer. Anal. (12), 1–19