124 SUBROUTINE zgtt02( TRANS, N, NRHS, DL, D, DU, X, LDX, B, LDB,
134 INTEGER ldb, ldx, n, nrhs
135 DOUBLE PRECISION resid
138 COMPLEX*16 b( ldb, * ), d( * ), dl( * ), du( * ),
145 DOUBLE PRECISION one, zero
146 parameter( one = 1.0d+0, zero = 0.0d+0 )
150 DOUBLE PRECISION anorm, bnorm, eps, xnorm
168 IF( n.LE.0 .OR. nrhs.EQ.0 )
174 IF(
lsame( trans,
'N' ) )
THEN
175 anorm =
zlangt(
'1', n, dl, d, du )
177 anorm =
zlangt(
'I', n, dl, d, du )
183 IF( anorm.LE.zero )
THEN
190 CALL
zlagtm( trans, n, nrhs, -one, dl, d, du,
x, ldx, one,
b,
196 IF( xnorm.LE.zero )
THEN
199 resid = max( resid, ( ( bnorm / anorm ) / xnorm ) / eps )
subroutine zlagtm(TRANS, N, NRHS, ALPHA, DL, D, DU, X, LDX, BETA, B, LDB)
ZLAGTM performs a matrix-matrix product of the form C = αAB+βC, where A is a tridiagonal matrix...
LOGICAL function lsame(CA, CB)
LSAME
DOUBLE PRECISION function dzasum(N, ZX, INCX)
DZASUM
subroutine zgtt02(TRANS, N, NRHS, DL, D, DU, X, LDX, B, LDB, RESID)
ZGTT02
DOUBLE PRECISION function zlangt(NORM, N, DL, D, DU)
ZLANGT returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value ...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
DOUBLE PRECISION function dlamch(CMACH)
DLAMCH