103 SUBROUTINE dget04( N, NRHS, X, LDX, XACT, LDXACT, RCOND, RESID )
111 INTEGER ldx, ldxact, n, nrhs
112 DOUBLE PRECISION rcond, resid
115 DOUBLE PRECISION x( ldx, * ), xact( ldxact, * )
121 DOUBLE PRECISION zero
122 parameter( zero = 0.0d+0 )
126 DOUBLE PRECISION diffnm, eps, xnorm
140 IF( n.LE.0 .OR. nrhs.LE.0 )
THEN
148 IF( rcond.LT.zero )
THEN
159 ix =
idamax( 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
167 $ resid = 1.0d0 / eps
169 resid = max( resid, ( diffnm / xnorm )*rcond )
172 IF( resid*eps.LT.1.0d0 )
173 $ resid = resid / eps
double precision function dlamch(CMACH)
DLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
integer function idamax(N, DX, INCX)
IDAMAX
subroutine dget04(N, NRHS, X, LDX, XACT, LDXACT, RCOND, RESID)
DGET04