139 SUBROUTINE zgbt02( TRANS, M, N, KL, KU, NRHS, A, LDA, X, LDX, B,
149 INTEGER kl, ku, lda, ldb, ldx, m, n, nrhs
150 DOUBLE PRECISION resid
153 COMPLEX*16 a( lda, * ),
b( ldb, * ), x( ldx, * )
159 DOUBLE PRECISION zero, one
160 parameter( zero = 0.0d+0, one = 1.0d+0 )
162 parameter( cone = ( 1.0d+0, 0.0d+0 ) )
165 INTEGER i1, i2,
j, kd, n1
166 DOUBLE PRECISION anorm, bnorm, eps, xnorm
183 IF( m.LE.0 .OR. n.LE.0 .OR. nrhs.LE.0 )
THEN
194 i1 = max( kd+1-
j, 1 )
195 i2 = min( kd+m-
j, kl+kd )
196 anorm = max( anorm,
dzasum( i2-i1+1, a( i1,
j ), 1 ) )
198 IF( anorm.LE.zero )
THEN
203 IF(
lsame( trans,
'T' ) .OR.
lsame( trans,
'C' ) )
THEN
212 CALL
zgbmv( trans, m, n, kl, ku, -cone, a, lda, x( 1,
j ), 1,
213 $ cone,
b( 1,
j ), 1 )
222 xnorm =
dzasum( n1, x( 1,
j ), 1 )
223 IF( xnorm.LE.zero )
THEN
226 resid = max( resid, ( ( bnorm / anorm ) / xnorm ) / eps )
subroutine zgbmv(TRANS, M, N, KL, KU, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
ZGBMV
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
logical function lsame(CA, CB)
LSAME
subroutine zgbt02(TRANS, M, N, KL, KU, NRHS, A, LDA, X, LDX, B, LDB, RESID)
ZGBT02
double precision function dzasum(N, ZX, INCX)
DZASUM
double precision function dlamch(CMACH)
DLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j