naginterfaces.library.blas.dtrsm

naginterfaces.library.blas.dtrsm(side, uplo, transa, diag, m, alpha, a, b)

dtrsm performs one of the matrix-matrix operations

$$\begin{array}{c}\begin{array}{cc}B\leftarrow \alpha A^{-1}B\text{,}& B\leftarrow \alpha A^{-T}B\text{,}\\ B\leftarrow \alpha B{A}^{-1}\text{or}& B\leftarrow \alpha B{A}^{-T}\text{,}\end{array}\end{array}$$

where $A$ is a real triangular matrix, $B$ is an $m\times n$ real matrix, and $\alpha$ is a real scalar. $A^{-T}$ denotes ${\left(A^T\right)}^{-1}$ or equivalently ${\left(A^{-1}\right)}^T$.

No test for singularity or near-singularity of $A$ is included in this function. Such tests must be performed before calling this function.

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

https://support.nag.com/numeric/nl/nagdoc_30.3/flhtml/f06/f06yjf.html

Parameters

sidestr, length 1

Specifies whether $B$ is operated on from the left or the right.

$side=\text{'L'}$

$B$ is pre-multiplied from the left.

$side=\text{'R'}$

$B$ is post-multiplied from the right.

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.

transastr, length 1

Specifies whether the operation involves $A^{-1}$ or $A^{-T}$.

$transa=\text{'N'}$

The operation involves $A^{-1}$.

$transa=\text{'T'}$ or $\text{'C'}$

The operation involves $A^{-T}$.

diagstr, length 1

Specifies whether $A$ has nonunit or unit diagonal elements.

$diag=\text{'N'}$

The diagonal elements are stored explicitly.

$diag=\text{'U'}$

The diagonal elements are assumed to be $1$, and are not referenced.

mint

$m$, the number of rows of the matrix $B$; the order of $A$ if $side=\text{'L'}$.

alphafloat

The scalar $\alpha$.

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

Note: the required extent for this argument in dimension 1 is determined as follows: if $side=\text{'L'}$: $m$; if $side=\text{'R'}$: $n$; otherwise: $0$.

Note: the required extent for this argument in dimension 2 is determined as follows: if $side=\text{'L'}$: $m$; if $side=\text{'R'}$: $n$; otherwise: $0$.

The triangular matrix $A$; $A$ is $m\times m$ if $side=\text{'L'}$, or $n\times n$ if $side=\text{'R'}$.

bfloat, array-like, shape $(m,n)$

The $m\times n$ matrix $B$.

If $alpha=0$, $b$ need not be set.

Returns

bfloat, ndarray, shape $(m,n)$

The updated matrix $B$.

Raises

NagValueError

(errno $1$)

On entry, error in parameter $side$.

Constraint: $side=\text{'L'}$ or $\text{'R'}$.

(errno $2$)

On entry, error in parameter $uplo$.

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

(errno $3$)

On entry, error in parameter $transa$.

Constraint: $transa=\text{'N'}$, $\text{'T'}$ or $\text{'C'}$.

(errno $4$)

On entry, error in parameter $diag$.

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

(errno $5$)

On entry, error in parameter $m$.

Constraint: $m\ge 0$.

(errno $6$)

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

Constraint: $n\ge 0$.

Notes

In the NAG Library the traditional C interface for this routine uses a different algorithmic base. Please contact NAG if you have any questions about compatibility.