Use of C99 extra long double math functions after r236148

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


I did a ULP test on clog.  The test code is attached.  (Not the cleanest 
code, I know, but it does the job.)  It needs the mpfr and unuran ports 
installed.

To my shock, I found that under certain circumstances, the ULP in the 
real part was huge.  The problem is when hypot(x,y) is close to 1, 
because then the real part of clog is close to zero.  I was seeing ULPs 
in the thousands.

I struggled to find a solution, and now I think I have the ULP down to 
about 2.  I am going to work on it more tomorrow to see if I can get ULP 
down even further.


More information about the freebsd-numerics mailing list