summaryrefslogtreecommitdiff
path: root/src/grandom.c
diff options
context:
space:
mode:
authorzimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>2018-03-29 13:17:23 +0000
committerzimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>2018-03-29 13:17:23 +0000
commit138bb31b5bc29f835dec02f5c092691806da5291 (patch)
treea856429513e7f5ecfbcfd72ae940c2c8e39250c7 /src/grandom.c
parentbae8a46ce2ea12035a35c51f0b44142ed9e724b2 (diff)
downloadmpfr-138bb31b5bc29f835dec02f5c092691806da5291.tar.gz
[tests/tgrandom.c] improve coverage
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@12514 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'src/grandom.c')
-rw-r--r--src/grandom.c3
1 files changed, 2 insertions, 1 deletions
diff --git a/src/grandom.c b/src/grandom.c
index 752312e0d..7a30136cb 100644
--- a/src/grandom.c
+++ b/src/grandom.c
@@ -104,7 +104,8 @@ mpfr_grandom (mpfr_ptr rop1, mpfr_ptr rop2, gmp_randstate_t rstate,
0 < x', y' <= 2^p', and we conclude by induction. */
if (mpz_sizeinbase (s, 2) <= 2 * tprec)
goto yeepee;
- /* Extend by 32 bits */
+ /* Extend by 32 bits: for tprec=12, the probability we get here
+ is 8191/13180825, i.e., about 0.000621 */
mpz_mul_2exp (xp, xp, 32);
mpz_mul_2exp (yp, yp, 32);
mpz_urandomb (x, rstate, 32);