140 SUBROUTINE claqsb( UPLO, N, KD, AB, LDAB, S, SCOND, AMAX, EQUED )
147 CHARACTER EQUED, UPLO
153 COMPLEX AB( LDAB, * )
160 parameter( one = 1.0e+0, thresh = 0.1e+0 )
164 REAL CJ, LARGE, SMALL
169 EXTERNAL lsame, slamch
185 small = slamch(
'Safe minimum' ) / slamch(
'Precision' )
188 IF( scond.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
THEN
197 IF( lsame( uplo,
'U' ) )
THEN
203 DO 10 i = max( 1, j-kd ), j
204 ab( kd+1+i-j, j ) = cj*s( i )*ab( kd+1+i-j, j )
213 DO 30 i = j, min( n, j+kd )
214 ab( 1+i-j, j ) = cj*s( i )*ab( 1+i-j, 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.