134 SUBROUTINE dlaqsy( UPLO, N, A, LDA, S, SCOND, AMAX, EQUED )
142 CHARACTER equed, uplo
144 DOUBLE PRECISION amax, scond
147 DOUBLE PRECISION a( lda, * ), s( * )
153 DOUBLE PRECISION one, thresh
154 parameter( one = 1.0d+0, thresh = 0.1d+0 )
158 DOUBLE PRECISION cj, large, small
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 )
205 a( i,
j ) = cj*s( i )*a( i,
j )
LOGICAL function lsame(CA, CB)
LSAME
subroutine dlaqsy(UPLO, N, A, LDA, S, SCOND, AMAX, EQUED)
DLAQSY scales a symmetric/Hermitian matrix, using scaling factors computed by spoequ.
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
DOUBLE PRECISION function dlamch(CMACH)
DLAMCH