1 files changed, 107 insertions, 0 deletions
diff --git a/mpc/tests/tfma.c b/mpc/tests/tfma.c
new file mode 100644
index 0000000000..81bc167fe4
--- /dev/null
+++ b/mpc/tests/tfma.c
@@ -0,0 +1,107 @@
+/* tfma -- test file for mpc_fma.
+
+Copyright (C) 2011, 2012 INRIA
+
+This file is part of GNU MPC.
+
+GNU MPC 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.
+
+GNU MPC 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 this program. If not, see http://www.gnu.org/licenses/ .
+*/
+
+#include "mpc-tests.h"
+
+static void
+cmpfma (mpc_srcptr a, mpc_srcptr b, mpc_srcptr c, mpc_rnd_t rnd)
+   /* computes a*b+c with the naive and fast functions using the rounding
+      mode rnd and compares the results and return values.                                          
+      In our current test suite, all input precisions are the same, and we
+      use this precision also for the result.
+   */
+{
+   mpc_t z, t;
+   int   inex_z, inex_t;
+
+   mpc_init2 (z, MPC_MAX_PREC (a));
+   mpc_init2 (t, MPC_MAX_PREC (a));
+
+   inex_z = mpc_fma_naive (z, a, b, c, rnd);
+   inex_t = mpc_fma (t, a, b, c, rnd);
+
+   if (mpc_cmp (z, t) != 0 || inex_z != inex_t) {
+      fprintf (stderr, "fma_naive and fma differ for rnd=(%s,%s)\n",
+               mpfr_print_rnd_mode(MPC_RND_RE(rnd)),
+               mpfr_print_rnd_mode(MPC_RND_IM(rnd)));
+      MPC_OUT (a);
+      MPC_OUT (b);
+      MPC_OUT (c);
+      MPC_OUT (z);
+      MPC_OUT (t);
+      if (inex_z != inex_t) {
+         fprintf (stderr, "inex_re (z): %s\n", MPC_INEX_STR (inex_z));
+         fprintf (stderr, "inex_re (t): %s\n", MPC_INEX_STR (inex_t));
+      }
+      exit (1);
+   }
+
+  mpc_clear (z);
+  mpc_clear (t);
+}
+
+
+static void
+check_random (void)
+{
+   mpfr_prec_t prec;
+   int rnd_re, rnd_im;
+   mpc_t a, b, c;
+
+   mpc_init2 (a, 1000);
+   mpc_init2 (b, 1000);
+   mpc_init2 (c, 1000);
+
+   for (prec = 2; prec < 1000; prec = (mpfr_prec_t) (prec * 1.1 + 1)) {
+      mpc_set_prec (a, prec);
+      mpc_set_prec (b, prec);
+      mpc_set_prec (c, prec);
+
+      test_default_random (a, -1024, 1024, 128, 0);
+      test_default_random (b, -1024, 1024, 128, 0);
+      test_default_random (c, -1024, 1024, 128, 0);
+
+      for (rnd_re = 0; rnd_re < 4; rnd_re ++)
+         for (rnd_im = 0; rnd_im < 4; rnd_im ++)
+            cmpfma (a, b, c, MPC_RND (rnd_re, rnd_im));
+   }
+
+   mpc_clear (a);
+   mpc_clear (b);
+   mpc_clear (c);
+}
+
+
+int
+main (void)
+{
+  DECL_FUNC (CCCC, f, mpc_fma);
+
+  test_start ();
+
+  check_random ();
+
+  data_check (f, "fma.dat");
+  tgeneric (f, 2, 1024, 1, 256);
+
+  test_end ();
+
+  return 0;
+}