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 )
184 xnorm =
scasum( n, x( 1,
j ), 1 )
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 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 b(3) integer i
subroutine cptt02(UPLO, N, NRHS, D, E, X, LDX, B, LDB, RESID)
CPTT02
real function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j