1 SUBROUTINE cget04( N, NRHS, X, LDX, XACT, LDXACT, RCOND, RESID )
8 INTEGER ldx, ldxact, n, nrhs
12 COMPLEX x( ldx, * ), xact( ldxact, * )
60 parameter( zero = 0.0e+0 )
64 REAL diffnm, eps, xnorm
70 EXTERNAL icamax, slamch
73 INTRINSIC abs, aimag,
max, real
79 cabs1( zdum ) = abs(
REAL( ZDUM ) ) + abs( aimag( zdum ) )
85 IF( n.LE.0 .OR. nrhs.LE.0 )
THEN
92 eps = slamch(
'Epsilon' )
93 IF( rcond.LT.zero )
THEN
104 ix = icamax( n, xact( 1, j ), 1 )
105 xnorm = cabs1( xact( ix, j ) )
108 diffnm =
max( diffnm, cabs1( x( i, j )-xact( i, j ) ) )
110 IF( xnorm.LE.zero )
THEN
114 resid =
max( resid, ( diffnm / xnorm )*rcond )
117 IF( resid*eps.LT.1.0 )
118 $ resid = resid / eps