Complex arg-trig functions
Stephen Montgomery-Smith
stephen at missouri.edu
Mon Aug 27 00:09:41 UTC 2012
On 08/26/2012 04:50 PM, Stephen Montgomery-Smith wrote:
> On 08/23/2012 10:07 PM, Stephen Montgomery-Smith wrote:
>> I just found out that the boost libraries implement the complex asin
>> function. I think their implementation is more faithful to the paper by
>> Hull et al than my implementation is. It does seem to have a BSD style
>> license. The only problem with it is that it is written in C++.
>>
> http://www.boost.org/doc/libs/1_51_0/boost/math/complex/asin.hpp
> Their asin seems to be about 10-15% faster than mine. Their error is
> slightly higher (4.5 ULP instead of 4ULP).
But the speed improvement are because they use their own log1p function.
When I replace it with the standard log1p, it goes 5% slower than mine.
More information about the freebsd-numerics
mailing list