161 SUBROUTINE cgbequb( M, N, KL, KU, AB, LDAB, R, C, ROWCND, COLCND,
170 INTEGER info, kl, ku, ldab, m, n
171 REAL amax, colcnd, rowcnd
175 COMPLEX ab( ldab, * )
182 parameter( one = 1.0e+0, zero = 0.0e+0 )
186 REAL bignum, rcmax, rcmin, smlnum, radix,
198 INTRINSIC abs, max, min, log,
REAL, aimag
204 cabs1( zdum ) = abs(
REAL( ZDUM ) ) + abs( aimag( zdum ) )
213 ELSE IF( n.LT.0 )
THEN
215 ELSE IF( kl.LT.0 )
THEN
217 ELSE IF( ku.LT.0 )
THEN
219 ELSE IF( ldab.LT.kl+ku+1 )
THEN
223 CALL
xerbla(
'CGBEQUB', -info )
229 IF( m.EQ.0 .OR. n.EQ.0 )
THEN
239 bignum = one / smlnum
253 DO 20 i = max(
j-ku, 1 ), min(
j+kl, m )
254 r( i ) = max( r( i ), cabs1( ab( kd+i-
j,
j ) ) )
258 IF( r( i ).GT.zero )
THEN
259 r( i ) = radix**int( log( r( i ) ) / logrdx )
268 rcmax = max( rcmax, r( i ) )
269 rcmin = min( rcmin, r( i ) )
273 IF( rcmin.EQ.zero )
THEN
278 IF( r( i ).EQ.zero )
THEN
288 r( i ) = one / min( max( r( i ), smlnum ), bignum )
293 rowcnd = max( rcmin, smlnum ) / min( rcmax, bignum )
306 DO 80 i = max(
j-ku, 1 ), min(
j+kl, m )
307 c(
j ) = max( c(
j ), cabs1( ab( kd+i-
j,
j ) )*r( i ) )
309 IF( c(
j ).GT.zero )
THEN
310 c(
j ) = radix**int( log( c(
j ) ) / logrdx )
319 rcmin = min( rcmin, c(
j ) )
320 rcmax = max( rcmax, c(
j ) )
323 IF( rcmin.EQ.zero )
THEN
328 IF( c(
j ).EQ.zero )
THEN
338 c(
j ) = one / min( max( c(
j ), smlnum ), bignum )
343 colcnd = max( rcmin, smlnum ) / min( rcmax, bignum )
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine cgbequb(M, N, KL, KU, AB, LDAB, R, C, ROWCND, COLCND, AMAX, INFO)
CGBEQUB
real function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j