140 SUBROUTINE cgeequ( M, N, A, LDA, R, C, ROWCND, COLCND, AMAX,
149 INTEGER INFO, LDA, M, N
150 REAL AMAX, COLCND, ROWCND
161 parameter ( one = 1.0e+0, zero = 0.0e+0 )
165 REAL BIGNUM, RCMAX, RCMIN, SMLNUM
176 INTRINSIC abs, aimag, max, min, real
182 cabs1( zdum ) = abs(
REAL( ZDUM ) ) + abs( AIMAG( zdum ) )
191 ELSE IF( n.LT.0 )
THEN
193 ELSE IF( lda.LT.max( 1, m ) )
THEN
197 CALL xerbla(
'CGEEQU', -info )
203 IF( m.EQ.0 .OR. n.EQ.0 )
THEN
212 smlnum = slamch(
'S' )
213 bignum = one / smlnum
225 r( i ) = max( r( i ), cabs1( a( i, j ) ) )
234 rcmax = max( rcmax, r( i ) )
235 rcmin = min( rcmin, r( i ) )
239 IF( rcmin.EQ.zero )
THEN
244 IF( r( i ).EQ.zero )
THEN
254 r( i ) = one / min( max( r( i ), smlnum ), bignum )
259 rowcnd = max( rcmin, smlnum ) / min( rcmax, bignum )
273 c( j ) = max( c( j ), cabs1( a( i, j ) )*r( i ) )
282 rcmin = min( rcmin, c( j ) )
283 rcmax = max( rcmax, c( j ) )
286 IF( rcmin.EQ.zero )
THEN
291 IF( c( j ).EQ.zero )
THEN
301 c( j ) = one / min( max( c( j ), smlnum ), bignum )
306 colcnd = max( rcmin, smlnum ) / min( rcmax, bignum )
subroutine cgeequ(M, N, A, LDA, R, C, ROWCND, COLCND, AMAX, INFO)
CGEEQU
subroutine xerbla(SRNAME, INFO)
XERBLA