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 )