svn commit: r251339 - in head/lib/msun: ld128 ld80

[Steve Kargl]

Author: kargl
Date: Mon Jun  3 19:13:44 2013
New Revision: 251339
URL: http://svnweb.freebsd.org/changeset/base/251339

Log:
  ld80/s_expl.c:
  
  * Use integral numerical constants, and let the compiler do the
    conversion to long double.
  
  ld128/s_expl.c:
  
  * Use integral numerical constants, and let the compiler do the
    conversion to long double.
  * Use the ENTERI/RETURNI macros, which are no-ops on ld128.  This
    however makes the ld80 and ld128 identical.
  
  Reviewed by:	bde (as part of larger diff)

Modified:
  head/lib/msun/ld128/s_expl.c
  head/lib/msun/ld80/s_expl.c

Modified: head/lib/msun/ld128/s_expl.c
==============================================================================
--- head/lib/msun/ld128/s_expl.c	Mon Jun  3 18:57:35 2013	(r251338)
+++ head/lib/msun/ld128/s_expl.c	Mon Jun  3 19:13:44 2013	(r251339)
@@ -228,7 +228,7 @@ expl(long double x)
 	double dr, fn, r2;
 
 	int k, n, n2;
-	uint32_t hx, ix;
+	uint16_t hx, ix;
 
 	/* Filter out exceptional cases. */
 	u.e = x;
@@ -248,6 +248,8 @@ expl(long double x)
 		return (1 + x);		/* 1 with inexact iff x != 0 */
 	}
 
+	ENTERI();
+
 	/* Reduce x to (k*ln2 + endpoint[n2] + r1 + r2). */
 	/* Use a specialized rint() to get fn.  Assume round-to-nearest. */
 	/* XXX assume no extra precision for the additions, as for trig fns. */
@@ -262,10 +264,11 @@ expl(long double x)
 	k = n >> LOG2_INTERVALS;
 	r1 = x - fn * L1;
 	r2 = fn * -L2;
+	r = r1 + r2;
 
 	/* Prepare scale factors. */
-	v.xbits.manh = 0;
-	v.xbits.manl = 0;
+	/* XXX sparc64 multiplication is so slow that scalbnl() is faster. */
+	v.e = 1;
 	if (k >= LDBL_MIN_EXP) {
 		v.xbits.expsign = BIAS + k;
 		twopk = v.e;
@@ -284,9 +287,9 @@ expl(long double x)
 	/* Scale by 2**k. */
 	if (k >= LDBL_MIN_EXP) {
 		if (k == LDBL_MAX_EXP)
-			return (t * 2.0L * 0x1p16383L);
-		return (t * twopk);
+			RETURNI(t * 2 * 0x1p16383L);
+		RETURNI(t * twopk);
 	} else {
-		return (t * twopkp10000 * twom10000);
+		RETURNI(t * twopkp10000 * twom10000);
 	}
 }

Modified: head/lib/msun/ld80/s_expl.c
==============================================================================
--- head/lib/msun/ld80/s_expl.c	Mon Jun  3 18:57:35 2013	(r251338)
+++ head/lib/msun/ld80/s_expl.c	Mon Jun  3 19:13:44 2013	(r251339)
@@ -278,7 +278,7 @@ expl(long double x)
 	r2 = fn * -L2;
 
 	/* Prepare scale factors. */
-	v.e = 1ULL << 63;
+	v.e = 1;
 	if (k >= LDBL_MIN_EXP) {
 		v.xbits.expsign = BIAS + k;
 		twopk = v.e;
@@ -296,7 +296,7 @@ expl(long double x)
 	/* Scale by 2**k. */
 	if (k >= LDBL_MIN_EXP) {
 		if (k == LDBL_MAX_EXP)
-			RETURNI(t * 2.0L * 0x1p16383L);
+			RETURNI(t * 2 * 0x1p16383L);
 		RETURNI(t * twopk);
 	} else {
 		RETURNI(t * twopkp10000 * twom10000);