diff options
author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2011-07-28 18:50:21 +0000 |
---|---|---|
committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2011-07-28 18:50:21 +0000 |
commit | 8d6b238ec8ab4b8521a8b1cc5ce3a310b6d05fef (patch) | |
tree | fdb1c7cf9dbd2c6b10675312f968c86a441f1726 /src | |
parent | 95736ac3703a0f7695e870045e0433f48907814f (diff) | |
download | mpfr-8d6b238ec8ab4b8521a8b1cc5ce3a310b6d05fef.tar.gz |
[src/mulders.c] MPFR coding conventions for r7752.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@7759 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'src')
-rw-r--r-- | src/mulders.c | 17 |
1 files changed, 10 insertions, 7 deletions
diff --git a/src/mulders.c b/src/mulders.c index fa9528234..8225448cd 100644 --- a/src/mulders.c +++ b/src/mulders.c @@ -67,7 +67,8 @@ mpfr_mulhigh_n_basecase (mpfr_limb_ptr rp, mpfr_limb_srcptr up, /* Put in rp[0..n] the n+1 low limbs of {up, n} * {vp, n}. Assume 2n limbs are allocated at rp. */ static void -mpfr_mullow_n_basecase (mp_ptr rp, mp_srcptr up, mp_srcptr vp, mp_size_t n) +mpfr_mullow_n_basecase (mpfr_limb_ptr rp, mpfr_limb_srcptr up, + mpfr_limb_srcptr vp, mp_size_t n) { mp_size_t i; @@ -116,7 +117,8 @@ mpfr_mulhigh_n (mpfr_limb_ptr rp, mpfr_limb_srcptr np, mpfr_limb_srcptr mp, /* Put in rp[0..n] the n+1 low limbs of {np, n} * {mp, n}. Assume 2n limbs are allocated at rp. */ void -mpfr_mullow_n (mp_ptr rp, mp_srcptr np, mp_srcptr mp, mp_size_t n) +mpfr_mullow_n (mpfr_limb_ptr rp, mpfr_limb_srcptr np, mpfr_limb_srcptr mp, + mp_size_t n) { mp_size_t k; @@ -239,10 +241,11 @@ mpfr_divhigh_n (mpfr_limb_ptr qp, mpfr_limb_ptr np, mpfr_limb_ptr dp, } #else /* below is the FoldDiv(K) algorithm from [1] */ mp_limb_t -mpfr_divhigh_n (mp_ptr qp, mp_ptr np, mp_ptr dp, mp_size_t n) +mpfr_divhigh_n (mpfr_limb_ptr qp, mpfr_limb_ptr np, mpfr_limb_ptr dp, + mp_size_t n) { mp_size_t k, r; - mp_ptr ip, tp, up; + mpfr_limb_ptr ip, tp, up; mp_limb_t qh = 0, cy, cc; int count; MPFR_TMP_DECL(marker); @@ -254,9 +257,9 @@ mpfr_divhigh_n (mp_ptr qp, mp_ptr np, mp_ptr dp, mp_size_t n) k = (n - 1) / K + 1; /* ceil(n/K) */ MPFR_TMP_MARK (marker); - ip = MPFR_TMP_ALLOC ((k + 1) * sizeof (mp_limb_t)); - tp = MPFR_TMP_ALLOC ((n + k) * sizeof (mp_limb_t)); - up = MPFR_TMP_ALLOC (2 * (k + 1) * sizeof (mp_limb_t)); + ip = MPFR_TMP_LIMBS_ALLOC (k + 1); + tp = MPFR_TMP_LIMBS_ALLOC (n + k); + up = MPFR_TMP_LIMBS_ALLOC (2 * (k + 1)); mpn_invert (ip, dp + n - (k + 1), k + 1, NULL); /* takes about 13% for n=1000 */ /* {ip, k+1} = floor((B^(2k+2)-1)/D - B^(k+1) where D = {dp+n-(k+1),k+1} */ for (r = n, cc = 0UL; r > 0;) |