floating point test suite



Also, there are functions float_approx_equal and bfloat_approx_equal.
The default for float_approx_equal_tolerance is about 8 * unit roundoff.

(%i1) float_approx_equal_tolerance;
(%o1) 8.881784197001254E-16

(%i2) float_approx_equal(float(2/3), 0.66666666666666666666666666666);
(%o2) true

(%i4) float_approx_equal(float(2/3), 0.66666666666);
(%o4) false

Source code comment:

;; Return true if $float converts both a and b to floats and
;; |a - b| <= float_approx_equal_tolerance * min(|a|, |b|).
;; In all other cases, return false.


Barton