DBL_EPSILON
if (d + (DBL_EPSILON * 3.0 / 4.0) == 1.0)
if (d + (DBL_EPSILON * 1.0 / 4.0) > 1.0)
test(fma, one, one, 0x1.0p-200, 1.0 + DBL_EPSILON,
test(fma, -one, one, -0x1.0p-200, -(1.0 + DBL_EPSILON),
test(fma, one, one, -0x1.0p-200, 1.0 - DBL_EPSILON / 2,
test(fma, -one, one, 0x1.0p-200, -1.0 + DBL_EPSILON / 2,
test(fma, 0x1.0p-100, 0x1.0p-100, 1.0, 1.0 + DBL_EPSILON,
test(fma, -0x1.0p-100, 0x1.0p-100, -1.0, -(1.0 + DBL_EPSILON),
-1.0 + DBL_EPSILON / 2, ALL_STD_EXCEPT, FE_INEXACT);
test(fma, 0x1.0p-100, -0x1.0p-100, 1.0, 1.0 - DBL_EPSILON / 2,
testboth(1.0, 2.0, 1.0 + DBL_EPSILON, 0, );
testboth(1.0, -INFINITY, 1.0 - DBL_EPSILON/2, 0, );
testboth(-1.0, 2.0, -1.0 + DBL_EPSILON/2, 0, );
testboth(-1.0, -INFINITY, -1.0 - DBL_EPSILON, 0, );
test(nexttoward(1.0, 1.0 + ldbl_eps), 1.0 + DBL_EPSILON, 0);
assert(fabs(sin(d_pi_odd[i])) < 2 * DBL_EPSILON);
assert(fabs(tan(d_pi_odd[i])) < 2 * DBL_EPSILON);
assert(fabs(sin(-d_pi_odd[i])) < 2 * DBL_EPSILON);
assert(fabs(tan(-d_pi_odd[i])) < 2 * DBL_EPSILON);
assert(fabs(sin(d_pi_odd[i] * 2)) < 2 * DBL_EPSILON);
assert(fabs(tan(d_pi_odd[i] * 2)) < 2 * DBL_EPSILON);
assert(fabs(sin(-d_pi_odd[i] * 2)) < 2 * DBL_EPSILON);
assert(fabs(tan(-d_pi_odd[i] * 2)) < 2 * DBL_EPSILON);