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

Chapter 13 Sparse Linear Equations 911
Matrix structure input by column DSLEIC
Name DSLEIC
Matrix structure input by column
Purpose This subprogram adds a list of indices in a single column to the set of known
nonzeros for the sparse matrix.
NOTE Matrix structure input by matrix supports five data types. Matrix structure
input by elements, columns, and finite elements supports only REAL*8
precision.
Usage VECLIB:
INTEGER*4 jcol, nzcol, jrowin(nzcol), ier
REAL*8 global(150)
CALL DSLEIC(jcol, nzcol, jrowin, global, ier)
VECLIB8:
INTEGER*8 jcol, nzcol, jrowin(nzcol), ier
REAL*8 global(150)
CALL DSLEIC(jcol, nzcol, jrowin, global, ier)
Input jcol Column index; 1 ≤ jcol ≤ neqns. Columns must be
presented in order from 1 to neqns, except that
columns with no entries can be skipped.
nzcol Number of nonzeros in column jcol of the matrix;
nzcol ≥ 0.
jrowin For a symmetric matrix, (LDL’ factorization), a list of
row indices for all nonzeros, in ascending order, in the
lower triangular part of column jcol of the matrix;
jcol ≤ jrowin(1) < jrowin(2) < ... < jrowin(nzcol) ≤
neqns.
For a nonsymmetric matrix (LU factorization), a list of
row indices for all nonzeros, in ascending order, in
column jcol of the matrix; 1 ≤ jrowin(1) < jrowin(2) < ...
< jrowin(nzcol) ≤ neqns
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: