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

462 HP MLIB User’s Guide
SBSCMM/DBSCMM/CBSCMM/ZBSCMM Block sparse column matrix-matrix multiply
SUBROUTINE ZBSCMM
INTEGER*4 transa, mb, n, kb, lb, ldb, ldc, lwork
INTEGER*4 descra(*), bindx(*), bpntrb(*), bpntre(*)
COMPLEX*16 alpha, beta
COMPLEX*16 val(*), b(ldb,*), c(ldc,*), work(*)
CALL ZBSCMM (transa, mb, n, kb, alpha, descra, val, bindx, bpntrb,
bpntre, lb, b, ldb, beta, c, ldc, work, lwork)
VECLIB8:
SUBROUTINE SBSCMM
INTEGER*8 transa, mb, n, kb, lb, ldb, ldc, lwork
INTEGER*8 descra(*), bindx(*), bpntrb(*), bpntre(*)
REAL*4 alpha, beta
REAL*4 val(*), b(ldb,*), c(ldc,*), work(*)
CALL SBSCMM (transa, mb, n, kb, alpha, descra, val, bindx, bpntrb,
bpntre, lb, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE DBSCMM
INTEGER*8 transa, mb, n, kb, lb, ldb, ldc, lwork
INTEGER*8 descra(*), bindx(*), bpntrb(*), bpntre(*)
REAL*8 alpha, beta
REAL*8 val(*), b(ldb,*), c(ldc,*), work(*)
CALL DBSCMM (transa, mb, n, kb, alpha, descra, val, bindx, bpntrb,
bpntre, lb, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE CBSCMM
INTEGER*8 transa, mb, n, kb, lb, ldb, ldc, lwork
INTEGER*8 descra(*), bindx(*), bpntrb(*), bpntre(*)
COMPLEX*8 alpha, beta
COMPLEX*8 val(*), b(ldb,*), c(ldc,*), work(*)
CALL CBSCMM (transa, mb, n, kb, alpha, descra, val, bindx, bpntrb,
bpntre, lb, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE ZBSCMM
INTEGER*8 transa, mb, n, kb, lb, ldb, ldc, lwork
INTEGER*8 descra(*), bindx(*), bpntrb(*), bpntre(*)
COMPLEX*16 alpha, beta
COMPLEX*16 val(*), b(ldb,*), c(ldc,*), work(*)
CALL ZBSCMM (transa, mb, n, kb, alpha, descra, val, bindx, bpntrb,
bpntre, lb, b, ldb, beta, c, ldc, work, lwork)
Input transa Indicates how to operate with the sparse matrix.
0: Operate with matrix