Eric Blake <ebb9 <at> byu.net> writes: > > Done as follows. FLT_EPSILON represents a 2-ulp error for 0.5;
Actually, FLT_EPSILON represents a 2-ulp error for 0.5f, and unless FLT_EVAL_METHOD>0, it's waaay more than 2-ulp of the double 0.5. I also botched FLT_MIN instead of DBL_MIN, now fixed. -- Eric Blake