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

Chapter 3 Basic Matrix Operations 223
Matrix-matrix multiply SGEMM/DGEMM/CGEMM/ZGEMM
VECLIB8:
CHARACTER*1 transa, transb
INTEGER*8 m, n, k, lda, ldb, ldc
REAL*4 alpha, beta, a(lda, *), b(ldb, *), c(ldc, n)
CALL SGEMM(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
CHARACTER*1 transa, transb
INTEGER*8 m, n, k, lda, ldb, ldc
REAL*8 alpha, beta, a(lda, *), b(ldb, *), c(ldc, n)
CALL DGEMM(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
CHARACTER*1 transa, transb
INTEGER*8 m, n, k, lda, ldb, ldc
COMPLEX*8 alpha, beta, a(lda, *), b(ldb, *), c(ldc, n)
CALL CGEMM(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
CHARACTER*1 transa, transb
INTEGER*8 m, n, k, lda, ldb, ldc
COMPLEX*16 alpha, beta, a(lda, *), b(ldb, *), c(ldc, n)
CALL ZGEMM(transa, transb, m, n, k, alpha, a, lda, b, ldb, beta, c, ldc)