141 SUBROUTINE slaqsb( UPLO, N, KD, AB, LDAB, S, SCOND, AMAX, EQUED )
149 CHARACTER equed, uplo
154 REAL ab( ldab, * ), s( * )
161 parameter( one = 1.0e+0, thresh = 0.1e+0 )
165 REAL cj, large, small
186 small =
slamch(
'Safe minimum' ) /
slamch(
'Precision' )
189 IF( scond.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
THEN
198 IF(
lsame( uplo,
'U' ) )
THEN
204 DO 10 i = max( 1,
j-kd ),
j
205 ab( kd+1+i-
j,
j ) = cj*s( i )*ab( kd+1+i-
j,
j )
214 DO 30 i =
j, min( n,
j+kd )
215 ab( 1+i-
j,
j ) = cj*s( i )*ab( 1+i-
j,
j )
logical function lsame(CA, CB)
LSAME
subroutine slaqsb(UPLO, N, KD, AB, LDAB, S, SCOND, AMAX, EQUED)
SLAQSB scales a symmetric/Hermitian band matrix, using scaling factors computed by spbequ...
real function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j