113 DOUBLE PRECISION ayb( n, nrhs ), berr( nrhs )
114 DOUBLE PRECISION res( n, nrhs )
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 )
subroutine dla_lin_berr(N, NZ, NRHS, RES, AYB, BERR)
DLA_LIN_BERR computes a component-wise relative backward error.
double precision function dlamch(CMACH)
DLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j