141 SUBROUTINE dlaqsb( UPLO, N, KD, AB, LDAB, S, SCOND, AMAX, EQUED )
149 CHARACTER equed, uplo
151 DOUBLE PRECISION amax, scond
154 DOUBLE PRECISION ab( ldab, * ), s( * )
160 DOUBLE PRECISION one, thresh
161 parameter( one = 1.0d+0, thresh = 0.1d+0 )
165 DOUBLE PRECISION cj, large, small
186 small =
dlamch(
'Safe minimum' ) /
dlamch(
'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 dlaqsb(UPLO, N, KD, AB, LDAB, S, SCOND, AMAX, EQUED)
DLAQSB scales a symmetric/Hermitian band matrix, using scaling factors computed by spbequ...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
DOUBLE PRECISION function dlamch(CMACH)
DLAMCH