122 SUBROUTINE sppt02( UPLO, N, NRHS, A, X, LDX, B, LDB, RWORK,
132 INTEGER ldb, ldx, n, nrhs
136 REAL a( * ),
b( ldb, * ), rwork( * ),
x( ldx, * )
143 parameter( zero = 0.0e+0, one = 1.0e+0 )
147 REAL anorm, bnorm, eps, xnorm
163 IF( n.LE.0 .OR. nrhs.LE.0 )
THEN
171 anorm =
slansp(
'1', uplo, n, a, rwork )
172 IF( anorm.LE.zero )
THEN
180 CALL
sspmv( uplo, n, -one, a,
x( 1,
j ), 1, one,
b( 1,
j ), 1 )
188 bnorm =
sasum( n,
b( 1,
j ), 1 )
189 xnorm =
sasum( n,
x( 1,
j ), 1 )
190 IF( xnorm.LE.zero )
THEN
193 resid = max( resid, ( ( bnorm / anorm ) / xnorm ) / eps )
REAL function slansp(NORM, UPLO, N, AP, WORK)
SLANSP returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a symmetric matrix supplied in packed form.
subroutine sspmv(UPLO, N, ALPHA, AP, X, INCX, BETA, Y, INCY)
SSPMV
REAL function slamch(CMACH)
SLAMCH
subroutine sppt02(UPLO, N, NRHS, A, X, LDX, B, LDB, RWORK, RESID)
SPPT02
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
REAL function sasum(N, SX, INCX)
SASUM