HP MLIB User's Guide Vol. 2 7th Ed.
1048 HP MLIB User’s Guide
DSEVIE Matrix structure input by finite element
Output ier Status response:
ier = 0 Normal return.
ier = −100 Incorrect processing path; DSEVIN
not called or matrix input already
finished.
ier = −101 Error in dynamic storage allocation.
ier = −107 Illegal value for nnode.
ier = −108 Illegal value for at least one entry in
nodlst.
ier = −109 Illegal value for matrix designator.
ier = −110 Attempt to add nonzero location to B,
which was not specified as having a
general sparse structure (B is a
diagonal matrix, an identity matrix,
or has same structure as A).
Notes Calls to DSEVIE can be intermixed with calls to DSEVI1. DSEVIC and
DSEVIM cannot be used if DSEVI1 or DSEVIE are used.
Actual character arguments in a subroutine call may be longer than
corresponding dummy arguments. Therefore, readability of the CALL
statement may be improved by coding the matrix argument as ’addleft’.
Example Rows and columns 345, 346, 347, and 989 form a small dense submatrix of A.
Add positions consisting of all pairs of numbers from this set to the list of
nonzeros in the matrix A.
INTEGER*4 NNODE,NODLST(10),IER
REAL*8 GLOBAL(150)
NNODE = 4
NODLST(1) = 345
NODLST(2) = 346
NODLST(3) = 347
NODLST(4) = 989
CALL DSEVIE (’A’,NNODE,NODLST,GLOBAL,IER)
IF ( IER .NE. 0 ) THEN
handle error condition
ENDIF