[<prev] [next>] [thread-next>] [day] [month] [year] [list]
Message-ID: <m3y3xcypfh.fsf@jhcloos.com>
Date: Sat, 11 Feb 2017 23:59:46 -0500
From: James Cloos <cloos@...loos.com>
To: musl@...ts.openwall.com
Subject: catrig
There was a recent post on one of the freebsd lists about Stephen
Montgomery-Smith's catrig routines:
http://faculty.missouri.edu/~stephen/software/#catrig
which is based on the paper "Implementing the complex arcsine and
arccosine functions using exception handling" by T. E. Hull, Thomas
F. Fairgrieve, and Ping Tak Peter Tang, published in ACM Transactions on
Mathematical Software, Volume 23 Issue 3, 1997, Pages 299-335,
http://dl.acm.org/citation.cfm?id=275324.
It looks to be the most accurate catrig code availalbe. Stephen writes
that tests indicate that both the real and imaginary parts of the
results are good to within 4 ULP, that catrigl.c supports both 80 bit and
128 bit long double but that the float version is occasionally slightly
worse that 4 ULP.
I don't know how much interest the musl community has in better accuracy
for libm, but in case there is interest Stephen's code looks like a good
choice.
And I see musl already uses freebsd libm code for some of the complex trig.
-JimC
--
James Cloos <cloos@...loos.com> OpenPGP: 0x997A9F17ED7DAEA6
Powered by blists - more mailing lists
Confused about mailing lists and their use?
Read about mailing lists on Wikipedia
and check out these
guidelines on proper formatting of your messages.
