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