/* 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;

}