122 SUBROUTINE ctrt06( RCOND, RCONDC, UPLO, DIAG, N, A, LDA, RWORK,
133 REAL rat, rcond, rcondc
144 parameter( zero = 0.0e+0, one = 1.0e+0 )
147 REAL anorm, bignum, eps, rmax, rmin
159 rmax = max( rcond, rcondc )
160 rmin = min( rcond, rcondc )
164 IF( rmin.LT.zero )
THEN
170 ELSE IF( rmin.GT.zero )
THEN
174 rat = rmax / rmin - one
176 ELSE IF( rmax.EQ.zero )
THEN
190 bignum = one /
slamch(
'Safe minimum' )
191 anorm =
clantr(
'M', uplo, diag, n, n, a, lda, rwork )
193 rat = rmax*( min( bignum / max( one, anorm ), one / eps ) )
subroutine ctrt06(RCOND, RCONDC, UPLO, DIAG, N, A, LDA, RWORK, RAT)
CTRT06
real function slamch(CMACH)
SLAMCH
real function clantr(NORM, UPLO, DIAG, M, N, A, LDA, WORK)
CLANTR 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.