diff options
author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2016-09-07 01:05:34 +0000 |
---|---|---|
committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2016-09-07 01:05:34 +0000 |
commit | bcee86c8cd3990252b780ab3f6759523383fe9d7 (patch) | |
tree | 9f42ba49e2384db5c9d540846ad614e63ae85f5e /src/zeta.c | |
parent | 82f48eeabd89b6f21def52fc487afbd5f1a7a9af (diff) | |
download | mpfr-bcee86c8cd3990252b780ab3f6759523383fe9d7.tar.gz |
[src/zeta.c] Added a FIXME: An overflow in gamma(s1) does not imply
that Zeta(s) will overflow.
[tests/tzeta.c] Added a test that triggers this bug with a 32-bit ABI.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@10805 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'src/zeta.c')
-rw-r--r-- | src/zeta.c | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/src/zeta.c b/src/zeta.c index 164a4fdf8..2cbe3d19f 100644 --- a/src/zeta.c +++ b/src/zeta.c @@ -430,6 +430,8 @@ mpfr_zeta (mpfr_t z, mpfr_srcptr s, mpfr_rnd_t rnd_mode) if (MPFR_IS_INF (y)) /* Zeta(s) < 0 for -4k-2 < s < -4k, Zeta(s) > 0 for -4k < s < -4k+2 */ { + /* FIXME: An overflow in gamma(s1) does not imply that + Zeta(s) will overflow. */ mpfr_div_2ui (s1, s, 2, MPFR_RNDN); /* s/4, exact */ mpfr_frac (s1, s1, MPFR_RNDN); /* exact, -1 < s1 < 0 */ overflow = (mpfr_cmp_si_2exp (s1, -1, -1) > 0) ? -1 : 1; |