108 DOUBLE PRECISION AYB( N, NRHS ), BERR( NRHS )
109 COMPLEX*16 RES( N, NRHS )
120 INTRINSIC abs, real, dimag, max
124 DOUBLE PRECISION DLAMCH
125 DOUBLE PRECISION SAFE1
131 cabs1( cdum ) = abs( dble( cdum ) ) + abs( dimag( cdum ) )
139 safe1 = dlamch(
'Safe minimum' )
145 IF (ayb(i,j) .NE. 0.0d+0)
THEN
146 tmp = (safe1 + cabs1(res(i,j)))/ayb(i,j)
147 berr(j) = max( berr(j), tmp )
subroutine zla_lin_berr(n, nz, nrhs, res, ayb, berr)
ZLA_LIN_BERR computes a component-wise relative backward error.