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

458 HP MLIB User’s Guide
SBELSM/DBELSM/CBELSM/ZBELSM Block Ellpack format triangular solve
SUBROUTINE ZBELSM
INTEGER*4 transa, mb, n, unitd, blda, maxbnz, lb, ldb, ldc, lwork
INTEGER*4 descra(*), bindx(*)
COMPLEX*16 alpha, beta
COMPLEX*16 val(*), b(ldb,*), c(ldc,*), work(*)
CALL ZBELSM (transa, mb, n, unitd, dv, alpha, descra, val, blda, bindx,
maxbnz, lb, b, ldb, beta, c, ldc, work, lwork)
VECLIB8:
SUBROUTINE SBELSM
INTEGER*8 transa, mb, n, unitd, blda, maxbnz, lb, ldb, ldc, lwork
INTEGER*8 descra(*), bindx(*)
REAL*4 alpha, beta
REAL*4 val(*), b(ldb,*), c(ldc,*), work(*)
CALL SBELSM (transa, mb, n, unitd, dv, alpha, descra, val, blda, bindx,
maxbnz, lb, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE DBELSM
INTEGER*8 transa, mb, n, unitd, blda, maxbnz, lb, ldb, ldc, lwork
INTEGER*8 descra(*), bindx(*)
REAL*8 alpha, beta
REAL*8 val(*), b(ldb,*), c(ldc,*), work(*)
CALL DBELSM (transa, mb, n, unitd, dv, alpha, descra, val, blda, bindx,
maxbnz, lb, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE CBELSM
INTEGER*8 transa, mb, n, unitd, blda, maxbnz, lb, ldb, ldc, lwork
INTEGER*8 descra(*), bindx(*)
COMPLEX*8 alpha, beta
COMPLEX*8 val(*), b(ldb,*), c(ldc,*), work(*)
CALL CBELSM (transa, mb, n, unitd, dv, alpha, descra, val, blda, bindx,
maxbnz, lb, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE ZBELSM
INTEGER*8 transa, mb, n, unitd, blda, maxbnz, lb, ldb, ldc, lwork
INTEGER*8 descra(*), bindx(*)
COMPLEX*16 alpha, beta
COMPLEX*16 val(*), b(ldb,*), c(ldc,*), work(*)
CALL ZBELSM (transa, mb, n, unitd, dv, alpha, descra, val, blda, bindx,
maxbnz, lb, b, ldb, beta, c, ldc, work, lwork)
Input transa Indicates how to operate with the sparse matrix.
0: Operate with matrix