112 SUBROUTINE stpt06( RCOND, RCONDC, UPLO, DIAG, N, AP, WORK, RAT )
122 REAL rat, rcond, rcondc
125 REAL ap( * ), work( * )
132 parameter( zero = 0.0e+0, one = 1.0e+0 )
135 REAL anorm, bignum, eps, rmax, rmin, smlnum
150 rmax = max( rcond, rcondc )
151 rmin = min( rcond, rcondc )
155 IF( rmin.LT.zero )
THEN
161 ELSE IF( rmin.GT.zero )
THEN
165 rat = rmax / rmin - one
167 ELSE IF( rmax.EQ.zero )
THEN
181 smlnum =
slamch(
'Safe minimum' )
182 bignum = one / smlnum
183 CALL
slabad( smlnum, bignum )
184 anorm =
slantp(
'M', uplo, diag, n, ap, work )
186 rat = rmax*( min( bignum / max( one, anorm ), one / eps ) )
REAL function slamch(CMACH)
SLAMCH
subroutine slabad(SMALL, LARGE)
SLABAD
subroutine stpt06(RCOND, RCONDC, UPLO, DIAG, N, AP, WORK, RAT)
STPT06
REAL function slantp(NORM, UPLO, DIAG, N, AP, WORK)
SLANTP returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a triangular matrix supplied in packed form.