lib/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2025-09-05 14:24:06 +03:00
Code Activity

Fix determination of lower precision in __mul

Browse Source
This commit is contained in:
Siddhesh Poyarekar
2013-02-15 23:56:20 +05:30
parent daaa7713e9
commit 2d0e0f29f8
2 changed files with 9 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
ChangeLog
View File

@@ -1,3 +1,8 @@
2013-02-15 Siddhesh Poyarekar <siddhesh@redhat.com>
* sysdeps/ieee754/dbl-64/mpa.c (__mul): Fix determination of
the lower precision input.
2013-02-15 Joseph Myers <joseph@codesourcery.com> 2013-02-15 Joseph Myers <joseph@codesourcery.com>
[BZ #13550] [BZ #13550]

5
sysdeps/ieee754/dbl-64/mpa.c
View File

@@ -612,6 +612,7 @@ __mul (const mp_no *x, const mp_no *y, mp_no *z, int p)
{ {
int i, j, k, ip, ip2; int i, j, k, ip, ip2;
double u, zk; double u, zk;
const mp_no *a;
/* Is z=0? */ /* Is z=0? */
if (__glibc_unlikely (X[0] * Y[0] == ZERO)) if (__glibc_unlikely (X[0] * Y[0] == ZERO))
@@ -626,9 +627,11 @@ __mul (const mp_no *x, const mp_no *y, mp_no *z, int p)
if (X[ip2] != ZERO || Y[ip2] != ZERO) if (X[ip2] != ZERO || Y[ip2] != ZERO)
break; break;
a = X[ip2] != ZERO ? y : x;
/* ... and here, at least one of them is still zero. */ /* ... and here, at least one of them is still zero. */
for (ip = ip2; ip > 0; ip--) for (ip = ip2; ip > 0; ip--)
if (X[ip] * Y[ip] != ZERO) if (a->d[ip] != ZERO)
break; break;
/* The product looks like this for p = 3 (as an example): /* The product looks like this for p = 3 (as an example):