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
subroutine sget04(N, NRHS, X, LDX, XACT, LDXACT, RCOND, RESID)
SGET04
real function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j