aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorSiddhesh Poyarekar <siddhesh@redhat.com>2013-02-16 00:15:57 +0530
committerSiddhesh Poyarekar <siddhesh@redhat.com>2013-02-16 00:15:57 +0530
commit22af19f9fb1909690dbe9840c0de093b656b4a68 (patch)
treea21f036eb56c8c5ca7ee00a6bc3f66d565e375ba
parent4709fe7602b56e9f6ee1ab6afb4067409a784f29 (diff)
downloadglibc-22af19f9fb1909690dbe9840c0de093b656b4a68.zip
glibc-22af19f9fb1909690dbe9840c0de093b656b4a68.tar.gz
glibc-22af19f9fb1909690dbe9840c0de093b656b4a68.tar.bz2
Don't require LIM to determine loop end in __sqr
-rw-r--r--ChangeLog3
-rw-r--r--sysdeps/ieee754/dbl-64/mpa.c14
2 files changed, 7 insertions, 10 deletions
diff --git a/ChangeLog b/ChangeLog
index 607433b..cfddab0 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,5 +1,8 @@
2013-02-16 Siddhesh Poyarekar <siddhesh@redhat.com>
+ * sysdeps/ieee754/dbl-64/mpa.c (__sqr): Avoid using LIM in
+ loop termination condition.
+
* sysdeps/ieee754/dbl-64/mpa.c (__mul): Use intermediate
variable to calculate EZ.
(__sqr): Likewise.
diff --git a/sysdeps/ieee754/dbl-64/mpa.c b/sysdeps/ieee754/dbl-64/mpa.c
index 9da4257..1f60fa5 100644
--- a/sysdeps/ieee754/dbl-64/mpa.c
+++ b/sysdeps/ieee754/dbl-64/mpa.c
@@ -747,12 +747,9 @@ __sqr (const mp_no *x, mp_no *y, int p)
long lim = k / 2;
if (k % 2 == 0)
- {
- yk += X[lim] * X[lim];
- lim--;
- }
+ yk += X[lim] * X[lim];
- for (i = k - p, j = p; i <= lim; i++, j--)
+ for (i = k - p, j = p; i < j; i++, j--)
yk2 += X[i] * X[j];
yk += 2.0 * yk2;
@@ -770,12 +767,9 @@ __sqr (const mp_no *x, mp_no *y, int p)
long lim = k / 2;
if (k % 2 == 0)
- {
- yk += X[lim] * X[lim];
- lim--;
- }
+ yk += X[lim] * X[lim];
- for (i = 1, j = k - 1; i <= lim; i++, j--)
+ for (i = 1, j = k - 1; i < j; i++, j--)
yk2 += X[i] * X[j];
yk += 2.0 * yk2;