113 REAL ayb( n, nrhs ), berr( nrhs )
114 COMPLEX res( n, nrhs )
125 INTRINSIC abs,
REAL, aimag, max
136 cabs1( cdum ) = abs(
REAL( CDUM ) ) + abs( aimag( cdum ) )
144 safe1 =
slamch(
'Safe minimum' )
150 IF (ayb(i,
j) .NE. 0.0)
THEN
151 tmp = (safe1 + cabs1(res(i,
j)))/ayb(i,
j)
152 berr(
j) = max( berr(
j), tmp )
real function slamch(CMACH)
SLAMCH
subroutine cla_lin_berr(N, NZ, NRHS, RES, AYB, BERR)
CLA_LIN_BERR computes a component-wise relative backward error.
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j