116 SUBROUTINE zptt02( UPLO, N, NRHS, D, E, X, LDX, B, LDB, RESID )
125 INTEGER ldb, ldx, n, nrhs
126 DOUBLE PRECISION resid
129 DOUBLE PRECISION d( * )
130 COMPLEX*16 b( ldb, * ), e( * ),
x( ldx, * )
136 DOUBLE PRECISION one, zero
137 parameter( one = 1.0d+0, zero = 0.0d+0 )
141 DOUBLE PRECISION anorm, bnorm, eps, xnorm
164 anorm =
zlanht(
'1', n, d, e )
169 IF( anorm.LE.zero )
THEN
176 CALL
zlaptm( 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 zptt02(UPLO, N, NRHS, D, E, X, LDX, B, LDB, RESID)
ZPTT02
DOUBLE PRECISION function dzasum(N, ZX, INCX)
DZASUM
DOUBLE PRECISION function zlanht(NORM, N, D, E)
ZLANHT 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
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
DOUBLE PRECISION function dlamch(CMACH)
DLAMCH
subroutine zlaptm(UPLO, N, NRHS, ALPHA, D, E, X, LDX, BETA, B, LDB)
ZLAPTM