F06 Chapter Contents
F06 Chapter Introduction
NAG Library Manual

# NAG Library Routine DocumentF06YPF (DSYRK)

Note:  before using this routine, please read the Users' Note for your implementation to check the interpretation of bold italicised terms and other implementation-dependent details.

## 1  Purpose

F06YPF (DSYRK) performs one of the symmetric rank-$k$ update operations
 $C←αAAT + βC or C←αATA + βC ,$
where $A$ is a real matrix, $C$ is an $n$ by $n$ real symmetric matrix, and $\alpha$ and $\beta$ are real scalars.

## 2  Specification

 SUBROUTINE F06YPF ( UPLO, TRANS, N, K, ALPHA, A, LDA, BETA, C, LDC)
 INTEGER N, K, LDA, LDC REAL (KIND=nag_wp) ALPHA, A(LDA,*), BETA, C(LDC,*) CHARACTER(1) UPLO, TRANS
The routine may be called by its BLAS name dsyrk.

None.

None.

## 5  Parameters

1:     UPLO – CHARACTER(1)Input
On entry: specifies whether the upper or lower triangular part of $C$ is stored.
${\mathbf{UPLO}}=\text{'U'}$
The upper triangular part of $C$ is stored.
${\mathbf{UPLO}}=\text{'L'}$
The lower triangular part of $C$ is stored.
Constraint: ${\mathbf{UPLO}}=\text{'U'}$ or $\text{'L'}$.
2:     TRANS – CHARACTER(1)Input
On entry: specifies the operation to be performed.
${\mathbf{TRANS}}=\text{'N'}$
$C←\alpha A{A}^{\mathrm{T}}+\beta C$.
${\mathbf{TRANS}}=\text{'T'}$ or $\text{'C'}$
$C←\alpha {A}^{\mathrm{T}}A+\beta C$.
Constraint: ${\mathbf{TRANS}}=\text{'N'}$, $\text{'T'}$ or $\text{'C'}$.
3:     N – INTEGERInput
On entry: $n$, the order of the matrix $C$; the number of rows of $A$ if ${\mathbf{TRANS}}=\text{'N'}$, or the number of columns of $A$ if ${\mathbf{TRANS}}=\text{'T'}$ or $\text{'C'}$.
Constraint: ${\mathbf{N}}\ge 0$.
4:     K – INTEGERInput
On entry: $k$, the number of columns of $A$ if ${\mathbf{TRANS}}=\text{'N'}$, or the number of rows of $A$ if ${\mathbf{TRANS}}=\text{'T'}$ or $\text{'C'}$.
Constraint: ${\mathbf{K}}\ge 0$.
5:     ALPHA – REAL (KIND=nag_wp)Input
On entry: the scalar $\alpha$.
6:     A(LDA,$*$) – REAL (KIND=nag_wp) arrayInput
Note: the second dimension of the array A must be at least $\mathrm{max}\phantom{\rule{0.125em}{0ex}}\left(1,{\mathbf{K}}\right)$ if ${\mathbf{TRANS}}=\text{'N'}$ and at least $\mathrm{max}\phantom{\rule{0.125em}{0ex}}\left(1,{\mathbf{N}}\right)$ if ${\mathbf{TRANS}}=\text{'T'}$ or $\text{'C'}$.
On entry: the matrix $A$; $A$ is $n$ by $k$ if ${\mathbf{TRANS}}=\text{'N'}$, or $k$ by $n$ if ${\mathbf{TRANS}}=\text{'T'}$ or $\text{'C'}$.
7:     LDA – INTEGERInput
On entry: the first dimension of the array A as declared in the (sub)program from which F06YPF (DSYRK) is called.
Constraints:
• if ${\mathbf{TRANS}}=\text{'N'}$, ${\mathbf{LDA}}\ge \mathrm{max}\phantom{\rule{0.125em}{0ex}}\left(1,{\mathbf{N}}\right)$;
• if ${\mathbf{TRANS}}=\text{'T'}$ or $\text{'C'}$, ${\mathbf{LDA}}\ge \mathrm{max}\phantom{\rule{0.125em}{0ex}}\left(1,{\mathbf{K}}\right)$.
8:     BETA – REAL (KIND=nag_wp)Input
On entry: the scalar $\beta$.
9:     C(LDC,$*$) – REAL (KIND=nag_wp) arrayInput/Output
Note: the second dimension of the array C must be at least $\mathrm{max}\phantom{\rule{0.125em}{0ex}}\left(1,{\mathbf{N}}\right)$.
On entry: the $n$ by $n$ symmetric matrix $C$.
• If ${\mathbf{UPLO}}=\text{'U'}$, the upper triangular part of $C$ must be stored and the elements of the array below the diagonal are not referenced.
• If ${\mathbf{UPLO}}=\text{'L'}$, the lower triangular part of $C$ must be stored and the elements of the array above the diagonal are not referenced.
On exit: the updated matrix $C$.
10:   LDC – INTEGERInput
On entry: the first dimension of the array C as declared in the (sub)program from which F06YPF (DSYRK) is called.
Constraint: ${\mathbf{LDC}}\ge \mathrm{max}\phantom{\rule{0.125em}{0ex}}\left(1,{\mathbf{N}}\right)$.

None.

Not applicable.