139 SUBROUTINE cgbt02( TRANS, M, N, KL, KU, NRHS, A, LDA, X, LDX, B,
149 INTEGER kl, ku, lda, ldb, ldx, m, n, nrhs
153 COMPLEX a( lda, * ),
b( ldb, * ),
x( ldx, * )
160 parameter( zero = 0.0e+0, one = 1.0e+0 )
162 parameter( cone = ( 1.0e+0, 0.0e+0 ) )
165 INTEGER i1, i2,
j, kd, n1
166 REAL 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,
scasum( 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
cgbmv( trans, m, n, kl, ku, -cone, a, lda,
x( 1,
j ), 1,
213 $ cone,
b( 1,
j ), 1 )
223 IF( xnorm.LE.zero )
THEN
226 resid = max( resid, ( ( bnorm/anorm )/xnorm )/eps )
LOGICAL function lsame(CA, CB)
LSAME
subroutine cgbmv(TRANS, M, N, KL, KU, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
CGBMV
REAL function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
REAL function scasum(N, CX, INCX)
SCASUM
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine cgbt02(TRANS, M, N, KL, KU, NRHS, A, LDA, X, LDX, B, LDB, RESID)
CGBT02