126 SUBROUTINE dlaqsp( UPLO, N, AP, S, SCOND, AMAX, EQUED )
134 CHARACTER equed, uplo
136 DOUBLE PRECISION amax, scond
139 DOUBLE PRECISION ap( * ), s( * )
145 DOUBLE PRECISION one, thresh
146 parameter( one = 1.0d+0, thresh = 0.1d+0 )
150 DOUBLE PRECISION cj, large, small
168 small =
dlamch(
'Safe minimum' ) /
dlamch(
'Precision' )
171 IF( scond.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
THEN
180 IF(
lsame( uplo,
'U' ) )
THEN
188 ap( jc+i-1 ) = cj*s( i )*ap( jc+i-1 )
200 ap( jc+i-
j ) = cj*s( i )*ap( jc+i-
j )
subroutine dlaqsp(UPLO, N, AP, S, SCOND, AMAX, EQUED)
DLAQSP scales a symmetric/Hermitian matrix in packed storage, using scaling factors computed by sppeq...
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