121 SUBROUTINE strt06( RCOND, RCONDC, UPLO, DIAG, N, A, LDA, WORK,
132 REAL rat, rcond, rcondc
135 REAL a( lda, * ), work( * )
142 parameter( zero = 0.0e+0, one = 1.0e+0 )
145 REAL anorm, bignum, eps, rmax, rmin, smlnum
160 rmax = max( rcond, rcondc )
161 rmin = min( rcond, rcondc )
165 IF( rmin.LT.zero )
THEN
171 ELSE IF( rmin.GT.zero )
THEN
175 rat = rmax / rmin - one
177 ELSE IF( rmax.EQ.zero )
THEN
191 smlnum =
slamch(
'Safe minimum' )
192 bignum = one / smlnum
193 CALL
slabad( smlnum, bignum )
194 anorm =
slantr(
'M', uplo, diag, n, n, a, lda, work )
196 rat = rmax*( min( bignum / max( one, anorm ), one / eps ) )
subroutine strt06(RCOND, RCONDC, UPLO, DIAG, N, A, LDA, WORK, RAT)
STRT06
real function slamch(CMACH)
SLAMCH
subroutine slabad(SMALL, LARGE)
SLABAD
real function slantr(NORM, UPLO, DIAG, M, N, A, LDA, WORK)
SLANTR returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a trapezoidal or triangular matrix.