127 SUBROUTINE zlaqhp( UPLO, N, AP, S, SCOND, AMAX, EQUED )
135 CHARACTER EQUED, UPLO
137 DOUBLE PRECISION AMAX, SCOND
140 DOUBLE PRECISION S( * )
147 DOUBLE PRECISION ONE, THRESH
148 parameter ( one = 1.0d+0, thresh = 0.1d+0 )
152 DOUBLE PRECISION CJ, LARGE, SMALL
156 DOUBLE PRECISION DLAMCH
157 EXTERNAL lsame, dlamch
173 small = dlamch(
'Safe minimum' ) / dlamch(
'Precision' )
176 IF( scond.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
THEN
185 IF( lsame( uplo,
'U' ) )
THEN
193 ap( jc+i-1 ) = cj*s( i )*ap( jc+i-1 )
195 ap( jc+j-1 ) = cj*cj*dble( ap( jc+j-1 ) )
205 ap( jc ) = cj*cj*dble( ap( jc ) )
207 ap( jc+i-j ) = cj*s( i )*ap( jc+i-j )
subroutine zlaqhp(UPLO, N, AP, S, SCOND, AMAX, EQUED)
ZLAQHP scales a Hermitian matrix stored in packed form.