Age  Commit message (Collapse)  Author  Lines 

both jn and yn functions had integer overflow issues for large
and small n
to handle these issues nm1 (== n1) is used instead of n and n
in the code and some loops are changed to make sure the iteration
counter does not overflow
(another solution could be to use larger integer type or even double
but that has more size and runtime cost, on x87 loading int64_t or
even uint32_t into an fpu register is more than two times slower than
loading int32_t, and using double for n slows down iteration logic)
yn(1,0) now returns inf
posix2008 specifies that on overflow and at +0 all y0,y1,yn functions
return inf, this is not consistent with math when n<0 odd integer in yn
(eg. when x>0, yn(1,x)>inf, but historically yn(1,0) seems to be
special cased and returned inf)
some threshold values in jnf and ynf were fixed that seems to be
incorrectly copypasted from the double version


zero, one, two, half are replaced by const literals
The policy was to use the f suffix for float consts (1.0f),
but don't use suffix for long double consts (these consts
can be exactly represented as double).


thanks to the hard work of Szabolcs Nagy (nsz), identifying the best
(from correctness and license standpoint) implementations from freebsd
and openbsd and cleaning them up! musl should now fully support c99
float and long double math functions, and has nearcomplete complex
math support. tgmath should also work (fully on gcccompatible
compilers, and mostly on any c99 compiler).
based largely on commit 0376d44a890fea261506f1fc63833e7a686dca19 from
nsz's libm git repo, with some additions (dummy versions of a few
missing long double complex functions, etc.) by me.
various cleanups still need to be made, including readding (if
they're correct) some asm functions that were dropped.
