fpequal
test("sign bits", fpequal(-0.0, -0.0) && !fpequal(0.0, -0.0));
test("NaN equality", fpequal(NAN, NAN) && NAN != NAN && vd != vd);
fpequal(d1 - d2, 0.0)
&& fpequal(-d1 + 0.0, 0.0)
&& fpequal(-d1 - d2, -0.0)
&& fpequal(-(d1 - d2), -0.0)
&& fpequal(-d1 - (-d2), 0.0));
fpequal(d / d, NAN) && fpequal(0.0 * d, NAN));