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

518HP MLIB User’s Guide