naginterfaces.library.lapacklin.dtrtri

naginterfaces.library.lapacklin.dtrtri(uplo, diag, a)

dtrtri computes the inverse of a real triangular matrix.

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

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

Parameters

uplostr, length 1

Specifies whether $A$ is upper or lower triangular.

$uplo=\text{'U'}$

$A$ is upper triangular.

$uplo=\text{'L'}$

$A$ is lower triangular.

diagstr, length 1

Indicates whether $A$ is a nonunit or unit triangular matrix.

$diag=\text{'N'}$

$A$ is a nonunit triangular matrix.

$diag=\text{'U'}$

$A$ is a unit triangular matrix; the diagonal elements are not referenced and are assumed to be $1$.

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

The $n\times n$ triangular matrix $A$.

Returns

afloat, ndarray, shape $(n,n)$

$A$ is overwritten by $A^{-1}$, using the same storage format as described above.

Raises

NagValueError

(errno $-1$)

On entry, error in parameter $uplo$.

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

(errno $-2$)

On entry, error in parameter $diag$.

Constraint: $diag=\text{'N'}$ or $\text{'U'}$.

(errno $-3$)

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

Constraint: $n\ge 0$.

Warns

NagAlgorithmicWarning

(errno $i>0$)

Element $⟨value⟩$ of the diagonal is exactly zero. $A$ is singular its inverse cannot be computed.

Notes

dtrtri forms the inverse of a real triangular matrix $A$. Note that the inverse of an upper (lower) triangular matrix is also upper (lower) triangular.

References

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