123 SUBROUTINE cppt02( UPLO, N, NRHS, A, X, LDX, B, LDB, RWORK,
133 INTEGER ldb, ldx, n, nrhs
138 COMPLEX a( * ),
b( ldb, * ),
x( ldx, * )
145 parameter( zero = 0.0e+0, one = 1.0e+0 )
147 parameter( cone = ( 1.0e+0, 0.0e+0 ) )
151 REAL anorm, bnorm, eps, xnorm
167 IF( n.LE.0 .OR. nrhs.LE.0 )
THEN
175 anorm =
clanhp(
'1', uplo, n, a, rwork )
176 IF( anorm.LE.zero )
THEN
184 CALL
chpmv( uplo, n, -cone, a,
x( 1,
j ), 1, cone,
b( 1,
j ),
195 IF( xnorm.LE.zero )
THEN
198 resid = max( resid, ( ( bnorm/anorm )/xnorm )/eps )
REAL function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine cppt02(UPLO, N, NRHS, A, X, LDX, B, LDB, RWORK, RESID)
CPPT02
subroutine chpmv(UPLO, N, ALPHA, AP, X, INCX, BETA, Y, INCY)
CHPMV
REAL function scasum(N, CX, INCX)
SCASUM
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
REAL function clanhp(NORM, UPLO, N, AP, WORK)
CLANHP 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.