diff options
author | pelissip <pelissip@280ebfd0-de03-0410-8827-d642c229c3f4> | 2004-11-19 14:38:28 +0000 |
---|---|---|
committer | pelissip <pelissip@280ebfd0-de03-0410-8827-d642c229c3f4> | 2004-11-19 14:38:28 +0000 |
commit | 527fd2f0289711eb3b1deac082d8407a8c243e59 (patch) | |
tree | 1cb67ecfb328b5c30149c5a3535b189358796416 /log.c | |
parent | b21fbe062ac1b61a48586c7ae871ef4c036ea252 (diff) | |
download | mpfr-527fd2f0289711eb3b1deac082d8407a8c243e59.tar.gz |
Improve the initial value of prec: 10% improvement (?) on average.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@3110 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'log.c')
-rw-r--r-- | log.c | 6 |
1 files changed, 3 insertions, 3 deletions
@@ -104,11 +104,11 @@ mpfr_log (mpfr_ptr r, mpfr_srcptr a, mp_rnd_t rnd_mode) q = MPFR_PREC (r); /* use initial precision about q+lg(q)+5 */ - p = q + 5 + MPFR_INT_CEIL_LOG2 (q); + p = q + 5 + 2*MPFR_INT_CEIL_LOG2 (q); /* % ~(mp_prec_t)BITS_PER_MP_LIMB ; m=q; while (m) { p++; m >>= 1; } */ - if (MPFR_LIKELY(p % BITS_PER_MP_LIMB != 0)) - p += BITS_PER_MP_LIMB - (p%BITS_PER_MP_LIMB); + /* if (MPFR_LIKELY(p % BITS_PER_MP_LIMB != 0)) + p += BITS_PER_MP_LIMB - (p%BITS_PER_MP_LIMB); */ TMP_MARK(marker); |