| /* 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 <stdio.h> |
| #include <stdlib.h> |
| |
| /* 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 <stdint.h> |
| #endif |
| #ifdef HAVE_INTTYPES_H |
| # include <inttypes.h> |
| #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 |