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

502 HP MLIB User’s Guide
SDIASM/DDIASM/CDIASM/ZDIASM Diagonal format triangular solve
SUBROUTINE ZDIASM
INTEGER*4 transa, m, n, unitd, lda, ndiag, ldb, ldc, lwork
INTEGER*4 descra(*), idiag(*)
COMPLEX*16 alpha, beta
COMPLEX*16 val(*), b(ldb,*), c(ldc,*), work(*)
CALL ZDIASM (transa, m, n, unitd, dv, alpha, descra, val, lda, idiag,
ndiag, b, ldb, beta, c, ldc, work, lwork)
VECLIB8:
SUBROUTINE SDIASM
INTEGER*8 transa, m, n, unitd, lda, ndiag, ldb, ldc, lwork
INTEGER*8 descra(*), idiag(*)
REAL*4 alpha, beta
REAL*4 val(*), b(ldb,*), c(ldc,*), work(*)
CALL SDIASM (transa, m, n, unitd, dv, alpha, descra, val, lda, idiag,
ndiag, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE DDIASM
INTEGER*8 transa, m, n, unitd, lda, ndiag, ldb, ldc, lwork
INTEGER*8 descra(*), idiag(*)
REAL*8 alpha, beta
REAL*8 val(*), b(ldb,*), c(ldc,*), work(*)
CALL DDIASM (transa, m, n, unitd, dv, alpha, descra, val, lda, idiag,
ndiag, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE CDIASM
INTEGER*8 transa, m, n, unitd, lda, ndiag, ldb, ldc, lwork
INTEGER*8 descra(*), idiag(*)
COMPLEX*8 alpha, beta
COMPLEX*8 val(*), b(ldb,*), c(ldc,*), work(*)
CALL CDIASM (transa, m, n, unitd, dv, alpha, descra, val, lda, idiag,
ndiag, b, ldb, beta, c, ldc, work, lwork)
SUBROUTINE ZDIASM
INTEGER*8 transa, m, n, unitd, lda, ndiag, ldb, ldc, lwork
INTEGER*8 descra(*), idiag(*)
COMPLEX*16 alpha, beta
COMPLEX*16 val(*), b(ldb,*), c(ldc,*), work(*)
CALL ZDIASM (transa, m, n, unitd, dv, alpha, descra, val, lda, idiag,
ndiag, b, ldb, beta, c, ldc, work, lwork)
Input transa Indicates how to operate with the sparse matrix.
0: Operate with matrix