HP MLIB User's Guide Vol. 2 7th Ed.
Chapter 15 Sparse Eigenvalues and Eigenvectors 1063
Return eigenvalue results DSEVRL
Output evalue List of eigenvalues.
ier Status response:
ier = 0 Normal return.
ier = −800 Incorrect processing path.
ier = −801 levalu not large enough.
ier = −803 Input error, indices out of range.
Example 1 Retrieve all of the final trust region eigenvalues.
INTEGER*4 NEVALS,NFOUND,IER
REAL*8 EVALUE(NEVALS),GLOBAL(150)
NEVALS = 100
LDEVCT = 1000
CALL DSEVRL (NEVALS,EVALUE,1,NFOUND,GLOBAL,IER)
IF ( IER .NE. 0 ) THEN
handle error condition
ENDIF
Example 2 Retrieve all of the computed eigenvalues.
INTEGER*4 NEVALS,NDISCD,NFOUND,IER
REAL*8 EVALUE(NEVALS),GLOBAL(150)
IF ( NDISCD .GT. 0 ) THEN
CALL DSEVRL (NEVALS,EVALUE,1,-NDISCD,GLOBAL,IER)
ELSE
CALL DSEVRL (NEVALS,EVALUE,1,NFOUND,GLOBAL,IER)
ENDIF
IF ( IER .NE. 0 ) THEN
handle error condition
ENDIF