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

912 HP MLIB User’s Guide
DSLEIC Matrix structure input by column
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 = −103 Illegal value for nzcol.
ier = −104 Illegal value for jcol, or out of order.
ier = −106 Illegal value for at least one entry in
jrowin or entries out of order.
Notes The entire matrix structure must be input with DSLEIC if it is used. Its use is
not compatible with DSELE1, DSLEIE, or DSLEIM.
If the matrix entries are available by column, using DSLEIC is more efficient
than using DSLEI1.
Example 1 Column 4519 of a symmetric matrix has entries in rows 1, 2735, 4519, 4520,
4521, 6000, 6002 and 6004. Add all five entries in the lower triangular part to
the list of nonzeros in the matrix. Columns 1 to 4518 already have been input to
the package using DSLEIC.
INTEGER*4 J,NZCOL,JROWIN(100),IER
REAL*8 GLOBAL(150)
J = 4519
NZCOL = 5
JROWIN(1) = 4520
JROWIN(2) = 4521
JROWIN(3) = 6000
JROWIN(4) = 6002
JROWIN(5) = 6004
CALL DSLEIC (J,NZCOL,JROWIN,GLOBAL,IER)
IF ( IER .NE. 0 ) THEN
handle error condition
ENDIF
Example 2 Column 4519 of a nonsymmetric matrix has entries in rows 1, 2735, 4519, 4520,
4521, 6000, 6002 and 6004. Add all seven entries to the list of nonzeros in the
matrix. Columns 1 to 4518 already have been input to the package using
DSLEIC.