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

962 HP MLIB User’s Guide
DSLEVM Matrix value Input by matrix
VALUES(1) = 1.0D0
VALUES(2) = 1.0D0
VALUES(3) = 1.0D0
VALUES(4) = 1.0D0
VALUES(5) = −1.0D0
VALUES(6) = −1.0D0
VALUES(7) = 1.0D0
VALUES(8) =
−1.0D0
VALUES(9) = 1.0D0
VALUES(10) =
−1.0D0
VALUES(11) = −1.0D0
VALUES(12) = −1.0D0
CALL DSLEVM (COLSTR,ROWIND,VALUES,GLOBAL,IER)
IF ( IER .NE. 0 ) THEN
handle error condition
ENDIF
DIAGVL = 4.0
DO 100 I = 1, NEQNS
CALL DSLEV1 (I,I,DIAGVL,GLOBAL,IER)
IF ( IER .NE. 0 ) THEN
handle error condition
ENDIF
100 CONTINUE
This matrix could also have been input strictly with DSLEVM by expanding
COLSTR, ROWIND, and VALUES to include the diagonal. This action would have
avoided the calls to DSLEV1. Including the diagonal entries in COLSTR and
ROWIND is also acceptable to DSLEIM.