summaryrefslogtreecommitdiff
path: root/src/math/modfl.c
AgeCommit message (Collapse)AuthorFilesLines
2014-10-31math: use the rounding idiom consistentlySzabolcs Nagy1-6/+4
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)
2014-04-11math: fix aliasing violation in long double wrappersSzabolcs Nagy1-1/+6
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.
2013-10-06math: remove an unused variable from modflSzabolcs Nagy1-1/+0
2013-09-05math: rewrite remainder functions (remainder, remquo, fmod, modf)Szabolcs Nagy1-83/+34
* 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
2012-07-02fix invalid implicit pointer conversion in ld64 modflRich Felker1-1/+1
2012-03-28math: fix modfl.c bugnsz1-1/+1
modfl(+-inf) was wrong on ld80 because the explicit msb was not taken into account during inf vs nan check
2012-03-19new modff.c code, fix nan handling in modflnsz1-2/+2
2012-03-13first commit of the new libm!Rich Felker1-0/+100
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.