naginterfaces.library.lapacklin.zpttrs

naginterfaces.library.lapacklin.zpttrs(uplo, d, e, b)

zpttrs computes the solution to a complex system of linear equations $AX=B$, where $A$ is an $n\times n$ Hermitian positive definite tridiagonal matrix and $X$ and $B$ are $n\times r$ matrices, using the $LD{L}^H$ factorization returned by zpttrf().

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

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

Parameters

uplostr, length 1

Specifies the form of the factorization as follows:

$uplo=\text{'U'}$

$A=U^{H}DU$.

$uplo=\text{'L'}$

$A=LD{L}^H$.

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

Must contain the $n$ diagonal elements of the diagonal matrix $D$ from the $LD{L}^H$ or $U^{H}DU$ factorization of $A$.

ecomplex, array-like, shape $(n-1)$

If $uplo=\text{'U'}$, $e$ must contain the $(n-1)$ superdiagonal elements of the unit upper bidiagonal matrix $U$ from the $U^{H}DU$ factorization of $A$.

If $uplo=\text{'L'}$, $e$ must contain the $(n-1)$ subdiagonal elements of the unit lower bidiagonal matrix $L$ from the $LD{L}^H$ factorization of $A$.

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

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

Returns

bcomplex, 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 $\text{n}$.

Constraint: $n\ge 0$.

(errno $-3$)

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

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

Notes

zpttrs should be preceded by a call to zpttrf(), which computes a modified Cholesky factorization of the matrix $A$ as

$$A=LD{L}^H\text{,}$$

where $L$ is a unit lower bidiagonal matrix and $D$ is a diagonal matrix, with positive diagonal elements. zpttrs then utilizes the factorization to solve the required equations. Note that the factorization may also be regarded as having the form $U^{H}DU$, where $U$ is a unit upper bidiagonal matrix.

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