fixed bug in erfc (forgot to compare result of cmp_si)

fixed infinite loop in eint (test for use of asympt was not tight enough)


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

1 files changed, 5 insertions, 1 deletions

diff --git a/eint.c b/eint.c
index 9153350a6..227d196f8 100644
--- a/eint.c
+++ b/eint.c
@@ -273,7 +273,11 @@ mpfr_eint (mpfr_ptr y, mpfr_srcptr x, mp_rnd_t rnd)
   MPFR_ZIV_INIT (loop, prec);            /* Initialize the ZivLoop controler */
   for (;;)                               /* Infinite loop */
     {
-      if (mpfr_cmp_d (x, (double) prec * LOG2) > 0)
+      /* We need that the smallest value of k!/x^k is smaller than 2^(-p).
+	 The minimum is obtained for x=k, and it is smaller than e*sqrt(x)/e^x
+	 for x>=1. */
+      if (MPFR_GET_EXP (x) > 0 && mpfr_cmp_d (x, ((double) prec +
+			    0.5 * (double) MPFR_GET_EXP (x)) * LOG2 + 1.0) > 0)
 	err = mpfr_eint_asympt (tmp, x);
       else
 	{