diff options
Diffstat (limited to 'lib/libm/src')
| -rw-r--r-- | lib/libm/src/b_log__D.c | 8 | ||||
| -rw-r--r-- | lib/libm/src/b_tgamma.c | 23 |
2 files changed, 20 insertions, 11 deletions
diff --git a/lib/libm/src/b_log__D.c b/lib/libm/src/b_log__D.c index 4c862a7b9d4..7e8c82046d5 100644 --- a/lib/libm/src/b_log__D.c +++ b/lib/libm/src/b_log__D.c @@ -1,4 +1,4 @@ -/* $OpenBSD: b_log__D.c,v 1.1 2008/06/11 20:53:27 martynas Exp $ */ +/* $OpenBSD: b_log__D.c,v 1.2 2008/06/12 22:15:36 martynas Exp $ */ /* * Copyright (c) 1992, 1993 * The Regents of the University of California. All rights reserved. @@ -381,8 +381,10 @@ __log__D(double x) q = u*v*(A1 + v*(A2 + v*(A3 + v*A4))); if (m | j) u1 = u + 513, u1 -= 513; - else - u1 = u, TRUNC(u1); + else { + u1 = u; + TRUNC(u1); + } u2 = (2.0*(f - F*u1) - u1*f) * g; u1 += m*logF_head[N] + logF_head[j]; diff --git a/lib/libm/src/b_tgamma.c b/lib/libm/src/b_tgamma.c index 74b36085474..cb46a0dbc16 100644 --- a/lib/libm/src/b_tgamma.c +++ b/lib/libm/src/b_tgamma.c @@ -1,4 +1,4 @@ -/* $OpenBSD: b_tgamma.c,v 1.1 2008/06/11 20:53:27 martynas Exp $ */ +/* $OpenBSD: b_tgamma.c,v 1.2 2008/06/12 22:15:36 martynas Exp $ */ /*- * Copyright (c) 1992, 1993 * The Regents of the University of California. All rights reserved. @@ -234,22 +234,26 @@ smaller_gam(double x) double t, d; struct Double r, xx; if (x < x0 + LEFT) { - t = x, TRUNC(t); + t = x; + TRUNC(t); d = (t+x)*(x-t); t *= t; - xx.a = (t + x), TRUNC(xx.a); + xx.a = (t + x); + TRUNC(xx.a); xx.b = x - xx.a; xx.b += t; xx.b += d; t = (one-x0); t += x; d = (one-x0); d -= t; d += x; x = xx.a + xx.b; } else { - xx.a = x, TRUNC(xx.a); + xx.a = x; + TRUNC(xx.a); xx.b = x - xx.a; t = x - x0; d = (-x0 -t); d += x; } r = ratfun_gam(t, d); - d = r.a/x, TRUNC(d); + d = r.a/x; + TRUNC(d); r.a -= d*xx.a; r.a -= d*xx.b; r.a += r.b; return (d + r.a/x); } @@ -269,17 +273,20 @@ ratfun_gam(double z, double c) /* return r.a + r.b = a0 + (z+c)^2*p/q, with r.a truncated to 26 bits. */ p = p/q; - t.a = z, TRUNC(t.a); /* t ~= z + c */ + t.a = z; + TRUNC(t.a); /* t ~= z + c */ t.b = (z - t.a) + c; t.b *= (t.a + z); q = (t.a *= t.a); /* t = (z+c)^2 */ TRUNC(t.a); t.b += (q - t.a); - r.a = p, TRUNC(r.a); /* r = P/Q */ + r.a = p; + TRUNC(r.a); /* r = P/Q */ r.b = p - r.a; t.b = t.b*p + t.a*r.b + a0_lo; t.a *= r.a; /* t = (z+c)^2*(P/Q) */ - r.a = t.a + a0_hi, TRUNC(r.a); + r.a = t.a + a0_hi; + TRUNC(r.a); r.b = ((a0_hi-r.a) + t.a) + t.b; return (r); /* r = a0 + t */ } |