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

918 HP MLIB User’s Guide
DSLEIM Matrix structure input by matrix
INTEGER*8 neqns, nnzero, colstr(neqns+1), rowind(nnzero), ier
REAL*8 global(150)
CALL CSLEIM(colstr, rowind, global, ier)
INTEGER*8 neqns, nnzero, colstr(neqns+1), rowind(nnzero), ier
REAL*8 global(150)
CALL ZSLEIM(colstr, rowind, global, ier)
Input colstr For a symmetric matrix, colstr(j) gives the index in
rowind of the first nonzero in the lower triangular part
of column j of the matrix.
For a nonsymmetric matrix, colstr(j) gives the index in
rowind of the first nonzero in column j of the matrix.
The nonzeros for column j are found, in ascending
order, in rowind(colstr(j)), rowind(colstr(j)+1), ...,
rowind(colstr(j+1)−1).
colstr(neqns+1) must be set to one greater than the
total number of nonzeros, nnzero listed in rowind.
rowind For a symmetric matrix (LDL’ factorization), a
column-by-column list of row indices for all nonzeros in
the lower triangle of the matrix, in ascending order
within each column.
For a nonsymmetric matrix (LU factorizations), a
column-by-column list of row indices for all nonzeros in
the matrix, in ascending order within each column.
Input of diagonal entries is optional.
Updated global Global communications array for this problem. This
array must be passed, untouched by the user, to
successive subroutines in this package.
Output ier Status response:
ier = 0 Normal return.
ier = −100 Incorrect processing path; DSLEIN
not called or matrix input already
finished.
ier = −101 Error in dynamic storage allocation.
ier = −106 Illegal value for at least one entry in
colstr (not increasing or negative) or
invalid entries in rowind (entries out