From f414520d3c9d661411db5d59a21daf9790cab3f4 Mon Sep 17 00:00:00 2001
From: Siddhesh Poyarekar <siddhesh@redhat.com>
Date: Thu, 14 Feb 2013 11:39:14 +0530
Subject: [PATCH] Use __sqr instead of __mul wherever possible

---
 ChangeLog                         | 6 ++++++
 sysdeps/ieee754/dbl-64/mpatan.c   | 2 +-
 sysdeps/ieee754/dbl-64/mpsqrt.c   | 2 +-
 sysdeps/ieee754/dbl-64/sincos32.c | 4 ++--
 4 files changed, 10 insertions(+), 4 deletions(-)

diff --git a/ChangeLog b/ChangeLog
index 1121ef5..97ad9da 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,5 +1,11 @@
 2013-02-14  Siddhesh Poyarekar  <siddhesh@redhat.com>
 
+	* sysdeps/ieee754/dbl-64/mpatan.c (__mpatan): Use __sqr
+	instead of __mul.
+	* sysdeps/ieee754/dbl-64/mpsqrt.c (__mpsqrt): Likewise.
+	* sysdeps/ieee754/dbl-64/sincos32.c (ss32): Likewise.
+	(cc32): Likewise.
+
 	* sysdeps/ieee754/dbl-64/mpa.c (__sqr): New function.
 	* sysdeps/ieee754/dbl-64/mpa.h (__sqr): Declare.
 	* sysdeps/ieee754/dbl-64/mpexp.c (__mpexp): use __sqr instead
diff --git a/sysdeps/ieee754/dbl-64/mpatan.c b/sysdeps/ieee754/dbl-64/mpatan.c
index db58680..0f5a24a 100644
--- a/sysdeps/ieee754/dbl-64/mpatan.c
+++ b/sysdeps/ieee754/dbl-64/mpatan.c
@@ -66,7 +66,7 @@ __mpatan(mp_no *x, mp_no *y, int p) {
     mptwoim1.d[0] = ONE;
 
 				 /* Reduce x m times */
-    __mul(x,x,&mpsm,p);
+    __sqr(x,&mpsm,p);
     if (m==0) __cpy(x,&mps,p);
     else {
       for (i=0; i<m; i++) {
diff --git a/sysdeps/ieee754/dbl-64/mpsqrt.c b/sysdeps/ieee754/dbl-64/mpsqrt.c
index 65df9fd..941a4e9 100644
--- a/sysdeps/ieee754/dbl-64/mpsqrt.c
+++ b/sysdeps/ieee754/dbl-64/mpsqrt.c
@@ -63,7 +63,7 @@ __mpsqrt(mp_no *x, mp_no *y, int p) {
 
   m=__mpsqrt_mp[p];
   for (i=0; i<m; i++) {
-    __mul(&mpu,&mpu,&mpt1,p);
+    __sqr(&mpu,&mpt1,p);
     __mul(&mpt1,&mpz,&mpt2,p);
     __sub(&mp3halfs,&mpt2,&mpt1,p);
     __mul(&mpu,&mpt1,&mpt2,p);
diff --git a/sysdeps/ieee754/dbl-64/sincos32.c b/sysdeps/ieee754/dbl-64/sincos32.c
index 6c5ffde..5a8f1bd 100644
--- a/sysdeps/ieee754/dbl-64/sincos32.c
+++ b/sysdeps/ieee754/dbl-64/sincos32.c
@@ -67,7 +67,7 @@ ss32(mp_no *x, mp_no *y, int p) {
 #endif
   for (i=1;i<=p;i++) mpk.d[i]=0;
 
-  __mul(x,x,&x2,p);
+  __sqr(x,&x2,p);
   __cpy(&oofac27,&gor,p);
   __cpy(&gor,&sum,p);
   for (a=27.0;a>1.0;a-=2.0) {
@@ -99,7 +99,7 @@ cc32(mp_no *x, mp_no *y, int p) {
 #endif
   for (i=1;i<=p;i++) mpk.d[i]=0;
 
-  __mul(x,x,&x2,p);
+  __sqr(x,&x2,p);
   mpk.d[1]=27.0;
   __mul(&oofac27,&mpk,&gor,p);
   __cpy(&gor,&sum,p);
-- 
2.7.4