Date: Sun, 29 Jul 2012 12:31:18 +1000 (EST) From: Bruce Evans <brde@optusnet.com.au> To: Stephen Montgomery-Smith <stephen@missouri.edu> Cc: freebsd-bugs@FreeBSD.org Subject: Re: bin/170206: complex arcsinh, log, etc. Message-ID: <20120729121443.M1008@besplex.bde.org> In-Reply-To: <201207281550.q6SFoBdC073014@freefall.freebsd.org> References: <201207281550.q6SFoBdC073014@freefall.freebsd.org>
next in thread | previous in thread | raw e-mail | index | archive | help
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
Want to link to this message? Use this URL: <https://mail-archive.FreeBSD.org/cgi/mid.cgi?20120729121443.M1008>