summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorvlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4>2016-09-02 07:07:31 +0000
committervlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4>2016-09-02 07:07:31 +0000
commit192e5406a88ed01405fafa5f5945f59f1f4bc687 (patch)
tree11826240a1e54657c0494249b45a2445e75db9df
parentc601294d79706976c334e4bf4799e787b69a9149 (diff)
downloadmpfr-192e5406a88ed01405fafa5f5945f59f1f4bc687.tar.gz
Removed trailing whitespace.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@10775 280ebfd0-de03-0410-8827-d642c229c3f4
-rw-r--r--src/sqrt.c10
1 files changed, 5 insertions, 5 deletions
diff --git a/src/sqrt.c b/src/sqrt.c
index ff335cd7c..332720e3f 100644
--- a/src/sqrt.c
+++ b/src/sqrt.c
@@ -519,15 +519,15 @@ mpn_rsqrtrem2 (mpfr_limb_ptr xp, mpfr_limb_srcptr ap)
mp_limb_t t1, t0, u1, u0, r0, r1, r2;
MPFR_STAT_STATIC_ASSERT (GMP_NUMB_BITS == 64);
-
+
xp[1] = mpn_rsqrtrem1 (ap[1]);
/* now we should compute x + (x/2) * (1 - a*x^2), where the upper ~40 bits of
a*x^2 should cancel with 1, thus we only need the following 40 bits */
-
+
/* xp[1] has 1+40 bits, with xp[1] <= 2^72/sqrt(ap[1]) */
umul_ppmm (t1, t0, xp[1], xp[1]);
-
+
/* now t1 has at most 18 bits, with least 16 bits being its fractional value */
u1 = ap[1] >> 48;
@@ -564,7 +564,7 @@ mpn_rsqrtrem2 (mpfr_limb_ptr xp, mpfr_limb_srcptr ap)
n := ap[3]*B^3 + ap[2]*B^2 + ap[1]*B + ap[0] = s^2 + x*B + r, with n < (s+1)^2
or equivalently x*B + r <= 2*s.
-
+
This code currently assumes GMP_NUMB_BITS = 64. */
static mp_limb_t
mpn_sqrtrem4 (mpfr_limb_ptr sp, mpfr_limb_ptr rp, mpfr_limb_srcptr ap)
@@ -582,7 +582,7 @@ mpn_sqrtrem4 (mpfr_limb_ptr sp, mpfr_limb_ptr rp, mpfr_limb_srcptr ap)
t1 = ap[3] >> 48;
t0 = (ap[3] << 16) | (ap[2] >> 48);
-
+
/* now t1:t0 is a (16+64)-bit approximation of a,
(x1*B+x0) * (t1*B+t0) = (x1*t1)*B^2 + (x1*t0+x0*t1)*B + x0*t0 */
r2 = x[1] * t1; /* r2 has 32 bits */