125 SUBROUTINE claqhp( UPLO, N, AP, S, SCOND, AMAX, EQUED )
132 CHARACTER EQUED, UPLO
145 parameter( one = 1.0e+0, thresh = 0.1e+0 )
149 REAL CJ, LARGE, SMALL
154 EXTERNAL lsame, slamch
170 small = slamch(
'Safe minimum' ) / slamch(
'Precision' )
173 IF( scond.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
THEN
182 IF( lsame( uplo,
'U' ) )
THEN
190 ap( jc+i-1 ) = cj*s( i )*ap( jc+i-1 )
192 ap( jc+j-1 ) = cj*cj*real( ap( jc+j-1 ) )
202 ap( jc ) = cj*cj*real( ap( jc ) )
204 ap( jc+i-j ) = cj*s( i )*ap( jc+i-j )
subroutine claqhp(uplo, n, ap, s, scond, amax, equed)
CLAQHP scales a Hermitian matrix stored in packed form.