110 SUBROUTINE stpt06( RCOND, RCONDC, UPLO, DIAG, N, AP, WORK, RAT )
119 REAL RAT, RCOND, RCONDC
122 REAL AP( * ), WORK( * )
129 parameter( zero = 0.0e+0, one = 1.0e+0 )
132 REAL ANORM, BIGNUM, EPS, RMAX, RMIN, SMLNUM
136 EXTERNAL slamch, slantp
146 eps = slamch(
'Epsilon' )
147 rmax = max( rcond, rcondc )
148 rmin = min( rcond, rcondc )
152 IF( rmin.LT.zero )
THEN
158 ELSE IF( rmin.GT.zero )
THEN
162 rat = rmax / rmin - one
164 ELSE IF( rmax.EQ.zero )
THEN
178 smlnum = slamch(
'Safe minimum' )
179 bignum = one / smlnum
180 CALL slabad( smlnum, bignum )
181 anorm = slantp(
'M', uplo, diag, n, ap, work )
183 rat = rmax*( min( bignum / max( one, anorm ), one / eps ) )
subroutine slabad(SMALL, LARGE)
SLABAD
subroutine stpt06(RCOND, RCONDC, UPLO, DIAG, N, AP, WORK, RAT)
STPT06