119 SUBROUTINE dpoequb( N, A, LDA, S, SCOND, AMAX, INFO )
128 DOUBLE PRECISION AMAX, SCOND
131 DOUBLE PRECISION A( lda, * ), S( * )
137 DOUBLE PRECISION ZERO, ONE
138 parameter ( zero = 0.0d+0, one = 1.0d+0 )
142 DOUBLE PRECISION SMIN, BASE, TMP
145 DOUBLE PRECISION DLAMCH
152 INTRINSIC max, min, sqrt, log, int
163 ELSE IF( lda.LT.max( 1, n ) )
THEN
167 CALL xerbla(
'DPOEQUB', -info )
180 tmp = -0.5d+0 / log( base )
189 smin = min( smin, s( i ) )
190 amax = max( amax, s( i ) )
193 IF( smin.LE.zero )
THEN
198 IF( s( i ).LE.zero )
THEN
209 s( i ) = base ** int( tmp * log( s( i ) ) )
214 scond = sqrt( smin ) / sqrt( amax )
subroutine dpoequb(N, A, LDA, S, SCOND, AMAX, INFO)
DPOEQUB
subroutine xerbla(SRNAME, INFO)
XERBLA