[Maxima] floating point test suite

Robert Dodier robert.dodier at gmail.com
Sat Aug 8 13:23:34 CDT 2009


On 8/7/09, Barton Willis <willisb at unk.edu> wrote:

> ;; 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.

Incidentally this definition is, I believe, unnecessarily strict for
quantities near zero. There are some tests which fail for some
Lisp implementations because e.g. 1e-16 isn't close enough to
1.1e-16 where these numbers are residual errors of some kind
(i.e. computed by subtracting much larger numbers which are
nearly equal).

FWIW

Robert Dodier


More information about the Maxima mailing list