Date: Sat, 28 Jul 2012 00:36:36 -0500 From: Stephen Montgomery-Smith <stephen@missouri.edu> To: Bruce Evans <brde@optusnet.com.au> Cc: freebsd-bugs@FreeBSD.org, FreeBSD-gnats-submit@FreeBSD.org, Stephen Montgomery-Smith <stephen@FreeBSD.org> Subject: Re: bin/170206: complex arcsinh, log, etc. Message-ID: <50137A64.3060002@missouri.edu> In-Reply-To: <20120728120632.M909@besplex.bde.org> References: <201207270247.q6R2lkeR021134@wilberforce.math.missouri.edu> <20120727233939.A7820@besplex.bde.org> <5012FF06.4030501@missouri.edu> <20120728120632.M909@besplex.bde.org>
next in thread | previous in thread | raw e-mail | index | archive | help
This is a multi-part message in MIME format. --------------090903080003030709020200 Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit 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. --------------090903080003030709020200--
Want to link to this message? Use this URL: <https://mail-archive.FreeBSD.org/cgi/mid.cgi?50137A64.3060002>