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 )