| Commit message (Expand) | Author | Age | Files | Lines |
* | Implement cabsl() in terms of hypotl(). | das | 2008-03-30 | 1 | -0/+20 |
* | Implement hypotl(). This is bde's conversion of fdlibm hypot(), with minor | das | 2008-03-30 | 1 | -0/+139 |
* | Use fabs[f]() instead of bit fiddling for setting absolute values. | bde | 2008-03-30 | 2 | -4/+4 |
* | Use the expression fabs(x+0.0)-fabs(y+0.0) instead of | bde | 2008-03-30 | 2 | -2/+2 |
* | Fix a missing mask in a hi+lo decomposition. Thus bug made the extra | bde | 2008-03-30 | 1 | -1/+1 |
* | Include math.h for the fmaf() prototype. | das | 2008-03-29 | 1 | -0/+2 |
* | Fix some rather obscene code that has ambiguous if...if...else... | das | 2008-03-29 | 6 | -10/+36 |
* | MI implementation of sqrtl(). This is very slow and should | das | 2008-03-02 | 3 | -4/+168 |
* | Fix and improve some magic numbers for the "medium size" case. | bde | 2008-02-28 | 2 | -4/+4 |
* | Inline __ieee754__rem_pio2f(). On amd64 (A64) and i386 (A64), this | bde | 2008-02-25 | 4 | -0/+15 |
* | Use a temporary array instead of the arg array y[] for calling | bde | 2008-02-25 | 2 | -8/+8 |
* | Change __ieee754_rem_pio2f() to return double instead of float so that | bde | 2008-02-25 | 5 | -27/+24 |
* | Fix some off-by-1 errors. | bde | 2008-02-25 | 3 | -6/+11 |
* | Optimize the 9pi/2 < |x| <= 2**19pi/2 case some more by avoiding an | bde | 2008-02-23 | 2 | -11/+8 |
* | Avoid using FP-to-integer conversion for !(amd64 || i386) too. Use the | bde | 2008-02-22 | 2 | -6/+4 |
* | Remove the "quick check no cancellation" optimization for | bde | 2008-02-22 | 1 | -12/+1 |
* | Optimize the 9pi/2 < |x| <= 2**19pi/2 case on amd64 and i386 by avoiding | bde | 2008-02-22 | 2 | -0/+18 |
* | Add an irint() function in inline asm for amd64 and i386. irint() is | bde | 2008-02-22 | 1 | -0/+30 |
* | Optimize the conversion to bits a little (by about 11 cycles or 16% | bde | 2008-02-22 | 1 | -5/+13 |
* | Optimize the fixup for +-0 by using better classification for this case | bde | 2008-02-22 | 1 | -2/+4 |
* | Fix rintl() on signaling NaNs and unsupported formats. | bde | 2008-02-22 | 1 | -5/+3 |
* | s/rcsid/__FBSDID/ | das | 2008-02-22 | 79 | -241/+158 |
* | Remove an unused variable. | das | 2008-02-22 | 3 | -3/+3 |
* | Merge cosmetic changes from e_rem_pio2.c 1.10 (convert to __FBSDID(); | bde | 2008-02-19 | 1 | -6/+4 |
* | Optimize for 3pi/4 <= |x| <= 9pi/4 in much the same way as for | bde | 2008-02-19 | 1 | -18/+56 |
* | Rearrange the polynomial evaluation for better parallelism. This | bde | 2008-02-19 | 2 | -9/+9 |
* | Add tgammaf() as a simple wrapper around tgamma(). | das | 2008-02-18 | 2 | -0/+44 |
* | Inline __ieee754__rem_pio2(). With gcc4-2, this gives an average | bde | 2008-02-18 | 4 | -13/+19 |
* | Use volatile hacks to make sure exp() generates an underflow | das | 2008-02-17 | 1 | -1/+2 |
* | Add implementations of sinl(), cosl(), and tanl(). | das | 2008-02-17 | 7 | -0/+367 |
* | Add more pi for long doubles. Also, avoid storing multiple copies | das | 2008-02-17 | 4 | -50/+154 |
* | Sigh, the weak reference for ceill(), floorl() and truncl() was in | bde | 2008-02-15 | 3 | -6/+10 |
* | Oops, the weak reference for ceill(), floorl() and truncl() was in the | bde | 2008-02-14 | 6 | -12/+12 |
* | Use the expression fabs(x+0.0)+fabs(y+0.0) instad of a+b (where a is | bde | 2008-02-14 | 2 | -8/+8 |
* | Fix the hi+lo decomposition for 2/(3ln2). The decomposition needs to | bde | 2008-02-14 | 1 | -2/+2 |
* | Use the expression (x+0.0)-(y+0.0) instead of x+y when mixing NaN arg(s). | bde | 2008-02-14 | 2 | -10/+8 |
* | s_ceill.c | bde | 2008-02-13 | 3 | -9/+6 |
* | On arches where long double is the same as double, alias ceil(), floor() | bde | 2008-02-13 | 3 | -0/+12 |
* | Fix the C version of ceill(x) for -1 < x <= -0 in all rounding modes. | bde | 2008-02-13 | 1 | -1/+1 |
* | Fix exp2*(x) on signaling NaNs by returning x+x as usual. | bde | 2008-02-13 | 2 | -2/+2 |
* | Rearrange the polynomial evaluation for better parallelism. This is | bde | 2008-02-13 | 1 | -3/+4 |
* | Fix remainder() and remainderf() in round-towards-minus-infinity mode | bde | 2008-02-12 | 2 | -8/+8 |
* | Use double precision for z and thus for the entire calculation of | bde | 2008-02-11 | 1 | -3/+4 |
* | As usual, use a minimax polynomial that is specialized for float | bde | 2008-02-09 | 1 | -7/+8 |
* | Fix a comment about coefficients and expand a related one. | bde | 2008-02-09 | 1 | -2/+2 |
* | Fix truncl() when the result should be -0.0L. When the result is +-0.0L, | bde | 2008-02-08 | 1 | -1/+2 |
* | Oops, fix the fix in rev.1.10. logb() and logbf() were broken on | bde | 2008-02-08 | 1 | -5/+4 |
* | Use a better method of scaling by 2**k. Instead of adding to the | bde | 2008-02-07 | 2 | -26/+16 |
* | Use a better method of scaling by 2**k. Instead of adding to the | bde | 2008-02-07 | 4 | -33/+28 |
* | As for the float trig functions and logf, use a minimax polynomial | bde | 2008-02-06 | 1 | -6/+7 |