Lines Matching refs:Ulp
536 /* returns ulp(x) for x a 'normal' double-precision number */
538 Ulp (double x)
546 /* as ulp(x) <= y = x/2^52 < 2*ulp(x),
547 we have x + ulp(x) <= x + y <= x + 2*ulp(x),
548 therefore o(x + y) = x + ulp(x) or x + 2*ulp(x) */
551 eps = eps - x; /* ulp(x) or 2*ulp(x) */
556 /* returns the number of ulp's between a and b,
560 ulp (double a, double b)
570 return (int) ((a - b) / Ulp (a));