From 7ccad68acf3805b35721e861fb9bc5a9e40effc8 Mon Sep 17 00:00:00 2001 From: Pedro Alvarez Date: Tue, 31 May 2016 10:43:16 +0100 Subject: Import http://www.multiprecision.org/mpc/download/mpc-1.0.3.tar.gz --- mpc/tests/tadd_fr.c | 72 +++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 72 insertions(+) create mode 100644 mpc/tests/tadd_fr.c (limited to 'mpc/tests/tadd_fr.c') diff --git a/mpc/tests/tadd_fr.c b/mpc/tests/tadd_fr.c new file mode 100644 index 0000000000..1f3b76a9a6 --- /dev/null +++ b/mpc/tests/tadd_fr.c @@ -0,0 +1,72 @@ +/* tadd_fr -- test file for mpc_add_fr. + +Copyright (C) 2008, 2010, 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 +#include "mpc-tests.h" + +static void +check_ternary_value (mpfr_prec_t prec_max, mpfr_prec_t step) +{ + mpfr_prec_t prec; + mpc_t z; + mpfr_t f; + + mpc_init2 (z, 2); + mpfr_init (f); + + for (prec = 2; prec < prec_max; prec += step) + { + mpc_set_prec (z, prec); + mpfr_set_prec (f, prec); + + mpc_set_ui (z, 1, MPC_RNDNN); + mpfr_set_ui (f, 1, GMP_RNDN); + if (mpc_add_fr (z, z, f, MPC_RNDNZ)) + { + printf ("Error in mpc_add_fr: 1+1 should be exact\n"); + exit (1); + } + + mpc_set_ui (z, 1, MPC_RNDNN); + mpc_mul_2ui (z, z, (unsigned long int) prec, MPC_RNDNN); + if (mpc_add_fr (z, z, f, MPC_RNDNN) == 0) + { + fprintf (stderr, "Error in mpc_add_fr: 2^prec+1 cannot be exact\n"); + exit (1); + } + } + mpc_clear (z); + mpfr_clear (f); +} + +int +main (void) +{ + DECL_FUNC (CCF, f, mpc_add_fr); + test_start (); + + check_ternary_value (1024, 1); + + data_check (f, "add_fr.dat"); + tgeneric (f, 2, 1024, 7, 10); + + test_end (); + return 0; +} -- cgit v1.2.1