diff options
| author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2020-07-10 10:29:34 +0000 |
|---|---|---|
| committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2020-07-10 10:29:34 +0000 |
| commit | 3a58a02190f6edc32ea63fe2a15bafe4280afd14 (patch) | |
| tree | 03b5f37eef586d6ccc2b30891b930d67ec12ba00 | |
| parent | 809ca254ba7a3c02c940d3e85ebe73b0fe4b503e (diff) | |
| download | mpfr-3a58a02190f6edc32ea63fe2a15bafe4280afd14.tar.gz | |
[src/jyn_asympt.c] Added an assertion.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@14062 280ebfd0-de03-0410-8827-d642c229c3f4
| -rw-r--r-- | src/jyn_asympt.c | 1 |
1 files changed, 1 insertions, 0 deletions
diff --git a/src/jyn_asympt.c b/src/jyn_asympt.c index 444be1ec8..9fa603246 100644 --- a/src/jyn_asympt.c +++ b/src/jyn_asympt.c @@ -54,6 +54,7 @@ FUNCTION (mpfr_ptr res, long n, mpfr_srcptr z, mpfr_rnd_t r) /* The terms of the asymptotic expansion grow like mu^(2k)/(8z)^(2k), where mu = 4n^2, thus we need mu < 8|z| so that it converges, i.e., n^2/2 < |z| */ + MPFR_ASSERTD (n >= 0); mpfr_set_ui (c, n, MPFR_RNDU); mpfr_mul_ui (c, c, n, MPFR_RNDU); mpfr_div_2ui (c, c, 1, MPFR_RNDU); |