diff options
author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2015-03-30 08:39:48 +0000 |
---|---|---|
committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2015-03-30 08:39:48 +0000 |
commit | fd8342b0954116d616cdf6cc084514c613097b88 (patch) | |
tree | 60ba1a62c91709ebf70a4c4edf2d72d1a076b880 | |
parent | abd8f071a89eec935a0a532658f264b013951f79 (diff) | |
download | mpfr-fd8342b0954116d616cdf6cc084514c613097b88.tar.gz |
[tests/tsum.c] Extend the check4 test even further.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/branches/new-sum@9351 280ebfd0-de03-0410-8827-d642c229c3f4
-rw-r--r-- | tests/tsum.c | 91 |
1 files changed, 50 insertions, 41 deletions
diff --git a/tests/tsum.c b/tests/tsum.c index 3a6cfe012..17921eed4 100644 --- a/tests/tsum.c +++ b/tests/tsum.c @@ -559,72 +559,81 @@ check3 (void) } /* Try to check MPFR_RNDN, tmd=2, rbit=0, sst=0 or 1, negative with: - * s * (2^n - 2^k) - 1 + 2^(-1) * i + 2^(-1) * j, with -1 <= i, j <= 1, + * s * (2^n - 2^k) + h + 2^(-1) * i + 2^(-1) * j, with -1 <= h, i, j <= 1, * s = -1 or 1, prec n-k. */ static void check4 (void) { - mpfr_t sum1, sum2, s1, s2, s3, t[5]; + mpfr_t sum1, sum2, s1, s2, s3, s4, t[5]; mpfr_ptr p[5]; - int i, j, k, n, s, prec, inex1, inex2; + int h, i, j, k, n, r, s, prec, inex1, inex2; - mpfr_inits2 (256, sum1, sum2, s1, s2, s3, (mpfr_ptr) 0); + mpfr_inits2 (256, sum1, sum2, s1, s2, s3, s4, (mpfr_ptr) 0); for (i = 0; i < 5; i++) { mpfr_init2 (t[i], 2); p[i] = t[i]; } - mpfr_set_si (t[0], -1, MPFR_RNDN); - for (i = -1; i <= 1; i++) + for (h = -1; h <= 1; h++) { - mpfr_set_si (t[1], i, MPFR_RNDN); - inex1 = mpfr_add (s1, t[0], t[1], MPFR_RNDN); - MPFR_ASSERTN (inex1 == 0); - for (j = -1; j <= 1; j++) + mpfr_set_si (t[0], h, MPFR_RNDN); + for (i = -1; i <= 1; i++) { - mpfr_set_si (t[2], j, MPFR_RNDN); - inex1 = mpfr_add (s2, s1, t[2], MPFR_RNDN); + mpfr_set_si (t[1], i, MPFR_RNDN); + inex1 = mpfr_add (s1, t[0], t[1], MPFR_RNDN); MPFR_ASSERTN (inex1 == 0); - for (s = -1; s <= 1; s += 2) - for (k = 1; k <= 240; k++) - { - mpfr_set_si_2exp (t[3], -s, k, MPFR_RNDN); - inex1 = mpfr_add (s3, s2, t[3], MPFR_RNDN); - MPFR_ASSERTN (inex1 == 0); - for (n = k + 2; n <= 255; n++) + for (j = -1; j <= 1; j++) + { + mpfr_set_si (t[2], j, MPFR_RNDN); + inex1 = mpfr_add (s2, s1, t[2], MPFR_RNDN); + MPFR_ASSERTN (inex1 == 0); + for (s = -1; s <= 1; s += 2) + for (k = 1; k <= 240; k++) { - prec = n - k; - mpfr_set_prec (sum1, prec); - mpfr_set_prec (sum2, prec); - mpfr_set_si_2exp (t[4], s, n, MPFR_RNDN); - inex1 = mpfr_add (sum1, s3, t[4], MPFR_RNDN); - inex2 = mpfr_sum (sum2, p, 5, MPFR_RNDN); - MPFR_ASSERTN (mpfr_check (sum1)); - MPFR_ASSERTN (mpfr_check (sum2)); - if (!(mpfr_equal_p (sum1, sum2) && - SAME_SIGN (inex1, inex2))) + mpfr_set_si_2exp (t[3], -s, k, MPFR_RNDN); + inex1 = mpfr_add (s3, s2, t[3], MPFR_RNDN); + MPFR_ASSERTN (inex1 == 0); + for (n = k + 2; n <= 255; n++) { - printf ("Error in check4 on prec = %d, " - "i = %d, j = %d, k = %d, n = %d\n", - prec, i, j, k, n); - printf ("Expected "); - mpfr_dump (sum1); - printf ("with inex = %d\n", inex1); - printf ("Got "); - mpfr_dump (sum2); - printf ("with inex = %d\n", inex2); - exit (1); + prec = n - k; + mpfr_set_prec (sum1, prec); + mpfr_set_prec (sum2, prec); + mpfr_set_si_2exp (t[4], s, n, MPFR_RNDN); + inex1 = mpfr_add (s4, s3, t[4], MPFR_RNDN); + MPFR_ASSERTN (inex1 == 0); + RND_LOOP (r) + { + inex1 = mpfr_set (sum1, s4, (mpfr_rnd_t) r); + inex2 = mpfr_sum (sum2, p, 5, (mpfr_rnd_t) r); + MPFR_ASSERTN (mpfr_check (sum1)); + MPFR_ASSERTN (mpfr_check (sum2)); + if (!(mpfr_equal_p (sum1, sum2) && + SAME_SIGN (inex1, inex2))) + { + printf ("Error in check4 on %s, prec = %d, " + "i = %d, j = %d, k = %d, n = %d\n", + mpfr_print_rnd_mode ((mpfr_rnd_t) r), + prec, i, j, k, n); + printf ("Expected "); + mpfr_dump (sum1); + printf ("with inex = %d\n", inex1); + printf ("Got "); + mpfr_dump (sum2); + printf ("with inex = %d\n", inex2); + exit (1); + } + } } } - } + } } } for (i = 0; i < 5; i++) mpfr_clear (t[i]); - mpfr_clears (sum1, sum2, s1, s2, s3, (mpfr_ptr) 0); + mpfr_clears (sum1, sum2, s1, s2, s3, s4, (mpfr_ptr) 0); } /* bug reported by Joseph S. Myers on 2013-10-27 |