This is the mail archive of the
libc-alpha@sourceware.org
mailing list for the glibc project.
[PATCH] Use mantissa_t in __mpexp
- From: Siddhesh Poyarekar <siddhesh at redhat dot com>
- To: libc-alpha at sourceware dot org
- Date: Tue, 2 Apr 2013 17:28:02 +0530
- Subject: [PATCH] Use mantissa_t in __mpexp
Hi,
Here's a minor cleanup in mpexp to use mantissa_t instead of double.
I have also adjusted HALFRAD accordingly so that it is integral on
x86_64. Verified that this works fine on x86_64 and does not regress
the testsuite. OK to commit?
Siddhesh
* sysdeps/ieee754/dbl-64/mpa.h (HALFRAD): Define as TWO23.
* sysdeps/ieee754/dbl-64/mpexp.c (__mpexp): Use mantissa_t.
diff --git a/sysdeps/ieee754/dbl-64/mpa.h b/sysdeps/ieee754/dbl-64/mpa.h
index ee81f5c..0d82115 100644
--- a/sysdeps/ieee754/dbl-64/mpa.h
+++ b/sysdeps/ieee754/dbl-64/mpa.h
@@ -100,6 +100,8 @@ extern const mp_no mptwo;
#define TWO19 TWOPOW (19) /* 2^19 */
#define TWO23 TWOPOW (23) /* 2^23 */
+#define HALFRAD TWO23
+
#define TWO57 0x1.0p57 /* 2^57 */
#define TWO71 0x1.0p71 /* 2^71 */
#define TWOM1032 0x1.0p-1032 /* 2^-1032 */
@@ -107,7 +109,6 @@ extern const mp_no mptwo;
#define HALF 0x1.0p-1 /* 1/2 */
#define MHALF -0x1.0p-1 /* -1/2 */
-#define HALFRAD 0x1.0p23 /* 2^23 */
int __acr (const mp_no *, const mp_no *, int);
void __cpy (const mp_no *, mp_no *, int);
diff --git a/sysdeps/ieee754/dbl-64/mpexp.c b/sysdeps/ieee754/dbl-64/mpexp.c
index a219401..06347b7 100644
--- a/sysdeps/ieee754/dbl-64/mpexp.c
+++ b/sysdeps/ieee754/dbl-64/mpexp.c
@@ -43,7 +43,7 @@ SECTION
__mpexp (mp_no *x, mp_no *y, int p)
{
int i, j, k, m, m1, m2, n;
- double b;
+ mantissa_t b;
static const int np[33] =
{
0, 0, 0, 0, 3, 3, 4, 4, 5, 4, 4, 5, 5, 5, 6, 6, 6, 6, 6, 6,
@@ -80,7 +80,7 @@ __mpexp (mp_no *x, mp_no *y, int p)
b = X[1];
m2 = 24 * EX;
for (; b < HALFRAD; m2--)
- b *= TWO;
+ b *= 2;
if (b == HALFRAD)
{
for (i = 2; i <= p; i++)