diff options
Diffstat (limited to 'tests/tdiv.c')
-rw-r--r-- | tests/tdiv.c | 17 |
1 files changed, 14 insertions, 3 deletions
diff --git a/tests/tdiv.c b/tests/tdiv.c index 2b328c5f3..430064330 100644 --- a/tests/tdiv.c +++ b/tests/tdiv.c @@ -1129,8 +1129,10 @@ test_20151023 (void) /* generate a random divisor of p bits */ mpfr_urandomb (d, RANDS); - /* generate a random quotient of GMP_NUMB_BITS bits */ - mpfr_urandomb (q0, RANDS); + /* generate a random non-zero quotient of GMP_NUMB_BITS bits */ + do + mpfr_urandomb (q0, RANDS); + while (mpfr_zero_p (q0)); /* zero-pad the quotient to p bits */ inex = mpfr_prec_round (q0, p, MPFR_RNDN); MPFR_ASSERTN(inex == 0); @@ -1152,7 +1154,16 @@ test_20151023 (void) MPFR_ASSERTN(inex == 0); mpfr_nextabove (n); mpfr_div (q, n, d, MPFR_RNDN); - MPFR_ASSERTN(mpfr_cmp (q, q0) == 0); + if (! mpfr_equal_p (q, q0)) + { + printf ("Error in test_20151023 for p=%ld, rnd=RNDN\n", + (long) p); + printf ("n="); mpfr_dump (n); + printf ("d="); mpfr_dump (d); + printf ("expected q0="); mpfr_dump (q0); + printf ("got q="); mpfr_dump (q); + exit (1); + } inex = mpfr_mul (n, d, q0, MPFR_RNDN); MPFR_ASSERTN(inex == 0); |