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)

Chapter 3 Basic Matrix Operations 223