126 INTEGER n, kl, ku, ncols, ldab, ldafb
129 COMPLEX ab( ldab, * ), afb( ldafb, * )
136 REAL amax, umax, rpvgrw
140 INTRINSIC abs, max, min,
REAL, aimag
146 cabs1( zdum ) = abs(
REAL( ZDUM ) ) + abs( aimag( zdum ) )
156 DO i = max(
j-ku, 1 ), min(
j+kl, n )
157 amax = max( cabs1( ab( kd+i-
j,
j ) ), amax )
159 DO i = max(
j-ku, 1 ),
j
160 umax = max( cabs1( afb( kd+i-
j,
j ) ), umax )
162 IF ( umax /= 0.0 )
THEN
163 rpvgrw = min( amax / umax, rpvgrw )
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
real function cla_gbrpvgrw(N, KL, KU, NCOLS, AB, LDAB, AFB, LDAFB)
CLA_GBRPVGRW computes the reciprocal pivot growth factor norm(A)/norm(U) for a general banded matrix...