110 DOUBLE PRECISION AYB( N, NRHS ), BERR( NRHS )
111 COMPLEX*16 RES( N, NRHS )
122 INTRINSIC abs, real, dimag, max
126 DOUBLE PRECISION DLAMCH
127 DOUBLE PRECISION SAFE1
133 cabs1( cdum ) = abs( dble( cdum ) ) + abs( dimag( cdum ) )
141 safe1 = dlamch(
'Safe minimum' )
147 IF (ayb(i,j) .NE. 0.0d+0)
THEN
148 tmp = (safe1 + cabs1(res(i,j)))/ayb(i,j)
149 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.