142 SUBROUTINE claqsb( UPLO, N, KD, AB, LDAB, S, SCOND, AMAX, EQUED )
150 CHARACTER equed, uplo
156 COMPLEX ab( ldab, * )
163 parameter( one = 1.0e+0, thresh = 0.1e+0 )
167 REAL cj, large, small
188 small =
slamch(
'Safe minimum' ) /
slamch(
'Precision' )
191 IF( scond.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
THEN
200 IF(
lsame( uplo,
'U' ) )
THEN
206 DO 10 i = max( 1,
j-kd ),
j
207 ab( kd+1+i-
j,
j ) = cj*s( i )*ab( kd+1+i-
j,
j )
216 DO 30 i =
j, min( n,
j+kd )
217 ab( 1+i-
j,
j ) = cj*s( i )*ab( 1+i-
j,
j )
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
subroutine claqsb(UPLO, N, KD, AB, LDAB, S, SCOND, AMAX, EQUED)
CLAQSB scales a symmetric/Hermitian band matrix, using scaling factors computed by spbequ...