127 SUBROUTINE claqsp( UPLO, N, AP, S, SCOND, AMAX, EQUED )
135 CHARACTER equed, uplo
148 parameter( one = 1.0e+0, thresh = 0.1e+0 )
152 REAL cj, large, small
170 small =
slamch(
'Safe minimum' ) /
slamch(
'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 )
subroutine claqsp(UPLO, N, AP, S, SCOND, AMAX, EQUED)
CLAQSP scales a symmetric/Hermitian matrix in packed storage, using scaling factors computed by sppeq...
LOGICAL function lsame(CA, CB)
LSAME
REAL function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j