127 SUBROUTINE zlaqsp( UPLO, N, AP, S, SCOND, AMAX, EQUED )
135 CHARACTER equed, uplo
137 DOUBLE PRECISION amax, scond
140 DOUBLE PRECISION s( * )
147 DOUBLE PRECISION one, thresh
148 parameter( one = 1.0d+0, thresh = 0.1d+0 )
152 DOUBLE PRECISION cj, large, small
170 small =
dlamch(
'Safe minimum' ) /
dlamch(
'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 )
202 ap( jc+i-
j ) = cj*s( i )*ap( jc+i-
j )
LOGICAL function lsame(CA, CB)
LSAME
subroutine zlaqsp(UPLO, N, AP, S, SCOND, AMAX, EQUED)
ZLAQSP scales a symmetric/Hermitian matrix in packed storage, using scaling factors computed by sppeq...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
DOUBLE PRECISION function dlamch(CMACH)
DLAMCH