On Sat, 28 Jul 2012, Stephen Montgomery-Smith wrote:
OK. This clog really seems to work.
x*x + y*y - 1 is computed with a ULP less than 0.8. The rest of the
errors seem to be due to the implementation of log1p. The ULP of the
final answer seems to be never bigger than a little over 2.
I really don't like this version. If we want bignums and slowness, we
can use lib[g]mp or even exec a mostly-interactive calculator (bc for
portability, closed-source bigware for unportability).
This version would have been useful in development to verify that the
doubled or tripled floating point was as exact as intended. Now writing
it in your favourite interactive calculator is easiest, except for the
problem of comparing the results.
Bruce
_______________________________________________
[email protected] mailing list
http://lists.freebsd.org/mailman/listinfo/freebsd-bugs
To unsubscribe, send any mail to "[email protected]"
