135 SUBROUTINE sgeequ( M, N, A, LDA, R, C, ROWCND, COLCND, AMAX,
143 INTEGER INFO, LDA, M, N
144 REAL AMAX, COLCND, ROWCND
147 REAL A( LDA, * ), C( * ), R( * )
154 parameter( one = 1.0e+0, zero = 0.0e+0 )
158 REAL BIGNUM, RCMAX, RCMIN, SMLNUM
168 INTRINSIC abs, max, min
177 ELSE IF( n.LT.0 )
THEN
179 ELSE IF( lda.LT.max( 1, m ) )
THEN
183 CALL xerbla(
'SGEEQU', -info )
189 IF( m.EQ.0 .OR. n.EQ.0 )
THEN
198 smlnum = slamch(
'S' )
199 bignum = one / smlnum
211 r( i ) = max( r( i ), abs( a( i, j ) ) )
220 rcmax = max( rcmax, r( i ) )
221 rcmin = min( rcmin, r( i ) )
225 IF( rcmin.EQ.zero )
THEN
230 IF( r( i ).EQ.zero )
THEN
240 r( i ) = one / min( max( r( i ), smlnum ), bignum )
245 rowcnd = max( rcmin, smlnum ) / min( rcmax, bignum )
259 c( j ) = max( c( j ), abs( a( i, j ) )*r( i ) )
268 rcmin = min( rcmin, c( j ) )
269 rcmax = max( rcmax, c( j ) )
272 IF( rcmin.EQ.zero )
THEN
277 IF( c( j ).EQ.zero )
THEN
287 c( j ) = one / min( max( c( j ), smlnum ), bignum )
292 colcnd = max( rcmin, smlnum ) / min( rcmax, bignum )
subroutine sgeequ(m, n, a, lda, r, c, rowcnd, colcnd, amax, info)
SGEEQU