126 INTEGER n, kl, ku, ncols, ldab, ldafb
129 COMPLEX*16 ab( ldab, * ), afb( ldafb, * )
136 DOUBLE PRECISION amax, umax, rpvgrw
140 INTRINSIC abs, max, min,
REAL, dimag
143 DOUBLE PRECISION cabs1
146 cabs1( zdum ) = abs( dble( zdum ) ) + abs( dimag( 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.0d+0 )
THEN
163 rpvgrw = min( amax / umax, rpvgrw )
DOUBLE PRECISION function zla_gbrpvgrw(N, KL, KU, NCOLS, AB, LDAB, AFB, LDAFB)
ZLA_GBRPVGRW computes the reciprocal pivot growth factor norm(A)/norm(U) for a general banded matrix...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j