103 SUBROUTINE cget04( N, NRHS, X, LDX, XACT, LDXACT, RCOND, RESID )
111 INTEGER ldx, ldxact, n, nrhs
115 COMPLEX x( ldx, * ), xact( ldxact, * )
122 parameter( zero = 0.0e+0 )
126 REAL diffnm, eps, xnorm
135 INTRINSIC abs, aimag, max, real
141 cabs1( zdum ) = abs(
REAL( ZDUM ) ) + abs( aimag( zdum ) )
147 IF( n.LE.0 .OR. nrhs.LE.0 )
THEN
155 IF( rcond.LT.zero )
THEN
166 ix =
icamax( n, xact( 1,
j ), 1 )
167 xnorm = cabs1( xact( ix,
j ) )
170 diffnm = max( diffnm, cabs1(
x( i,
j )-xact( i,
j ) ) )
172 IF( xnorm.LE.zero )
THEN
176 resid = max( resid, ( diffnm / xnorm )*rcond )
179 IF( resid*eps.LT.1.0 )
180 $ resid = resid / eps
REAL function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine cget04(N, NRHS, X, LDX, XACT, LDXACT, RCOND, RESID)
CGET04
INTEGER function icamax(N, CX, INCX)
ICAMAX