113 DOUBLE PRECISION AYB( n, nrhs ), BERR( nrhs )
114 DOUBLE PRECISION RES( n, nrhs )
128 DOUBLE PRECISION DLAMCH
129 DOUBLE PRECISION SAFE1
137 safe1 = dlamch(
'Safe minimum' )
143 IF (ayb(i,j) .NE. 0.0d+0)
THEN
144 tmp = (safe1+abs(res(i,j)))/ayb(i,j)
145 berr(j) = max( berr(j), tmp )
double precision function dlamch(CMACH)
DLAMCH
subroutine dla_lin_berr(N, NZ, NRHS, RES, AYB, BERR)
DLA_LIN_BERR computes a component-wise relative backward error.