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

Chapter 3 Basic Matrix Operations 373
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 , the
subprograms do not reference A, X, or Y.
N Number of columns in matrix A, n > 0. If , 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
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 REAL or COMPLEX 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.
m 0≤
n 0≤