diff options
| author | enge <enge@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2012-10-01 11:34:45 +0000 |
|---|---|---|
| committer | enge <enge@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2012-10-01 11:34:45 +0000 |
| commit | a9b6f7b27268cd2e009abc4f5bfa15cc208ac936 (patch) | |
| tree | 2252bb31e47dabf8fe4c33497faeabd0f4322a9f | |
| parent | 25728a7a94d5cb941fe576cc606e1732907b8935 (diff) | |
| download | mpc-a9b6f7b27268cd2e009abc4f5bfa15cc208ac936.tar.gz | |
log10.c, log10.dat: alternative patch using mpfr_log10
git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@1276 211d60ee-9f03-0410-a15a-8952a2c7a4e4
| -rw-r--r-- | src/log10.c | 19 | ||||
| -rw-r--r-- | tests/log10.dat | 4 |
2 files changed, 8 insertions, 15 deletions
diff --git a/src/log10.c b/src/log10.c index 6b7e687..a4ff9af 100644 --- a/src/log10.c +++ b/src/log10.c @@ -77,16 +77,6 @@ mpc_log10_aux (mpc_ptr rop, mpc_srcptr op, mpc_rnd_t rnd, int flag, int nb) ret = mpfr_set (mpc_imagref (rop), mpc_imagref (tmp), MPC_RND_IM (rnd)); break; - case 3: /* real <- log(y) */ - mpfr_log (mpc_realref (tmp), mpc_imagref (op), MPC_RND_RE (rnd)); - mpfr_div (mpc_realref (tmp), mpc_realref (tmp), log10, MPFR_RNDN); - ok = mpfr_can_round (mpc_realref (tmp), prec - 2, MPFR_RNDN, - MPFR_RNDZ, MPC_PREC_RE(rop) + - (MPC_RND_RE (rnd) == MPFR_RNDN)); - if (ok) - ret = mpfr_set (mpc_realref (rop), mpc_realref (tmp), - MPC_RND_RE (rnd)); - break; } prec += prec / 2; mpc_set_prec (tmp, prec); @@ -197,17 +187,16 @@ mpc_log10 (mpc_ptr rop, mpc_srcptr op, mpc_rnd_t rnd) { if (im_cmp > 0) { - inex_re = mpc_log10_aux (rop, op, rnd, 0, 3); + inex_re = mpfr_log10 (mpc_realref (rop), mpc_imagref (op), MPC_RND_RE (rnd)); inex_im = mpc_log10_aux (rop, op, rnd, 1, 2); /* division by 2 does not change the ternary flag */ mpfr_div_2ui (mpc_imagref (rop), mpc_imagref (rop), 1, MPFR_RNDN); } else { - ww->re[0] = *mpc_realref (op); - ww->im[0] = *mpc_imagref (op); - MPFR_CHANGE_SIGN (ww->im); - inex_re = mpc_log10_aux (rop, ww, rnd, 0, 3); + w [0] = *mpc_imagref (op); + MPFR_CHANGE_SIGN (w); + inex_re = mpfr_log10 (mpc_realref (rop), w, MPC_RND_RE (rnd)); invrnd = MPC_RND (0, INV_RND (MPC_RND_IM (rnd))); inex_im = mpc_log10_aux (rop, op, invrnd, 1, 2); /* division by 2 does not change the ternary flag */ diff --git a/tests/log10.dat b/tests/log10.dat index fe49e28..9316c79 100644 --- a/tests/log10.dat +++ b/tests/log10.dat @@ -173,3 +173,7 @@ # log10(3+I) has an exact real part (from Joseph S. Myers) 0 + 53 0.5 53 0x8f168ee8415e7p-54 2 3 2 1 N N + +# bug found by Mickael Gastineau on 29 Sep 2012 +0 + 2 +0 2 0.75 2 +0 2 1 N N +0 - 2 +0 2 -0.75 2 +0 2 -1 N N |