naginterfaces.library.blas.dsyrk

naginterfaces.library.blas.dsyrk(uplo, trans, k, alpha, a, beta, c)

dsyrk performs one of the symmetric rank-$k$ update operations

$$C\leftarrow \alpha A{A}^T+\beta C\text{or}C\leftarrow \alpha A^{T}A+\beta C\text{,}$$

where $A$ is a real matrix, $C$ is an $n\times n$ real symmetric matrix, and $\alpha$ and $\beta$ are real scalars.

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

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

Parameters

uplostr, length 1

Specifies whether the upper or lower triangular part of $C$ is stored.

$uplo=\text{'U'}$

The upper triangular part of $C$ is stored.

$uplo=\text{'L'}$

The lower triangular part of $C$ is stored.

transstr, length 1

Specifies the operation to be performed.

$trans=\text{'N'}$

$C\leftarrow \alpha A{A}^T+\beta C$.

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

$C\leftarrow \alpha A^{T}A+\beta C$.

kint

$k$, the number of columns of $A$ if $trans=\text{'N'}$, or the number of rows of $A$ if $trans=\text{'T'}$ or $\text{'C'}$.

alphafloat

The scalar $\alpha$.

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

Note: the required extent for this argument in dimension 1 is determined as follows: if $trans=\text{'N'}$: $n$; if $trans\text{in}(\text{'T'},\text{'C'})$: $k$; otherwise: $0$.

Note: the required extent for this argument in dimension 2 is determined as follows: if $trans=\text{'N'}$: $k$; if $trans\text{in}(\text{'T'},\text{'C'})$: $n$; otherwise: $0$.

The matrix $A$; $A$ is $n\times k$ if $trans=\text{'N'}$, or $k\times n$ if $trans=\text{'T'}$ or $\text{'C'}$.

betafloat

The scalar $\beta$.

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

The $n\times n$ symmetric matrix $C$.

Returns

cfloat, ndarray, shape $(n,n)$

The updated matrix $C$.

Raises

NagValueError

(errno $1$)

On entry, error in parameter $uplo$.

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

(errno $2$)

On entry, error in parameter $trans$.

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

(errno $3$)

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

Constraint: $n\ge 0$.

(errno $4$)

On entry, error in parameter $k$.

Constraint: $k\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.