103 SUBROUTINE sget04( N, NRHS, X, LDX, XACT, LDXACT, RCOND, RESID )
111 INTEGER ldx, ldxact, n, nrhs
115 REAL x( ldx, * ), xact( ldxact, * )
122 parameter( zero = 0.0e+0 )
126 REAL diffnm, eps, xnorm
140 IF( n.LE.0 .OR. nrhs.LE.0 )
THEN
148 IF( rcond.LT.zero )
THEN
159 ix =
isamax( n, xact( 1,
j ), 1 )
160 xnorm = abs( xact( ix,
j ) )
163 diffnm = max( diffnm, abs(
x( i,
j )-xact( i,
j ) ) )
165 IF( xnorm.LE.zero )
THEN
169 resid = max( resid, ( diffnm / xnorm )*rcond )
172 IF( resid*eps.LT.1.0 )
173 $ resid = resid / eps
INTEGER function isamax(N, SX, INCX)
ISAMAX
REAL function slamch(CMACH)
SLAMCH
subroutine sget04(N, NRHS, X, LDX, XACT, LDXACT, RCOND, RESID)
SGET04
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j