126 SUBROUTINE cgbt01( M, N, KL, KU, A, LDA, AFAC, LDAFAC, IPIV, WORK,
135 INTEGER kl, ku, lda, ldafac, m, n
140 COMPLEX a( lda, * ), afac( ldafac, * ), work( * )
147 parameter( zero = 0.0e+0, one = 1.0e+0 )
150 INTEGER i, i1, i2, il, ip, iw,
j, jl, ju, jua, kd, lenj
162 INTRINSIC cmplx, max, min, real
169 IF( m.LE.0 .OR. n.LE.0 )
178 i1 = max( kd+1-
j, 1 )
179 i2 = min( kd+m-
j, kl+kd )
181 $ anorm = max( anorm,
scasum( i2-i1+1, a( i1,
j ), 1 ) )
191 ju = min( kl+ku,
j-1 )
193 lenj = min( m,
j ) -
j + ju + 1
195 CALL
ccopy( lenj, afac( kd-ju,
j ), 1, work, 1 )
196 DO 20 i = lenj + 1, ju + jl + 1
203 DO 30 i = min( m-1,
j ),
j - ju, -1
208 CALL
caxpy( il, t, afac( kd+1, i ), 1, work( iw+1 ),
213 work( iw ) = work( ip )
223 $ CALL
caxpy( jua+jl+1, -cmplx( one ), a( ku+1-jua,
j ), 1,
224 $ work( ju+1-jua ), 1 )
228 resid = max( resid,
scasum( ju+jl+1, work, 1 ) )
234 IF( anorm.LE.zero )
THEN
238 resid = ( ( resid /
REAL( N ) ) / anorm ) / eps
subroutine caxpy(N, CA, CX, INCX, CY, INCY)
CAXPY
real function scasum(N, CX, INCX)
SCASUM
subroutine cgbt01(M, N, KL, KU, A, LDA, AFAC, LDAFAC, IPIV, WORK, RESID)
CGBT01
real function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine ccopy(N, CX, INCX, CY, INCY)
CCOPY