135 SUBROUTINE zlaqhe( UPLO, N, A, LDA, S, SCOND, AMAX, EQUED )
143 CHARACTER equed, uplo
145 DOUBLE PRECISION amax, scond
148 DOUBLE PRECISION s( * )
149 COMPLEX*16 a( lda, * )
155 DOUBLE PRECISION one, thresh
156 parameter( one = 1.0d+0, thresh = 0.1d+0 )
160 DOUBLE PRECISION cj, large, small
181 small =
dlamch(
'Safe minimum' ) /
dlamch(
'Precision' )
184 IF( scond.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
THEN
193 IF(
lsame( uplo,
'U' ) )
THEN
200 a( i,
j ) = cj*s( i )*a( i,
j )
202 a(
j,
j ) = cj*cj*dble( a(
j,
j ) )
210 a(
j,
j ) = cj*cj*dble( a(
j,
j ) )
212 a( i,
j ) = cj*s( i )*a( i,
j )
logical function lsame(CA, CB)
LSAME
subroutine zlaqhe(UPLO, N, A, LDA, S, SCOND, AMAX, EQUED)
ZLAQHE scales a Hermitian matrix.
double precision function dlamch(CMACH)
DLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j