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

372 HP MLIB User’s Guide
F_SGEMVT/F_DGEMVT/F_CGEMVT/F_ZGEMVT Multiple matrix-vector multiply
Name F_SGEMVT/F_DGEMVT/F_CGEMVT/F_ZGEMVT
Multiple matrix-vector multiply
Purpose F_xGEMVT combines a matrix-vector and a transposed matrix-vector multiply.
Specifically, F_xGEMVT routines perform the following operations:
1. Multiply a vector y by a general matrix A
T
.
2. Scale the resulting vector by β and store the result in the vector operand x.
3. Multiply the matrix by the resultant vector x.
4. Scale the resulting vector by α and store it in the vector operand w.
If m or n is less than or equal to zero, this function returns immediately.
Usage VECLIB
INTEGER*4 INCX, INCY, LDA, M, N
REAL*4 ALPHA, BETA, A( LDA, * ), X( * ), Y( * ), W( * ), Z( * )
SUBROUTINE F_SGEMVT (M, N, ALPHA, A, LDA, X, INCX, Y, INCY,
BETA, W, INCW, Z, INCZ)
INTEGER*4 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*4 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*4 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)
VECLIB8
INTEGER*8 INCX, INCY, LDA, M, N
REAL*4 ALPHA, BETA, A( LDA, * ), X( * ), Y( * ), W( * ), Z( * )
SUBROUTINE F_SGEMVT (M, N, ALPHA, A, LDA, X, INCX, Y, INCY,
BETA, W, INCW, Z, INCZ)
x βA
T
yz+←
w αAx←