126 INTEGER n, kl, ku, ncols, ldab, ldafb
129 DOUBLE PRECISION ab( ldab, * ), afb( ldafb, * )
136 DOUBLE PRECISION amax, umax, rpvgrw
139 INTRINSIC abs, max, min
149 DO i = max(
j-ku, 1 ), min(
j+kl, n )
150 amax = max( abs( ab( kd+i-
j,
j)), amax )
152 DO i = max(
j-ku, 1 ),
j
153 umax = max( abs( afb( kd+i-
j,
j ) ), umax )
155 IF ( umax /= 0.0d+0 )
THEN
156 rpvgrw = min( amax / umax, rpvgrw )
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
double precision function dla_gbrpvgrw(N, KL, KU, NCOLS, AB, LDAB, AFB, LDAFB)
DLA_GBRPVGRW computes the reciprocal pivot growth factor norm(A)/norm(U) for a general banded matrix...