[li2.c] fixed bug when x is near 0

[tli2.c] added new test


git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@6512 280ebfd0-de03-0410-8827-d642c229c3f4

1 files changed, 3 insertions, 0 deletions

diff --git a/li2.c b/li2.c
index 4f0fa54a4..b65d6ab25 100644
--- a/li2.c
+++ b/li2.c
@@ -276,6 +276,8 @@ mpfr_li2 (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode)
         {
           mpfr_ui_sub (u, 1, x, MPFR_RNDN);
           mpfr_log (u, u, MPFR_RNDU);
+          if (MPFR_IS_ZERO(u))
+            goto next_m;
           mpfr_neg (u, u, MPFR_RNDN);    /* u = -log(1-x) */
           expo_l = MPFR_GET_EXP (u);
           k = li2_series (s, u, MPFR_RNDU);
@@ -292,6 +294,7 @@ mpfr_li2 (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode)
           if (MPFR_CAN_ROUND (s, (mp_exp_t) m - err, yp, rnd_mode))
             break;
 
+        next_m:
           MPFR_ZIV_NEXT (loop, m);
           mpfr_set_prec (u, m);
           mpfr_set_prec (s, m);