116 SUBROUTINE cptt02( UPLO, N, NRHS, D, E, X, LDX, B, LDB, RESID )
125 INTEGER ldb, ldx, n, nrhs
130 COMPLEX b( ldb, * ), e( * ),
x( ldx, * )
137 parameter( one = 1.0e+0, zero = 0.0e+0 )
141 REAL anorm, bnorm, eps, xnorm
164 anorm =
clanht(
'1', n, d, e )
169 IF( anorm.LE.zero )
THEN
176 CALL
claptm( uplo, n, nrhs, -one, d, e,
x, ldx, one,
b, ldb )
185 IF( xnorm.LE.zero )
THEN
188 resid = max( resid, ( ( bnorm / anorm ) / xnorm ) / eps )
subroutine claptm(UPLO, N, NRHS, ALPHA, D, E, X, LDX, BETA, B, LDB)
CLAPTM
REAL function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine cptt02(UPLO, N, NRHS, D, E, X, LDX, B, LDB, RESID)
CPTT02
REAL function scasum(N, CX, INCX)
SCASUM
REAL function clanht(NORM, N, D, E)
CLANHT returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex Hermitian tridiagonal matrix.
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j