131 SUBROUTINE zlaqhe( UPLO, N, A, LDA, S, SCOND, AMAX, EQUED )
138 CHARACTER EQUED, UPLO
140 DOUBLE PRECISION AMAX, SCOND
143 DOUBLE PRECISION S( * )
144 COMPLEX*16 A( LDA, * )
150 DOUBLE PRECISION ONE, THRESH
151 parameter( one = 1.0d+0, thresh = 0.1d+0 )
155 DOUBLE PRECISION CJ, LARGE, SMALL
159 DOUBLE PRECISION DLAMCH
160 EXTERNAL lsame, dlamch
176 small = dlamch(
'Safe minimum' ) / dlamch(
'Precision' )
179 IF( scond.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
THEN
188 IF( lsame( uplo,
'U' ) )
THEN
195 a( i, j ) = cj*s( i )*a( i, j )
197 a( j, j ) = cj*cj*dble( a( j, j ) )
205 a( j, j ) = cj*cj*dble( a( j, j ) )
207 a( i, j ) = cj*s( i )*a( i, j )
subroutine zlaqhe(uplo, n, a, lda, s, scond, amax, equed)
ZLAQHE scales a Hermitian matrix.