126 INTEGER n, kl, ku, ncols, ldab, ldafb
129 REAL ab( ldab, * ), afb( ldafb, * )
136 REAL 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.0 )
THEN
156 rpvgrw = min( amax / umax, rpvgrw )