135 SUBROUTINE zlaqsy( 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
178 small =
dlamch(
'Safe minimum' ) /
dlamch(
'Precision' )
181 IF( scond.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
THEN
190 IF(
lsame( uplo,
'U' ) )
THEN
197 a( i,
j ) = cj*s( i )*a( i,
j )
207 a( i,
j ) = cj*s( i )*a( i,
j )
subroutine zlaqsy(UPLO, N, A, LDA, S, SCOND, AMAX, EQUED)
ZLAQSY scales a symmetric/Hermitian matrix, using scaling factors computed by spoequ.
logical function lsame(CA, CB)
LSAME
double precision function dlamch(CMACH)
DLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j