diff options
Diffstat (limited to 'tests/tsi_op.c')
| -rw-r--r-- | tests/tsi_op.c | 146 |
1 files changed, 146 insertions, 0 deletions
diff --git a/tests/tsi_op.c b/tests/tsi_op.c new file mode 100644 index 0000000..1422ca5 --- /dev/null +++ b/tests/tsi_op.c @@ -0,0 +1,146 @@ +/* Test file for mpfr_add_si, mpfr_sub_si, mpfr_si_sub, mpfr_mul_si, + mpfr_div_si, mpfr_si_div + +Copyright 2004, 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013 Free Software Foundation, Inc. +Contributed by the AriC and Caramel projects, INRIA. + +This file is part of the GNU MPFR Library. + +The GNU MPFR Library 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. + +The GNU MPFR Library 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 the GNU MPFR Library; see the file COPYING.LESSER. If not, see +http://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., +51 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA. */ + +#include <stdio.h> +#include <stdlib.h> + +#include "mpfr-test.h" + +#define ERROR1(s, i, z, exp) \ +{\ + printf("Error for "s" and i=%d\n", i);\ + printf("Expected %s\n", exp);\ + printf("Got "); mpfr_out_str (stdout, 16, 0, z, MPFR_RNDN);\ + putchar ('\n');\ + exit(1);\ +} + +const struct { + const char * op1; + long int op2; + const char * res_add; + const char * res_sub; + const char * res_mul; + const char * res_div; +} tab[] = { + {"10", 0x1, "11", "0F", "10", "10"}, + {"1", -1, "0", "2", "-1", "-1"}, + {"17.42", -0x17, "0.42", "2E.42", "-216.ee", "-1.02de9bd37a6f4"}, + {"-1024.0", -0x16, "-103A", "-100E", "16318", "bb.d1745d1745d0"} +}; + +static void +check_invert (void) +{ + mpfr_t x; + mpfr_init2 (x, MPFR_PREC_MIN); + + mpfr_set_ui (x, 0xC, MPFR_RNDN); + mpfr_si_sub (x, -1, x, MPFR_RNDD); /* -0001 - 1100 = - 1101 --> -1 0000 */ + if (mpfr_cmp_si (x, -0x10) ) + { + printf ("Special rounding error\n"); + exit (1); + } + mpfr_clear (x); +} + +#define TEST_FUNCTION mpfr_add_si +#define TEST_FUNCTION_NAME "mpfr_add_si" +#define INTEGER_TYPE long +#define RAND_FUNCTION(x) mpfr_random2(x, MPFR_LIMB_SIZE (x), 1, RANDS) +#define test_generic_ui test_generic_add_si +#include "tgeneric_ui.c" + +#define TEST_FUNCTION mpfr_sub_si +#define TEST_FUNCTION_NAME "mpfr_sub_si" +#define INTEGER_TYPE long +#define RAND_FUNCTION(x) mpfr_random2(x, MPFR_LIMB_SIZE (x), 1, RANDS) +#define test_generic_ui test_generic_sub_si +#include "tgeneric_ui.c" + +#define TEST_FUNCTION mpfr_mul_si +#define TEST_FUNCTION_NAME "mpfr_mul_si" +#define INTEGER_TYPE long +#define RAND_FUNCTION(x) mpfr_random2(x, MPFR_LIMB_SIZE (x), 1, RANDS) +#define test_generic_ui test_generic_mul_si +#include "tgeneric_ui.c" + +#define TEST_FUNCTION mpfr_div_si +#define TEST_FUNCTION_NAME "mpfr_div_si" +#define INTEGER_TYPE long +#define RAND_FUNCTION(x) mpfr_random2(x, MPFR_LIMB_SIZE (x), 1, RANDS) +#define test_generic_ui test_generic_div_si +#include "tgeneric_ui.c" + + +int +main (int argc, char *argv[]) +{ + mpfr_t x, z; + int y; + int i; + + tests_start_mpfr (); + mpfr_inits2 (53, x, z, (mpfr_ptr) 0); + for(i = 0 ; i < numberof (tab) ; i++) + { + mpfr_set_str (x, tab[i].op1, 16, MPFR_RNDN); + y = tab[i].op2; + mpfr_add_si (z, x, y, MPFR_RNDZ); + if (mpfr_cmp_str (z, tab[i].res_add, 16, MPFR_RNDN)) + ERROR1("add_si", i, z, tab[i].res_add); + mpfr_sub_si (z, x, y, MPFR_RNDZ); + if (mpfr_cmp_str (z, tab[i].res_sub, 16, MPFR_RNDN)) + ERROR1("sub_si", i, z, tab[i].res_sub); + mpfr_si_sub (z, y, x, MPFR_RNDZ); + mpfr_neg (z, z, MPFR_RNDZ); + if (mpfr_cmp_str (z, tab[i].res_sub, 16, MPFR_RNDN)) + ERROR1("si_sub", i, z, tab[i].res_sub); + mpfr_mul_si (z, x, y, MPFR_RNDZ); + if (mpfr_cmp_str (z, tab[i].res_mul, 16, MPFR_RNDN)) + ERROR1("mul_si", i, z, tab[i].res_mul); + mpfr_div_si (z, x, y, MPFR_RNDZ); + if (mpfr_cmp_str (z, tab[i].res_div, 16, MPFR_RNDN)) + ERROR1("div_si", i, z, tab[i].res_div); + } + mpfr_set_str1 (x, "1"); + mpfr_si_div (z, 1024, x, MPFR_RNDN); + if (mpfr_cmp_str1 (z, "1024")) + ERROR1("si_div", i, z, "1024"); + mpfr_si_div (z, -1024, x, MPFR_RNDN); + if (mpfr_cmp_str1 (z, "-1024")) + ERROR1("si_div", i, z, "-1024"); + + mpfr_clears (x, z, (mpfr_ptr) 0); + + check_invert (); + + test_generic_add_si (2, 200, 17); + test_generic_sub_si (2, 200, 17); + test_generic_mul_si (2, 200, 17); + test_generic_div_si (2, 200, 17); + + tests_end_mpfr (); + return 0; +} |