forgot two files

git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@9793 280ebfd0-de03-0410-8827-d642c229c3f4

1 files changed, 107 insertions, 0 deletions

diff --git a/tests/tfmma.c b/tests/tfmma.c
new file mode 100644
index 000000000..651d145ff
--- /dev/null
+++ b/tests/tfmma.c
@@ -0,0 +1,107 @@
+/* Test file for mpfr_fmma and mpfr_fmms.
+
+Copyright 2016 Free Software Foundation, Inc.
+Contributed by the AriC and Caramel projects, INRIA.
+
+This file is part of the GNU MPFR Library.
+
+The GNU MPFR Library is free software; you can redistribute it and/or modify
+it under the terms of the GNU Lesser General Public License as published by
+the Free Software Foundation; either version 3 of the License, or (at your
+option) any later version.
+
+The GNU MPFR Library is distributed in the hope that it will be useful, but
+WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
+or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
+License for more details.
+
+You should have received a copy of the GNU Lesser General Public License
+along with the GNU MPFR Library; see the file COPYING.LESSER.  If not, see
+http://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc.,
+51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA. */
+
+#include "mpfr-test.h"
+
+/* check both mpfr_fmma and mpfr_fmms */
+static void
+random_test (mpfr_t a, mpfr_t b, mpfr_t c, mpfr_t d, mpfr_rnd_t rnd)
+{
+  mpfr_t ref, res, ab, cd;
+  int inex_ref, inex_res;
+  mpfr_prec_t p = MPFR_PREC(a);
+
+  mpfr_init2 (res, p);
+  mpfr_init2 (ref, p);
+  mpfr_init2 (ab, mpfr_get_prec (a) +  mpfr_get_prec (b));
+  mpfr_init2 (cd, mpfr_get_prec (c) +  mpfr_get_prec (d));
+
+  /* first check fmma */
+  inex_res = mpfr_fmma (res, a, b, c, d, rnd);
+  mpfr_mul (ab, a, b, rnd);
+  mpfr_mul (cd, c, d, rnd);
+  inex_ref = mpfr_add (ref, ab, cd, rnd);
+  if (inex_res != inex_ref)
+    {
+      mpfr_printf ("mpfr_fmma failed for p=%u a=%Rf b=%Rf c=%Rf d=%Rf rnd=%d\n",
+                   p, a, b, c, d, rnd);
+      printf ("expected inex %d, got %d\n", inex_ref, inex_res);
+      exit (1);
+    }
+  if (mpfr_nan_p (res))
+    MPFR_ASSERTN (mpfr_nan_p (ref));
+  else
+    MPFR_ASSERTN (mpfr_cmp (res, ref) == 0);
+
+  /* then check fmms */
+  inex_res = mpfr_fmms (res, a, b, c, d, rnd);
+  mpfr_mul (ab, a, b, rnd);
+  mpfr_mul (cd, c, d, rnd);
+  inex_ref = mpfr_sub (ref, ab, cd, rnd);
+  if (inex_res != inex_ref)
+    {
+      mpfr_printf ("mpfr_fmms failed for p=%u a=%Rf b=%Rf c=%Rf d=%Rf rnd=%d\n",
+                   p, a, b, c, d, rnd);
+      printf ("expected inex %d, got %d\n", inex_ref, inex_res);
+      exit (1);
+    }
+  if (mpfr_nan_p (res))
+    MPFR_ASSERTN (mpfr_nan_p (ref));
+  else
+    MPFR_ASSERTN (mpfr_cmp (res, ref) == 0);
+
+  mpfr_clear (ab);
+  mpfr_clear (cd);
+  mpfr_clear (res);
+  mpfr_clear (ref);
+}
+
+static void
+random_tests (void)
+{
+  mpfr_prec_t p;
+  mpfr_rnd_t r;
+  mpfr_t a, b, c, d;
+
+  for (p = MPFR_PREC_MIN; p <= 4096; p++)
+    {
+      mpfr_inits2 (p, a, b, c, d, (mpfr_ptr) 0);
+      mpfr_urandomb (a, RANDS);
+      mpfr_urandomb (b, RANDS);
+      mpfr_urandomb (c, RANDS);
+      mpfr_urandomb (d, RANDS);
+      for (r = 0; r < 4; r++)
+        random_test (a, b, c, d, r);
+      mpfr_clears (a, b, c, d, (mpfr_ptr) 0);
+    }
+}
+
+int
+main (int argc, char *argv[])
+{
+  tests_start_mpfr ();
+
+  random_tests ();
+
+  tests_end_mpfr ();
+  return 0;
+}