naginterfaces.library.lapacklin.dsptrs

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

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

$$AX=B\text{,}$$

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

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

https://support.nag.com/numeric/nl/nagdoc_30.3/flhtml/f07/f07pef.html

Parameters

uplostr, length 1

Specifies how $A$ has been factorized.

$uplo=\text{'U'}$

$A=PUD{U}^T{P}^T$, where $U$ is upper triangular.

$uplo=\text{'L'}$

$A=PLD{L}^T{P}^T$, where $L$ is lower triangular.

nint

$n$, the order of the matrix $A$.

apfloat, array-like, shape $(n\times (n+1)/2)$

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

ipivint, array-like, shape $\left(n\right)$

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

bfloat, array-like, shape $(n,\text{nrhs})$

The $n\times r$ right-hand side matrix $B$.

Returns

bfloat, ndarray, shape $(n,\text{nrhs})$

The $n\times r$ solution matrix $X$.

Raises

NagValueError

(errno $-1$)

On entry, error in parameter $uplo$.

Constraint: $uplo=\text{'U'}$ or $\text{'L'}$.

(errno $-2$)

On entry, error in parameter $n$.

Constraint: $n\ge 0$.

(errno $-3$)

On entry, error in parameter $\text{nrhs}$.

Constraint: $\text{nrhs}\ge 0$.

Notes

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

If $uplo=\text{'U'}$, $A=PUD{U}^T{P}^T$, where $P$ is a permutation matrix, $U$ is an upper triangular matrix and $D$ is a symmetric block diagonal matrix with $1\times 1$ and $2\times 2$ blocks; the solution $X$ is computed by solving $PUDY=B$ and then $U^T{P}^{T}X=Y$.

If $uplo=\text{'L'}$, $A=PLD{L}^T{P}^T$, where $L$ is a lower triangular matrix; the solution $X$ is computed by solving $PLDY=B$ and then $L^T{P}^{T}X=Y$.

References

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