109 SUBROUTINE dpoequ( N, A, LDA, S, SCOND, AMAX, INFO )
117 DOUBLE PRECISION AMAX, SCOND
120 DOUBLE PRECISION A( LDA, * ), S( * )
126 DOUBLE PRECISION ZERO, ONE
127 parameter( zero = 0.0d+0, one = 1.0d+0 )
131 DOUBLE PRECISION SMIN
137 INTRINSIC max, min, sqrt
146 ELSE IF( lda.LT.max( 1, n ) )
THEN
150 CALL xerbla(
'DPOEQU', -info )
169 smin = min( smin, s( i ) )
170 amax = max( amax, s( i ) )
173 IF( smin.LE.zero )
THEN
178 IF( s( i ).LE.zero )
THEN
189 s( i ) = one / sqrt( s( i ) )
194 scond = sqrt( smin ) / sqrt( amax )
subroutine dpoequ(n, a, lda, s, scond, amax, info)
DPOEQU