153 SUBROUTINE dgbequ( M, N, KL, KU, AB, LDAB, R, C, ROWCND, COLCND,
162 INTEGER info, kl, ku, ldab, m, n
163 DOUBLE PRECISION amax, colcnd, rowcnd
166 DOUBLE PRECISION ab( ldab, * ), c( * ), r( * )
172 DOUBLE PRECISION one, zero
173 parameter( one = 1.0d+0, zero = 0.0d+0 )
177 DOUBLE PRECISION bignum, rcmax, rcmin, smlnum
187 INTRINSIC abs, max, min
196 ELSE IF( n.LT.0 )
THEN
198 ELSE IF( kl.LT.0 )
THEN
200 ELSE IF( ku.LT.0 )
THEN
202 ELSE IF( ldab.LT.kl+ku+1 )
THEN
206 CALL
xerbla(
'DGBEQU', -info )
212 IF( m.EQ.0 .OR. n.EQ.0 )
THEN
222 bignum = one / smlnum
234 DO 20 i = max(
j-ku, 1 ), min(
j+kl, m )
235 r( i ) = max( r( i ), abs( ab( kd+i-
j,
j ) ) )
244 rcmax = max( rcmax, r( i ) )
245 rcmin = min( rcmin, r( i ) )
249 IF( rcmin.EQ.zero )
THEN
254 IF( r( i ).EQ.zero )
THEN
264 r( i ) = one / min( max( r( i ), smlnum ), bignum )
269 rowcnd = max( rcmin, smlnum ) / min( rcmax, bignum )
283 DO 80 i = max(
j-ku, 1 ), min(
j+kl, m )
284 c(
j ) = max( c(
j ), abs( ab( kd+i-
j,
j ) )*r( i ) )
293 rcmin = min( rcmin, c(
j ) )
294 rcmax = max( rcmax, c(
j ) )
297 IF( rcmin.EQ.zero )
THEN
302 IF( c(
j ).EQ.zero )
THEN
312 c(
j ) = one / min( max( c(
j ), smlnum ), bignum )
317 colcnd = max( rcmin, smlnum ) / min( rcmax, bignum )
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine dgbequ(M, N, KL, KU, AB, LDAB, R, C, ROWCND, COLCND, AMAX, INFO)
DGBEQU
double precision function dlamch(CMACH)
DLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j