114 SUBROUTINE zpoequ( N, A, LDA, S, SCOND, AMAX, INFO )
123 DOUBLE PRECISION AMAX, SCOND
126 DOUBLE PRECISION S( * )
127 COMPLEX*16 A( lda, * )
133 DOUBLE PRECISION ZERO, ONE
134 parameter ( zero = 0.0d+0, one = 1.0d+0 )
138 DOUBLE PRECISION SMIN
144 INTRINSIC dble, max, min, sqrt
153 ELSE IF( lda.LT.max( 1, n ) )
THEN
157 CALL xerbla(
'ZPOEQU', -info )
171 s( 1 ) = dble( a( 1, 1 ) )
175 s( i ) = dble( a( i, i ) )
176 smin = min( smin, s( i ) )
177 amax = max( amax, s( i ) )
180 IF( smin.LE.zero )
THEN
185 IF( s( i ).LE.zero )
THEN
196 s( i ) = one / sqrt( s( i ) )
201 scond = sqrt( smin ) / sqrt( amax )
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine zpoequ(N, A, LDA, S, SCOND, AMAX, INFO)
ZPOEQU