diff options
Diffstat (limited to 'tests/tpow_si_tmpl.c')
-rw-r--r-- | tests/tpow_si_tmpl.c | 98 |
1 files changed, 0 insertions, 98 deletions
diff --git a/tests/tpow_si_tmpl.c b/tests/tpow_si_tmpl.c deleted file mode 100644 index 142294d..0000000 --- a/tests/tpow_si_tmpl.c +++ /dev/null @@ -1,98 +0,0 @@ -/* tpow_si_tmpl.c -- test file for mpc_pow_si. - -Copyright (C) 2009, 2010, 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 <limits.h> /* for CHAR_BIT */ - -#include "templates.h" -#include "mpc-tests.h" - -static void -compare_mpc_pow (mpfr_prec_t pmax, int iter, unsigned long nbits) - /* copied from tpow_ui.c and replaced unsigned by signed */ -{ - mpfr_prec_t p; - mpc_t x, y, z, t; - long n; - int i, inex_pow, inex_pow_si; - mpc_rnd_t rnd; - - mpc_init3 (y, sizeof (unsigned long) * CHAR_BIT, MPFR_PREC_MIN); - for (p = MPFR_PREC_MIN; p <= pmax; p++) - for (i = 0; i < iter; i++) - { - mpc_init2 (x, p); - mpc_init2 (z, p); - mpc_init2 (t, p); - mpc_urandom (x, rands); - n = (signed long) gmp_urandomb_ui (rands, nbits); - mpc_set_si (y, n, MPC_RNDNN); - for (rnd = 0; rnd < 16; rnd ++) - { - inex_pow = mpc_pow (z, x, y, rnd); - inex_pow_si = mpc_pow_si (t, x, n, rnd); - if (mpc_cmp (z, t) != 0) - { - printf ("mpc_pow and mpc_pow_si differ for x="); - mpc_out_str (stdout, 10, 0, x, MPC_RNDNN); - printf (" n=%li\n", n); - printf ("mpc_pow gives "); - mpc_out_str (stdout, 10, 0, z, MPC_RNDNN); - printf ("\nmpc_pow_si gives "); - mpc_out_str (stdout, 10, 0, t, MPC_RNDNN); - printf ("\n"); - exit (1); - } - if (inex_pow != inex_pow_si) - { - printf ("mpc_pow and mpc_pow_si give different flags for x="); - mpc_out_str (stdout, 10, 0, x, MPC_RNDNN); - printf (" n=%li\n", n); - printf ("mpc_pow gives %d\n", inex_pow); - printf ("mpc_pow_si gives %d\n", inex_pow_si); - exit (1); - } - } - mpc_clear (x); - mpc_clear (z); - mpc_clear (t); - } - mpc_clear (y); -} - -#define MPC_FUNCTION_CALL \ - P[0].mpc_inex = mpc_pow_si (P[1].mpc, P[2].mpc, P[3].si, P[4].mpc_rnd) -#define MPC_FUNCTION_CALL_REUSE_OP1 \ - P[0].mpc_inex = mpc_pow_si (P[1].mpc, P[1].mpc, P[3].si, P[4].mpc_rnd) - -#include "data_check.tpl" - -int -main (void) -{ - test_start (); - - data_check_template ("pow_si.dsc", "pow_si.dat"); - - compare_mpc_pow (100, 5, 19); - - test_end (); - - return 0; -} |