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

Chapter 3 Basic Matrix Operations 285
Rank-2k update SSYR2K/DSYR2K/CHER2K/CSYR2K/ZHER2K/ZSYR2K
CHARACTER*1 uplo, trans
INTEGER*4 n, k, lda, ldb, ldc
REAL*4 beta
COMPLEX*8 alpha, a(lda, *), b(ldb, *), c(ldc, n)
CALL CHER2K(uplo, trans, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
CHARACTER*1 uplo, trans
INTEGER*4 n, k, lda, ldb, ldc
COMPLEX*8 alpha, beta, a(lda, *), b(ldb, *), c(ldc, n)
CALL CSYR2K(uplo, trans, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
CHARACTER*1 uplo, trans
INTEGER*4 n, k, lda, ldb, ldc
REAL*8 beta
COMPLEX*16 alpha, a(lda, *), b(ldb, *), c(ldc, n)
CALL ZHER2K(uplo, trans, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
CHARACTER*1 uplo, trans
INTEGER*4 n, k, lda, ldb, ldc
COMPLEX*16 alpha, beta, a(lda, *), b(ldb, *), c(ldc, n)
CALL ZSYR2K(uplo, trans, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
VECLIB8:
CHARACTER*1 uplo, trans
INTEGER*8 n, k, lda, ldb, ldc
REAL*4 alpha, beta, a(lda, *), b(ldb, *), c(ldc, n)
CALL SSYR2K(uplo, trans, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
CHARACTER*1 uplo, trans
INTEGER*8 n, k, lda, ldb, ldc
REAL*8 alpha, beta, a(lda, *), b(ldb, *), c(ldc, n)
CALL DSYR2K(uplo, trans, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
CHARACTER*1 uplo, trans
INTEGER*8 n, k, lda, ldb, ldc
REAL*4 beta
COMPLEX*8 alpha, a(lda, *), b(ldb, *), c(ldc, n)
CALL CHER2K(uplo, trans, n, k, alpha, a, lda, b, ldb, beta, c, ldc)
CHARACTER*1 uplo, trans
INTEGER*8 n, k, lda, ldb, ldc
COMPLEX*8 alpha, beta, a(lda, *), b(ldb, *), c(ldc, n)
CALL CSYR2K(uplo, trans, n, k, alpha, a, lda, b, ldb, beta, c, ldc)