113 REAL ayb( n, nrhs ), berr( nrhs )
137 safe1 =
slamch(
'Safe minimum' )
143 IF (ayb(i,
j) .NE. 0.0)
THEN
144 tmp = (safe1+abs(res(i,
j)))/ayb(i,
j)
145 berr(
j) = max( berr(
j), tmp )
subroutine sla_lin_berr(N, NZ, NRHS, RES, AYB, BERR)
SLA_LIN_BERR computes a component-wise relative backward error.
real function slamch(CMACH)
SLAMCH
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j