Use of C99 extra long double math functions after r236148

Stephen Montgomery-Smith stephen at missouri.edu
Sun Aug 12 23:12:13 UTC 2012


So now that people pointed out that casinh is available in mpc, I no 
longer need the interface to pari.  I used mpc with 300 bit precision.

I tested out the casinh program.  After a small tweek, the worst case 
ULP is about 2.5 or 3.  Edge cases close to I do very well, with a ULP 
of about 0.5.  I am very pleased with how well it performs.

I looked at Peter Jeremy's code for catanh, and I notice that he hasn't 
yet worked to get optimal ULP, since he has concentrated on the handling 
of the edge cases.

If Peter is OK with me butting in, I'll could try to work on getting 
optimal ULP for catanh.  I think it will be easier to analyze than 
casinh, since no csqrts need to be involved.  I anticipate that the hard 
case will be when z is close to the imaginary axis, and moderately large 
(like 1e-5 + 10*I).




More information about the freebsd-numerics mailing list