113 DOUBLE PRECISION AYB( n, nrhs ), BERR( nrhs )
114 COMPLEX*16 RES( n, nrhs )
125 INTRINSIC abs,
REAL, DIMAG, MAX
129 DOUBLE PRECISION DLAMCH
130 DOUBLE PRECISION SAFE1
136 cabs1( cdum ) = abs( dble( cdum ) ) + abs( dimag( cdum ) )
144 safe1 = dlamch(
'Safe minimum' )
150 IF (ayb(i,j) .NE. 0.0d+0)
THEN
151 tmp = (safe1 + cabs1(res(i,j)))/ayb(i,j)
152 berr(j) = max( berr(j), tmp )
double precision function dlamch(CMACH)
DLAMCH
subroutine zla_lin_berr(N, NZ, NRHS, RES, AYB, BERR)
ZLA_LIN_BERR computes a component-wise relative backward error.