#
0ce946cf |
|
28-Oct-2014 |
Szabolcs Nagy <nsz@port70.net> |
math: use the rounding idiom consistently the idiomatic rounding of x is n = x + toint - toint; where toint is either 1/EPSILON (x is non-negative) or 1.5/EPSILON (x may be negative and nearest rounding mode is assumed) and EPSILON is according to the evaluation precision (the type of toint is not very important, because single precision float can represent the 1/EPSILON of ieee binary128). in case of FLT_EVAL_METHOD!=0 this avoids a useless store to double or float precision, and the long double code became cleaner with 1/LDBL_EPSILON instead of ifdefs for toint. __rem_pio2f and __rem_pio2 functions slightly changed semantics: on i386 a double-rounding is avoided so close to half-way cases may get evaluated differently eg. as sin(pi/4-eps) instead of cos(pi/4+eps)
|
#
73c870ed |
|
11-Apr-2014 |
Szabolcs Nagy <nsz@port70.net> |
math: fix aliasing violation in long double wrappers modfl and sincosl were passing long double* instead of double* to the wrapped double precision functions (on archs where long double and double have the same size). This is fixed now by using temporaries (this is not optimized to a single branch so the generated code is a bit bigger). Found by Morten Welinder.
|
#
b9cd1d4d |
|
05-Oct-2013 |
Szabolcs Nagy <nsz@port70.net> |
math: remove an unused variable from modfl
|
#
ee2ee92d |
|
02-Sep-2013 |
Szabolcs Nagy <nsz@port70.net> |
math: rewrite remainder functions (remainder, remquo, fmod, modf) * results are exact * modfl follows truncl (raises inexact flag spuriously now) * modf and modff only had cosmetic cleanup * remainder is just a wrapper around remquo now * using iterative shift+subtract for remquo and fmod * ld80 and ld128 are supported as well
|
#
fe2df4ad |
|
02-Jul-2012 |
Rich Felker <dalias@aerifal.cx> |
fix invalid implicit pointer conversion in ld64 modfl
|
#
9f58d060 |
|
28-Mar-2012 |
nsz <nsz@port70.net> |
math: fix modfl.c bug modfl(+-inf) was wrong on ld80 because the explicit msb was not taken into account during inf vs nan check
|
#
75483499 |
|
19-Mar-2012 |
nsz <nsz@port70.net> |
new modff.c code, fix nan handling in modfl
|
#
b69f695a |
|
12-Mar-2012 |
Rich Felker <dalias@aerifal.cx> |
first commit of the new libm! 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 near-complete complex math support. tgmath should also work (fully on gcc-compatible 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 re-adding (if they're correct) some asm functions that were dropped.
|