102 SUBROUTINE sstech( N, A, B, EIG, TOL, WORK, INFO )
114 REAL a( * ),
b( * ), eig( * ), work( * )
121 parameter( zero = 0.0e0 )
124 INTEGER bpnt, count, i, isub,
j, numl, numu, tpnt
125 REAL emin, eps, lower, mx, tuppr, unflep, upper
148 IF( tol.LT.zero )
THEN
156 unflep =
slamch(
'Safe minimum' ) / eps
163 mx = max( mx, abs( eig( i ) ) )
165 eps = max( eps*mx, unflep )
175 DO 30
j = 2, n + 1 - i
176 IF( work(
j ).LT.emin )
THEN
181 IF( isub.NE.n+1-i )
THEN
182 work( isub ) = work( n+1-i )
196 upper = work( tpnt ) + eps
197 lower = work( bpnt ) - eps
204 tuppr = work( bpnt+1 ) + eps
211 lower = work( bpnt ) - eps
217 CALL
sstect( n, a,
b, lower, numl )
218 CALL
sstect( n, a,
b, upper, numu )
220 IF( count.NE.bpnt-tpnt+1 )
THEN
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine sstech(N, A, B, EIG, TOL, WORK, INFO)
SSTECH
subroutine sstect(N, A, B, SHIFT, NUM)
SSTECT
real function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j