Multiple matrix-vector multiply

F_SGEMVT/F_DGEMVT/F_CGEMVT/F_ZGEMVT

 

INTEGER*8

INCX, INCY, LDA, M, N

 

REAL*8

ALPHA, BETA, A( LDA, * ), X( * ), Y( * ), W( * ), Z( * )

 

SUBROUTINE F_DGEMVT (M, N, ALPHA, A, LDA, X, INCX, Y, INCY,

 

BETA, W, INCW, Z, INCZ)

 

INTEGER*8

INCX, INCY, LDA, M, N

 

COMPLEX*8

ALPHA, BETA, A( LDA, * ), X( * ), Y( * ), W( * ), Z( * )

 

SUBROUTINE F_CGEMVT (M, N, ALPHA, A, LDA, X, INCX, Y, INCY,

 

BETA, W, INCW, Z, INCZ)

 

INTEGER*8

INCX, INCY, LDA, M, N

 

COMPLEX*16

ALPHA, BETA, A( LDA, * ), X( * ), Y( * ), W( * ), Z( * )

 

SUBROUTINE F_ZGEMVT (M, N, ALPHA, A, LDA, X, INCX, Y, INCY,

 

BETA, W, INCW, Z, INCZ)

Input

M

Number of rows in matrix A, m > 0. If m 0 , the

 

 

subprograms do not reference A, X, or Y.

 

N

Number of columns in matrix A, n > 0. If n 0 , the

 

 

subprograms do not reference A, X, or Y.

 

ALPHA

REAL or COMPLEX scalar ALPHA.

 

A

REAL or COMPLEX array, dimension (LDA, N).

 

LDA

Leading dimension of array A. If lda < 1 or lda < m, an

 

 

error flag is set and passed to the error handler.

 

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

REAL or COMPLEX 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.

Chapter 3 Basic Matrix Operations 373