mirror of
https://github.com/freebsd/freebsd-src
synced 2024-10-06 16:40:47 +00:00
e38f230827
Undo the 80-bit "stub" implementation of the 128-bit long double tgammal(3) function. The latest (as of Feb 2024) version of the src/contrib/arm-optimised-routines library includes a standalone, full 128-bit replacement. This needs a small bit of wrapping to fit it in, but is otherwise a drop-in replacement. Testing this is hard, as most maths packages blow up as soon as their 80-bit floating-point capability is exceeded. With 128-bit tgammal(), this is easy to do, and this is the range that needs to be checked the most carefully. Using my copy of Maple, I was able to check that the output was within a few ULP of the correct answer, right up to the point of 128-bit over- and underflow. Additionally, the results are no worse, and indeed better than the 80-bit version. Steve Kargl sent me his libm testing code, which I used to verify that the excpetions for certain key values were correct. Tested in this case were +-Inf, +-NaN, +-1 and +-0. Differential Revision: https://reviews.freebsd.org/D44168 Reviewed by: theraven, andrew, imp |
||
---|---|---|
.. | ||
b_tgammal.c | ||
e_lgammal_r.c | ||
e_powl.c | ||
e_rem_pio2l.h | ||
invtrig.c | ||
invtrig.h | ||
k_cosl.c | ||
k_cospil.h | ||
k_expl.h | ||
k_sinl.c | ||
k_sinpil.h | ||
k_tanl.c | ||
s_cexpl.c | ||
s_cospil.c | ||
s_erfl.c | ||
s_exp2l.c | ||
s_expl.c | ||
s_logl.c | ||
s_nanl.c | ||
s_sinpil.c | ||
s_tanpil.c |