106 DOUBLE PRECISION FUNCTION dsxt1( IJOB, D1, N1, D2, N2, ABSTOL,
116 DOUBLE PRECISION abstol, ulp, unfl
119 DOUBLE PRECISION d1( * ), d2( * )
125 DOUBLE PRECISION zero
126 parameter( zero = 0.0d0 )
130 DOUBLE PRECISION temp1, temp2
133 INTRINSIC abs, max, min
142 IF( d2(
j ).LT.d1( i ) .AND.
j.LT.n2 )
THEN
147 temp2 = abs( d2(
j )-d1( i ) )
149 $ temp2 = temp2 / max( unfl, abstol+ulp*abs( d1( i ) ) )
151 temp2 = min( abs( d2(
j )-d1( i ) ),
152 $ abs( d1( i )-d2(
j-1 ) ) )
154 $ temp2 = temp2 / max( unfl, abstol+ulp*abs( d1( i ) ) )
156 temp1 = max( temp1, temp2 )
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
DOUBLE PRECISION function dsxt1(IJOB, D1, N1, D2, N2, ABSTOL, ULP, UNFL)
DSXT1