naginterfaces.library.lapacklin.zpotri

naginterfaces.library.lapacklin.zpotri(uplo, a)

zpotri computes the inverse of a complex Hermitian positive definite matrix $A$, where $A$ has been factorized by zpotrf().

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

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

Parameters

uplostr, length 1

Specifies how $A$ has been factorized.

$uplo=\text{'U'}$

$A=U^{H}U$, where $U$ is upper triangular.

$uplo=\text{'L'}$

$A=L{L}^H$, where $L$ is lower triangular.

acomplex, array-like, shape $(n,n)$

The upper triangular matrix $U$ if $uplo=\text{'U'}$ or the lower triangular matrix $L$ if $uplo=\text{'L'}$, as returned by zpotrf().

Returns

acomplex, ndarray, shape $(n,n)$

$U$ is overwritten by the upper triangle of $A^{-1}$ if $uplo=\text{'U'}$; $L$ is overwritten by the lower triangle of $A^{-1}$ if $uplo=\text{'L'}$.

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

Warns

NagAlgorithmicWarning

(errno $i>0$)

Diagonal element $⟨value⟩$ of the Cholesky factor is zero; the Cholesky factor is singular and the inverse of $A$ cannot be computed.

Notes

zpotri is used to compute the inverse of a complex Hermitian positive definite matrix $A$, the function must be preceded by a call to zpotrf(), which computes the Cholesky factorization of $A$.

If $uplo=\text{'U'}$, $A=U^{H}U$ and $A^{-1}$ is computed by first inverting $U$ and then forming $\left(U^{-1}\right)U^{-H}$.

If $uplo=\text{'L'}$, $A=L{L}^H$ and $A^{-1}$ is computed by first inverting $L$ and then forming $L^{-H}\left(L^{-1}\right)$.

References

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