Developer Reference

  • 2020.2
  • 07/15/2020
  • Public Content
Contents

p?her2

Performs a rank-2 update of a distributed Hermitian matrix.

Syntax

call pcher2
(
uplo
,
n
,
alpha
,
x
,
ix
,
jx
,
descx
,
incx
,
y
,
iy
,
jy
,
descy
,
incy
,
a
,
ia
,
ja
,
desca
)
call pzher2
(
uplo
,
n
,
alpha
,
x
,
ix
,
jx
,
descx
,
incx
,
y
,
iy
,
jy
,
descy
,
incy
,
a
,
ia
,
ja
,
desca
)
Include Files
  • mkl_pblas.h
Description
The
p?her2
routines perform a distributed matrix-vector operation defined as
sub(
A
) :=
alpha
*sub(
x
)*conj(sub(
y
)')+ conj(
alpha
)*sub(
y
)*conj(sub(
x
)') + sub(
A
),
where:
alpha
is a scalar,
sub(
A
)
is a
n
-by-
n
distributed Hermitian matrix,
sub(
A
)=
A
(
ia
:
ia
+
n
-1,
ja
:
ja
+
n
-1)
,
sub(
x
)
and
sub(
y
)
are distributed vectors.
sub(
x
)
denotes
X
(
ix
,
jx
:
jx
+
n
-1)
if
incx
=
m_x
, and
X
(
ix
:
ix
+
n
-1,
jx
)
if
incx
= 1,
sub(
y
)
denotes
Y
(
iy
,
jy
:
jy
+
n
-1)
if
incy
=
m_y
, and
Y
(
iy
:
iy
+
n
-1,
jy
)
if
incy
= 1
.
Input Parameters
uplo
(global)
CHARACTER*1
.
Specifies whether the upper or lower triangular part of the distributed Hermitian matrix
sub(
A
)
is used:
If
uplo
=
'U'
or
'u'
, then the upper triangular part of the
sub(
A
)
is used.
If
uplo
=
'L'
or
'l'
, then the low triangular part of the
sub(
A
)
is used.
n
(global)
INTEGER
.
Specifies the order of the distributed matrix
sub(
A
)
,
n
0.
alpha
(global)
COMPLEX
for
pcher2
DOUBLE COMPLEX
for
pzher2
Specifies the scalar
alpha
.
x
(local)
COMPLEX
for
pcher2
DOUBLE COMPLEX
for
pzher2
Array, size at least
(
jx
-1)*
m_x
+
ix
+(
n
-1)*abs(
incx
))
.
This array contains the entries of the distributed vector
sub(
x
)
.
ix
,
jx
(global)
INTEGER
.
The row and column indices in the distributed matrix
X
indicating the first row and the first column of the submatrix
sub(
x
)
, respectively.
descx
(global and local)
INTEGER
array of dimension 9. The array descriptor of the distributed matrix
X
.
incx
(global)
INTEGER
.
Specifies the increment for the elements of
sub(
x
)
. Only two values are supported, namely 1 and
m_x
.
incx
must not be zero.
y