160 SUBROUTINE dgbequb( M, N, KL, KU, AB, LDAB, R, C, ROWCND, COLCND,
169 INTEGER info, kl, ku, ldab, m, n
170 DOUBLE PRECISION amax, colcnd, rowcnd
173 DOUBLE PRECISION ab( ldab, * ), c( * ), r( * )
179 DOUBLE PRECISION one, zero
180 parameter( one = 1.0d+0, zero = 0.0d+0 )
184 DOUBLE PRECISION bignum, rcmax, rcmin, smlnum, radix, logrdx
194 INTRINSIC abs, max, min, log
203 ELSE IF( n.LT.0 )
THEN
205 ELSE IF( kl.LT.0 )
THEN
207 ELSE IF( ku.LT.0 )
THEN
209 ELSE IF( ldab.LT.kl+ku+1 )
THEN
213 CALL
xerbla(
'DGBEQUB', -info )
219 IF( m.EQ.0 .OR. n.EQ.0 )
THEN
229 bignum = one / smlnum
243 DO 20 i = max( j-ku, 1 ), min( j+kl, m )
244 r( i ) = max( r( i ), abs( ab( kd+i-j, j ) ) )
248 IF( r( i ).GT.zero )
THEN
249 r( i ) = radix**int( log( r( i ) ) / logrdx )
258 rcmax = max( rcmax, r( i ) )
259 rcmin = min( rcmin, r( i ) )
263 IF( rcmin.EQ.zero )
THEN
268 IF( r( i ).EQ.zero )
THEN
278 r( i ) = one / min( max( r( i ), smlnum ), bignum )
283 rowcnd = max( rcmin, smlnum ) / min( rcmax, bignum )
296 DO 80 i = max( j-ku, 1 ), min( j+kl, m )
297 c( j ) = max( c( j ), abs( ab( kd+i-j, j ) )*r( i ) )
299 IF( c( j ).GT.zero )
THEN
300 c( j ) = radix**int( log( c( j ) ) / logrdx )
309 rcmin = min( rcmin, c( j ) )
310 rcmax = max( rcmax, c( j ) )
313 IF( rcmin.EQ.zero )
THEN
318 IF( c( j ).EQ.zero )
THEN
328 c( j ) = one / min( max( c( j ), smlnum ), bignum )
333 colcnd = max( rcmin, smlnum ) / min( rcmax, bignum )