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

518 HP MLIB User’s Guide
SJADSM/DJADSM/CJADSM/ZJADSM Jagged diagonal format triangular solve
SUBROUTINE ZJADSM
INTEGER*4 transa, m, n, unitd, maxnz, ldb, ldc, lwork
INTEGER*4 descra(*), indx(*), pntr(*), iperm(*)
COMPLEX*16 alpha, beta
COMPLEX*16 val(*), b(ldb,*), c(ldc,*), work(*)
CALL ZJADSM (transa, m, n, unitd, dv, alpha, descra, val, pntr, iperm,
indx, maxnz, b, ldb, beta, c, ldc, work, lwork)
VECLIB8:
SUBROUTINE SJADSM
INTEGER*8 transa, m, n, unitd, maxnz, ldb, ldc, lwork
INTEGER*8 descra(*), indx(*), pntr(*), iperm(*)
REAL*4 alpha, beta
REAL*4 val(*), b(ldb,*), c(ldc,*), work(*)
CALL SJADSM (transa, m, n, unitd, dv, alpha, descra, val, pntr, iperm,
indx, maxnz, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE DJADSM
INTEGER*8 transa, m, n, unitd, maxnz, ldb, ldc, lwork
INTEGER*8 descra(*), indx(*), pntr(*), iperm(*)
REAL*8 alpha, beta
REAL*8 val(*), b(ldb,*), c(ldc,*), work(*)
CALL DJADSM (transa, m, n, unitd, dv, alpha, descra, val, pntr, iperm,
indx, maxnz, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE CJADSM
INTEGER*8 transa, m, n, unitd, maxnz, ldb, ldc, lwork
INTEGER*8 descra(*), indx(*), pntr(*), iperm(*)
COMPLEX*8 alpha, beta
COMPLEX*8 val(*), b(ldb,*), c(ldc,*), work(*)
CALL CJADSM (transa, m, n, unitd, dv, alpha, descra, val, pntr, iperm,
indx, maxnz, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE ZJADSM
INTEGER*8 transa, m, n, unitd, maxnz, ldb, ldc, lwork
INTEGER*8 descra(*), indx(*), pntr(*), iperm(*)
COMPLEX*16 alpha, beta
COMPLEX*16 val(*), b(ldb,*), c(ldc,*), work(*)
CALL ZJADSM (transa, m, n, unitd, dv, alpha, descra, val, pntr, iperm,
indx, maxnz, b, ldb, beta, c, ldc, work, lwork)
Input transa Indicates how to operate with the sparse matrix.
0: Operate with matrix