123 SUBROUTINE zppt02( UPLO, N, NRHS, A, X, LDX, B, LDB, RWORK,
133 INTEGER ldb, ldx, n, nrhs
134 DOUBLE PRECISION resid
137 DOUBLE PRECISION rwork( * )
138 COMPLEX*16 a( * ),
b( ldb, * ), x( ldx, * )
144 DOUBLE PRECISION zero, one
145 parameter( zero = 0.0d+0, one = 1.0d+0 )
147 parameter( cone = ( 1.0d+0, 0.0d+0 ) )
151 DOUBLE PRECISION anorm, bnorm, eps, xnorm
167 IF( n.LE.0 .OR. nrhs.LE.0 )
THEN
175 anorm =
zlanhp(
'1', uplo, n, a, rwork )
176 IF( anorm.LE.zero )
THEN
184 CALL
zhpmv( uplo, n, -cone, a, x( 1,
j ), 1, cone,
b( 1,
j ),
194 xnorm =
dzasum( n, x( 1,
j ), 1 )
195 IF( xnorm.LE.zero )
THEN
198 resid = max( resid, ( ( bnorm / anorm ) / xnorm ) / eps )
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
double precision function zlanhp(NORM, UPLO, N, AP, WORK)
ZLANHP 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 matrix supplied in packed form.
double precision function dzasum(N, ZX, INCX)
DZASUM
subroutine zppt02(UPLO, N, NRHS, A, X, LDX, B, LDB, RWORK, RESID)
ZPPT02
double precision function dlamch(CMACH)
DLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine zhpmv(UPLO, N, ALPHA, AP, X, INCX, BETA, Y, INCY)
ZHPMV