[tests/tsum.c] New tests.

git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/branches/new-sum@9339 280ebfd0-de03-0410-8827-d642c229c3f4

1 files changed, 77 insertions, 0 deletions

diff --git a/tests/tsum.c b/tests/tsum.c
index fffda07a7..108d710ff 100644
--- a/tests/tsum.c
+++ b/tests/tsum.c
@@ -471,6 +471,82 @@ check2 (void)
   mpfr_clears (sum1, sum2, s1, s2, s3, s4, (mpfr_ptr) 0);
 }
 
+/* t[i] = (2^17 - 1) * 2^(17*(i-8)) for 0 <= i <= 16.
+ * t[17] = 2^(17*9+1) * j for 4 <= j <= 7.
+ * t[18] = 2^(-17*8) * k for -3 <= k <= 3.
+ * prec = 17*9+4
+ */
+static void
+check3 (void)
+{
+  mpfr_t sum1, sum2, s1, s2, s3, t[19];
+  mpfr_ptr p[19];
+  int i, s, j, k, r, inex1, inex2;
+  int prec = 17*9+4;
+
+  mpfr_init2 (sum1, prec);
+  mpfr_init2 (sum2, prec);
+  mpfr_init2 (s1, 17*17);
+  mpfr_init2 (s2, 17*17+4);
+  mpfr_init2 (s3, 17*17+5);
+  mpfr_set_ui (s1, 0, MPFR_RNDN);
+  for (i = 0; i < 19; i++)
+    {
+      mpfr_init2 (t[i], 20);
+      p[i] = t[i];
+      if (i < 17)
+        {
+          mpfr_set_ui_2exp (t[i], 0x1ffff, 17*(i-8), MPFR_RNDN);
+          inex1 = mpfr_add (s1, s1, t[i], MPFR_RNDN);
+          MPFR_ASSERTN (inex1 == 0);
+        }
+    }
+
+  for (s = 1; s >= -1; s -= 2)
+    {
+      for (j = 4; j <= 7; j++)
+        {
+          mpfr_set_si_2exp (t[17], s*j, 17*9+1, MPFR_RNDN);
+          inex1 = mpfr_add (s2, s1, t[17], MPFR_RNDN);
+          MPFR_ASSERTN (inex1 == 0);
+          for (k = -3; k <= 3; k++)
+            {
+              mpfr_set_si_2exp (t[18], k, -17*8, MPFR_RNDN);
+              inex1 = mpfr_add (s3, s2, t[18], MPFR_RNDN);
+              MPFR_ASSERTN (inex1 == 0);
+              RND_LOOP (r)
+                {
+                  inex1 = mpfr_set (sum1, s3, (mpfr_rnd_t) r);
+                  inex2 = mpfr_sum (sum2, p, 19, (mpfr_rnd_t) r);
+                  MPFR_ASSERTN (mpfr_check (sum1));
+                  MPFR_ASSERTN (mpfr_check (sum2));
+                  if (! mpfr_equal_p (sum1, sum2) || inex1 != inex2)
+                    {
+                      printf ("Error in check3 on %s, "
+                              "s = %d, j = %d, k = %d\n",
+                              mpfr_print_rnd_mode ((mpfr_rnd_t) r),
+                              s, j, k);
+                      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 < 17; i++)
+        mpfr_neg (t[i], t[i], MPFR_RNDN);
+      mpfr_neg (s1, s1, MPFR_RNDN);
+    }
+
+  for (i = 0; i < 19; i++)
+    mpfr_clear (t[i]);
+  mpfr_clears (sum1, sum2, s1, s2, s3, (mpfr_ptr) 0);
+}
+
 /* bug reported by Joseph S. Myers on 2013-10-27
    https://sympa.inria.fr/sympa/arc/mpfr/2013-10/msg00015.html */
 static void
@@ -652,6 +728,7 @@ main (void)
   check_more_special ();
   check1 ();
   check2 ();
+  check3 ();
   bug20131027 ();
   generic_tests ();
   check_extreme ();