mirror of
https://review.haiku-os.org/buildtools
synced 2025-02-12 08:47:41 +01:00
108 lines
2.7 KiB
C
108 lines
2.7 KiB
C
/* 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;
|
|
}
|