bin/170206: complex arcsinh, log, etc.
Stephen Montgomery-Smith
stephen at missouri.edu
Sat Jul 28 05:36:43 UTC 2012
Yes, everywhere I said "double precision" I meant "doubled precision."
I think the papers by Hull et al were perfectly happy with a ULP of
around 4.
I have been trying to do a little better, but like you I am noticing
that log1p isn't that good either.
I have tried some other things. I am attaching this example which gets
a ULP a little over 2. I simulate high precision arithmetic by
expanding everything out into integers. I certainly didn't aim for a
speedy program.
More information about the freebsd-bugs
mailing list