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

 

 

xi, i = 1,..., n, is stored in an array X() with increment

 

 

argument incx. If incx > 0 then xiis stored in

 

 

X (1 + (i - 1) x incx). If incx < 0 then xi 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

 

 

yi, i = 1,..., n, is stored in an array Y() with increment

 

 

argument incy. If incy > 0 then yi is stored in

 

 

Y(1 + (i - 1) x incy). If incy < 0 then yi 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 wi, i = 1,..., n, is stored in an array W() with

 

 

increment argument incw. If incw > 0 then wi is stored

 

 

in W (1 + (i - 1) x incw). If incw < 0 then wi 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

 

 

zi, i = 1,..., n, is stored in an array Z() with increment

 

 

argument incw. If incw > 0 then zi is stored in

 

 

Z(1 + (i - 1) x incw). If incw < 0 then ziis 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

 

 

u1v1T and u2v2T , then the transpose of the updated

 

 

matrix is multiplied by vector y.

 

 

ˆ

T

T

 

 

A A + u1v1 + u

2v2

 

 

ˆ

T

 

 

 

x ← β A

y + z

 

370HP MLIB User’s Guide