diff options
author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2012-05-03 14:44:57 +0000 |
---|---|---|
committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2012-05-03 14:44:57 +0000 |
commit | 0220288c9f63d5f22caf4de140057374441afa8d (patch) | |
tree | 70f23a5bcd8dcdf2cd327a241c82165090a6224f | |
parent | 11a523d0f04a98ca2fce87ca98bb594fdd37a624 (diff) | |
download | mpfr-0220288c9f63d5f22caf4de140057374441afa8d.tar.gz |
[src/gamma.c] Reverted changeset r8180, which seems incorrect.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/branches/3.1@8184 280ebfd0-de03-0410-8827-d642c229c3f4
-rw-r--r-- | src/gamma.c | 12 |
1 files changed, 0 insertions, 12 deletions
diff --git a/src/gamma.c b/src/gamma.c index ed2c90f91..a346f8b1b 100644 --- a/src/gamma.c +++ b/src/gamma.c @@ -377,18 +377,6 @@ mpfr_gamma (mpfr_ptr gamma, mpfr_srcptr x, mpfr_rnd_t rnd_mode) mpfr_mul (GammaTrial, tmp2, xp, MPFR_RNDN); /* Pi*(2-x), error (1+u)^2 */ err_g = MPFR_GET_EXP(GammaTrial); mpfr_sin (GammaTrial, GammaTrial, MPFR_RNDN); /* sin(Pi*(2-x)) */ - /* if tmp is +Inf, there is an underflow, since the - Pi*(x-1)/sin(Pi*(2-x)) term is larger than 1 in absolute value. - The sign is that of -sin(Pi*(2-x)). */ - if (mpfr_inf_p (tmp)) - { - int sgn = mpfr_sgn (GammaTrial); - MPFR_ZIV_FREE (loop); - MPFR_GROUP_CLEAR (group); - mpz_clear (fact); - MPFR_SAVE_EXPO_FREE (expo); - return mpfr_underflow (gamma, (rnd_mode == MPFR_RNDN) ? MPFR_RNDZ : rnd_mode, -sgn); - } err_g = err_g + 1 - MPFR_GET_EXP(GammaTrial); /* let g0 the true value of Pi*(2-x), g the computed value. We have g = g0 + h with |h| <= |(1+u^2)-1|*g. |