Add support for MPFR_FAST_COMPUTE_IS_SMALL_INPUT

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

1 files changed, 3 insertions, 0 deletions

diff --git a/expm1.c b/expm1.c
index 841d369f9..f2cef1036 100644
--- a/expm1.c
+++ b/expm1.c
@@ -60,6 +60,9 @@ mpfr_expm1 (mpfr_ptr y, mpfr_srcptr x , mp_rnd_t rnd_mode)
 	}
     }
 
+  /* exp(x)-1 = x +x^2/2 + ... so the error is < 2^(2*EXP(x)-1) */
+  MPFR_FAST_COMPUTE_IF_SMALL_INPUT (y, x, -MPFR_GET_EXP (x)+1,1,rnd_mode,);
+
   MPFR_SAVE_EXPO_MARK (expo);
   /* General case */
   {