NAG Library Manual, Mark 26
NAG AD Library Manual, Mark 26
NAG C Library Manual, Mark 26
F06 (blas) Chapter Contents
F06 (blas) Chapter Introduction
NAG Library Routine Document
f06srf (zher2)
Keyword Search:
NAG Library Manual, Mark 26
NAG AD Library Manual, Mark 26
NAG C Library Manual, Mark 26
F06 (blas) Chapter Contents
F06 (blas) Chapter Introduction
▸
▿
Contents
1
Purpose
2
Specification
3
Description
4
References
5
Arguments
6
Error Indicators and Warnings
7
Accuracy
8
Parallelism and Performance
9
Further Comments
10
Example
© The Numerical Algorithms Group Ltd. 2018
1
Purpose
f06srf (zher2)
computes the rank-2 update of a complex Hermitian matrix.
2
Specification
Fortran Interface
Subroutine f06srf (
uplo
,
n
,
alpha
,
x
,
incx
,
y
,
incy
,
a
,
lda
)
Integer, Intent (In)
::
n
,
incx
,
incy
,
lda
Complex (Kind=nag_wp), Intent (In)
::
alpha
,
x(*)
,
y(*)
Complex (Kind=nag_wp), Intent (Inout)
::
a(lda,*)
Character (1), Intent (In)
::
uplo
C Header Interface
#include <nagmk26.h>
void
f06srf_ (
const char *
uplo
,
const Integer *
n
,
const Complex *
alpha
,
const Complex
x
[]
,
const Integer *
incx
,
const Complex
y
[]
,
const Integer *
incy
,
Complex
a
[]
,
const Integer *
lda
,
const Charlen
length_uplo
)
The routine may be called by its BLAS name
zher2
.
3
Description
f06srf (zher2)
performs the Hermitian rank-2 update operation
A
←
α
x
y
H
+
α
-
y
x
H
+
A
,
where
A
is an
n
by
n
complex Hermitian matrix,
x
and
y
are
n
-element complex vectors, and
α
is a complex scalar.
4
References
None.
5
Arguments
1:
uplo
– Character(1)
Input
On entry
: specifies whether the upper or lower triangular part of
A
is stored.
uplo
=
'U'
The upper triangular part of
A
is stored.
uplo
=
'L'
The lower triangular part of
A
is stored.
Constraint
:
uplo
=
'U'
or
'L'
.
2:
n
– Integer
Input
On entry
:
n
, the order of the matrix
A
.
Constraint
:
n
≥
0
.
3:
alpha
– Complex (Kind=nag_wp)
Input
On entry
: the scalar
α
.
4:
x
*
– Complex (Kind=nag_wp) array
Input
Note:
the dimension of the array
x
must be at least
max
1
,
1
+
n
-
1
×
incx
.
On entry
: the
n
-element vector
x
.
If
incx
>
0
,
x
i
must be stored in
x
1
+
i
-
1
×
incx
, for
i
=
1
,
2
,
…
,
n
.
If
incx
<
0
,
x
i
must be stored in
x
1
-
n
-
i
×
incx
, for
i
=
1
,
2
,
…
,
n
.
Intermediate elements of
x
are not referenced.
5:
incx
– Integer
Input
On entry
: the increment in the subscripts of
x
between successive elements of
x
.
Constraint
:
incx
≠
0
.
6:
y
*
– Complex (Kind=nag_wp) array
Input
Note:
the dimension of the array
y
must be at least
max
1
,
1
+
n
-
1
×
incy
.
On entry
: the
n
-element vector
y
.
If
incy
>
0
,
y
i
must be stored in
y
1
+
i
-
1
×
incy
, for
i
=
1
,
2
,
…
,
n
.
If
incy
<
0
,
y
i
must be stored in
y
1
-
n
-
i
×
incy
, for
i
=
1
,
2
,
…
,
n
.
Intermediate elements of
y
are not referenced.
7:
incy
– Integer
Input
On entry
: the increment in the subscripts of
y
between successive elements of
y
.
Constraint
:
incy
≠
0
.
8:
a
lda
*
– Complex (Kind=nag_wp) array
Input/Output
Note:
the second dimension of the array
a
must be at least
max
1
,
n
.
On entry
: the
n
by
n
Hermitian matrix
A
.
If
uplo
=
'U'
, the upper triangular part of
A
must be stored and the elements of the array below the diagonal are not referenced.
If
uplo
=
'L'
, the lower triangular part of
A
must be stored and the elements of the array above the diagonal are not referenced.
On exit
: the updated matrix
A
. The imaginary parts of the diagonal elements are set to zero.
9:
lda
– Integer
Input
On entry
: the first dimension of the array
a
as declared in the (sub)program from which
f06srf (zher2)
is called.
Constraint
:
lda
≥
max
1
,
n
.
6
Error Indicators and Warnings
None.
7
Accuracy
Not applicable.
8
Parallelism and Performance
f06srf (zher2)
is not threaded in any implementation.
9
Further Comments
None.
10
Example
None.
NAG Library Manual, Mark 26
NAG AD Library Manual, Mark 26
NAG C Library Manual, Mark 26
F06 (blas) Chapter Contents
F06 (blas) Chapter Introduction
© The Numerical Algorithms Group Ltd. 2018