X-Git-Url: https://oss.titaniummirror.com/gitweb?a=blobdiff_plain;f=mpfr%2Ftests%2Ftget_sj.c;fp=mpfr%2Ftests%2Ftget_sj.c;h=51c293f1d21060011d6e40fbaa1c121d9094d5b8;hb=6fed43773c9b0ce596dca5686f37ac3fc0fa11c0;hp=0000000000000000000000000000000000000000;hpb=27b11d56b743098deb193d510b337ba22dc52e5c;p=msp430-gcc.git diff --git a/mpfr/tests/tget_sj.c b/mpfr/tests/tget_sj.c new file mode 100644 index 00000000..51c293f1 --- /dev/null +++ b/mpfr/tests/tget_sj.c @@ -0,0 +1,265 @@ +/* Test file for mpfr_get_sj and mpfr_get_uj. + +Copyright 2004, 2005, 2006, 2007, 2008, 2009 Free Software Foundation, Inc. +Contributed by the Arenaire and Cacao 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 2.1 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.LIB. If not, write to +the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, +MA 02110-1301, USA. */ + +#ifdef HAVE_CONFIG_H +# include "config.h" /* for a build within gmp */ +#endif + +#include +#include + +/* The ISO C99 standard specifies that in C++ implementations the + INTMAX_MAX, ... macros should only be defined if explicitly requested. */ +#if defined __cplusplus +# define __STDC_LIMIT_MACROS +# define __STDC_CONSTANT_MACROS +#endif + +#ifdef HAVE_STDINT_H +# include +#endif +#ifdef HAVE_INTTYPES_H +# include +#endif + +#include "mpfr-test.h" + +#ifndef _MPFR_H_HAVE_INTMAX_T + +int +main (void) +{ + return 0; +} + +#else + +static void +check_sj (intmax_t s, mpfr_ptr x) +{ + mpfr_t y; + int i; + + mpfr_init2 (y, MPFR_PREC (x)); + + for (i = -1; i <= 1; i++) + { + int rnd; + + mpfr_set_si_2exp (y, i, -2, GMP_RNDN); + mpfr_add (y, y, x, GMP_RNDN); + for (rnd = 0; rnd < GMP_RND_MAX; rnd++) + { + intmax_t r; + + if (rnd == GMP_RNDZ && i < 0 && s >= 0) + continue; + if (rnd == GMP_RNDZ && i > 0 && s <= 0) + continue; + if (rnd == GMP_RNDD && i < 0) + continue; + if (rnd == GMP_RNDU && i > 0) + continue; + /* rint (y) == x == s */ + r = mpfr_get_sj (y, (mp_rnd_t) rnd); + if (r != s) + { + printf ("Error in check_sj for y = "); + mpfr_out_str (stdout, 2, 0, y, GMP_RNDN); + printf (" in %s\n", mpfr_print_rnd_mode ((mp_rnd_t) rnd)); + printf ("Got %jd instead of %jd.\n", r, s); + exit (1); + } + } + } + + mpfr_clear (y); +} + +static void +check_uj (uintmax_t u, mpfr_ptr x) +{ + mpfr_t y; + int i; + + mpfr_init2 (y, MPFR_PREC (x)); + + for (i = -1; i <= 1; i++) + { + int rnd; + + mpfr_set_si_2exp (y, i, -2, GMP_RNDN); + mpfr_add (y, y, x, GMP_RNDN); + for (rnd = 0; rnd < GMP_RND_MAX; rnd++) + { + uintmax_t r; + + if (rnd == GMP_RNDZ && i < 0) + continue; + if (rnd == GMP_RNDD && i < 0) + continue; + if (rnd == GMP_RNDU && i > 0) + continue; + /* rint (y) == x == u */ + r = mpfr_get_uj (y, (mp_rnd_t) rnd); + if (r != u) + { + printf ("Error in check_uj for y = "); + mpfr_out_str (stdout, 2, 0, y, GMP_RNDN); + printf (" in %s\n", mpfr_print_rnd_mode ((mp_rnd_t) rnd)); + printf ("Got %ju instead of %ju.\n", r, u); + exit (1); + } + } + } + + mpfr_clear (y); +} + +static void +check_erange (void) +{ + mpfr_t x; + uintmax_t dl; + intmax_t d; + + /* Test for ERANGE flag + correct behaviour if overflow */ + + mpfr_init2 (x, 256); + mpfr_set_uj (x, UINTMAX_MAX, GMP_RNDN); + mpfr_clear_erangeflag (); + dl = mpfr_get_uj (x, GMP_RNDN); + if (dl != UINTMAX_MAX || mpfr_erangeflag_p ()) + { + printf ("ERROR for get_uj + ERANGE + UINTMAX_MAX (1)\n"); + exit (1); + } + mpfr_add_ui (x, x, 1, GMP_RNDN); + dl = mpfr_get_uj (x, GMP_RNDN); + if (dl != UINTMAX_MAX || !mpfr_erangeflag_p ()) + { + printf ("ERROR for get_uj + ERANGE + UINTMAX_MAX (2)\n"); + exit (1); + } + mpfr_set_sj (x, -1, GMP_RNDN); + mpfr_clear_erangeflag (); + dl = mpfr_get_uj (x, GMP_RNDN); + if (dl != 0 || !mpfr_erangeflag_p ()) + { + printf ("ERROR for get_uj + ERANGE + -1 \n"); + exit (1); + } + mpfr_set_sj (x, INTMAX_MAX, GMP_RNDN); + mpfr_clear_erangeflag (); + d = mpfr_get_sj (x, GMP_RNDN); + if (d != INTMAX_MAX || mpfr_erangeflag_p ()) + { + printf ("ERROR for get_sj + ERANGE + INTMAX_MAX (1)\n"); + exit (1); + } + mpfr_add_ui (x, x, 1, GMP_RNDN); + d = mpfr_get_sj (x, GMP_RNDN); + if (d != INTMAX_MAX || !mpfr_erangeflag_p ()) + { + printf ("ERROR for get_sj + ERANGE + INTMAX_MAX (2)\n"); + exit (1); + } + mpfr_set_sj (x, INTMAX_MIN, GMP_RNDN); + mpfr_clear_erangeflag (); + d = mpfr_get_sj (x, GMP_RNDN); + if (d != INTMAX_MIN || mpfr_erangeflag_p ()) + { + printf ("ERROR for get_sj + ERANGE + INTMAX_MIN (1)\n"); + exit (1); + } + mpfr_sub_ui (x, x, 1, GMP_RNDN); + d = mpfr_get_sj (x, GMP_RNDN); + if (d != INTMAX_MIN || !mpfr_erangeflag_p ()) + { + printf ("ERROR for get_sj + ERANGE + INTMAX_MIN (2)\n"); + exit (1); + } + + mpfr_clear (x); +} + +int +main (void) +{ + mp_prec_t prec; + mpfr_t x, y; + intmax_t s; + uintmax_t u; + + tests_start_mpfr (); + + for (u = UINTMAX_MAX, prec = 0; u != 0; u /= 2, prec++) + { } + + mpfr_init2 (x, prec + 4); + mpfr_init2 (y, prec + 4); + + mpfr_set_ui (x, 0, GMP_RNDN); + check_sj (0, x); + check_uj (0, x); + + mpfr_set_si_2exp (x, 1, prec, GMP_RNDN); + mpfr_sub_ui (x, x, 1, GMP_RNDN); /* UINTMAX_MAX */ + + mpfr_div_ui (y, x, 2, GMP_RNDZ); + mpfr_trunc (y, y); /* INTMAX_MAX */ + for (s = INTMAX_MAX; s != 0; s /= 17) + { + check_sj (s, y); + mpfr_div_ui (y, y, 17, GMP_RNDZ); + mpfr_trunc (y, y); + } + + mpfr_div_ui (y, x, 2, GMP_RNDZ); + mpfr_trunc (y, y); /* INTMAX_MAX */ + mpfr_neg (y, y, GMP_RNDN); + if (INTMAX_MIN + INTMAX_MAX != 0) + mpfr_sub_ui (y, y, 1, GMP_RNDN); /* INTMAX_MIN */ + for (s = INTMAX_MIN; s != 0; s /= 17) + { + check_sj (s, y); + mpfr_div_ui (y, y, 17, GMP_RNDZ); + mpfr_trunc (y, y); + } + + for (u = UINTMAX_MAX; u != 0; u /= 17) + { + check_uj (u, x); + mpfr_div_ui (x, x, 17, GMP_RNDZ); + mpfr_trunc (x, x); + } + + mpfr_clear (x); + mpfr_clear (y); + + check_erange (); + + tests_end_mpfr (); + return 0; +} + +#endif