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

1028 HP MLIB User’s Guide
DSEVE1 One-call usage
warnng = 1xx Interval had to be expanded because
one or both endpoints were very close
to eigenvalues.
Notes 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 pbtype argument as ’vibration’ or
’buckling’.
Example In a small (order 6) sparse generalized eigenvalue problem, the matrices A and
B are as follows:
The eigenproblem is to compute the two eigenvalues nearest zero.
INTEGER*4 MSGLVL,OUTPUT,ACOLST(7),AROWIN(13),IER
REAL*8 GLOBAL(150),AVALUE(13),BVALUE(13)
MSGLVL = 1
OUTPUT = 6
ACOLST(1) = 1
ACOLST(2) = 4
ACOLST(3) = 7
ACOLST(4) = 9
ACOLST(5) = 12
ACOLST(6) = 13
ACOLST(7) = 14
AROWIN(1) = 1
AROWIN(2) = 3
AROWIN(3) = 4
AROWIN(4) = 2
AROWIN(5) = 3
AROWIN(6) = 5
A =
4.0 0.0 1.0 1.0 0.0 0.0
0.0 4.0 1.0 0.0 1.0 0.0
1.0 1.0 4.0 0.0 1.0 0.0
1.0 0.0 0.0 4.0 1.0 1.0
0.0 1.0 1.0 1.0 4.0 0.0
0.0 0.0 0.0 1.0 0.0 1.0
B =
1.0 0.0 1.0 1.0 0.0 0.0
0.0 1.0 1.0 0.0 1.0 0.0
1.0 1.0 1.0 0.0 1.0 0.0
1.0 0.0 0.0 1.0 1.0 1.0
0.0 1.0 1.0 1.0 1.0 0.0
0.0 0.0 0.0 1.0 0.0 1.0