HP MLIB User's Guide Vol. 1 7th Ed.

370 HP MLIB User’s Guide
F_SGEMVER/F_DGEMVER/F_CGEMVER/F_ZGEMVER Multiple matrix-vector multiply, rank 2 update
U1 REAL or COMPLEX array, dimension M.
V1 REAL or COMPLEX array, dimension N.
U2 REAL or COMPLEX array, dimension M.
V2 REAL or COMPLEX array, dimension N.
ALPHA The scalar ALPHA.
INCX Increment for the array x. A vector x having component
x
i
, i = 1,..., n, is stored in an array X() with increment
argument incx. If incx > 0 then x
i
is stored in
X (1 + (i - 1) x incx). If incx < 0 then x
i
is stored in
X (1 + (N - i) x |incx|). incx = 0 is an illegal value.
Y REAL or COMPLEX array, minimum length
(N - 1) x |incy| + 1.
INCY Increment for the array y. A vector y having component
y
i
, i = 1,..., n, is stored in an array Y() with increment
argument incy. If incy > 0 then y
i
is stored in
Y(1 + (i - 1) x incy). If incy < 0 then y
i
is stored in
Y(1 + (N - i) x |incy|). incy = 0 is an illegal value.
BETA The scalar BETA.
INCW Increment for the array w. A vector w having
component w
i
, i = 1,..., n, is stored in an array W() with
increment argument incw. If incw > 0 then w
i
is stored
in W (1 + (i - 1) x incw). If incw < 0 then w
i
is stored in
W (1 + (N - i) x |incw|). incw = 0 is an illegal value.
Z REAL or COMPLEX array, minimum length
(N - 1) x |incz| + 1.
INCZ Increment for the array y. A vector z having component
z
i
, i = 1,..., n, is stored in an array Z() with increment
argument incw. If incw > 0 then z
i
is stored in
Z(1 + (i - 1) x incw). If incw < 0 then z
i
is stored in
Z(1 + (N - i) x |incw|). incw = 0 is an illegal value.
Output X Contains the result after A is first updated by
and , then the transpose of the updated
matrix is multiplied by vector y.
u
1
v
1
T
u
2
v
2
T
A
ˆ
Au
1
v
1
T
u
2
v
2
T
++←
x βA
ˆ
T
yz+←