--- /dev/null
+/* Test file for exceptions.
+
+Copyright 2001, 2002, 2003, 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. */
+
+#include <stdio.h>
+#include <stdlib.h>
+
+#include "mpfr-test.h"
+
+#define ERROR(s) do { printf(s"\n"); exit(1); } while(0)
+
+/* Test powerof2 */
+static void
+check_powerof2 (void)
+{
+ mpfr_t x;
+
+ mpfr_init (x);
+ mpfr_set_ui (x, 1, GMP_RNDN);
+ MPFR_ASSERTN (mpfr_powerof2_raw (x));
+ mpfr_set_ui (x, 3, GMP_RNDN);
+ MPFR_ASSERTN (!mpfr_powerof2_raw (x));
+ mpfr_clear (x);
+}
+
+/* Test default rounding mode */
+static void
+check_default_rnd (void)
+{
+ int r;
+ mp_rnd_t t;
+ for(r = 0 ; r < GMP_RND_MAX ; r++)
+ {
+ mpfr_set_default_rounding_mode ((mp_rnd_t) r);
+ t = (mpfr_get_default_rounding_mode) ();
+ if ((mp_rnd_t) r != t)
+ ERROR("ERROR in setting / getting default rounding mode (1)");
+ }
+ mpfr_set_default_rounding_mode ((mp_rnd_t) 4);
+ if (mpfr_get_default_rounding_mode() != GMP_RNDD)
+ ERROR("ERROR in setting / getting default rounding mode (2)");
+ mpfr_set_default_rounding_mode((mp_rnd_t) -1);
+ if (mpfr_get_default_rounding_mode() != GMP_RNDD)
+ ERROR("ERROR in setting / getting default rounding mode (3)");
+}
+
+static void
+check_emin_emax (void)
+{
+ mp_exp_t old_emin, old_emax;
+
+ old_emin = mpfr_get_emin ();
+ old_emax = mpfr_get_emax ();
+
+ /* Check the functions not the macros ! */
+ if ((mpfr_set_emin)(MPFR_EMIN_MIN) != 0)
+ ERROR("set_emin failed!");
+ if ((mpfr_get_emin)() != MPFR_EMIN_MIN)
+ ERROR("get_emin FAILED!");
+ if ((mpfr_set_emin)(MPFR_EMIN_MIN-1) == 0)
+ ERROR("set_emin failed! (2)");
+
+ if ((mpfr_set_emax)(MPFR_EMAX_MAX) != 0)
+ ERROR("set_emax failed!");
+ if ((mpfr_get_emax)() != MPFR_EMAX_MAX)
+ ERROR("get_emax FAILED!");
+ if ((mpfr_set_emax)(MPFR_EMAX_MAX+1) == 0)
+ ERROR("set_emax failed! (2)");
+
+ if ((mpfr_get_emin_min) () != MPFR_EMIN_MIN)
+ ERROR ("get_emin_min");
+ if ((mpfr_get_emin_max) () != MPFR_EMIN_MAX)
+ ERROR ("get_emin_max");
+ if ((mpfr_get_emax_min) () != MPFR_EMAX_MIN)
+ ERROR ("get_emax_min");
+ if ((mpfr_get_emax_max) () != MPFR_EMAX_MAX)
+ ERROR ("get_emax_max");
+
+ set_emin (old_emin);
+ set_emax (old_emax);
+}
+
+static void
+check_set_get_prec (void)
+{
+ mpfr_t x;
+
+ mpfr_init2 (x, 17);
+ if (mpfr_get_prec (x) != 17 || (mpfr_get_prec)(x) != 17)
+ ERROR ("mpfr_get_prec");
+ mpfr_clear (x);
+}
+
+static void
+mpfr_set_double_range (void)
+{
+ mpfr_set_default_prec (54);
+ if (mpfr_get_default_prec () != 54)
+ ERROR ("get_default_prec failed (1)");
+ mpfr_set_default_prec (53);
+ if ((mpfr_get_default_prec) () != 53)
+ ERROR ("get_default_prec failed (2)");
+
+ /* in double precision format, the unbiased exponent is between 0 and
+ 2047, where 0 is used for subnormal numbers, and 2047 for special
+ numbers (infinities, NaN), and the bias is 1023, thus "normal" numbers
+ have an exponent between -1022 and 1023, corresponding to numbers
+ between 2^(-1022) and previous(2^(1024)).
+ (The smallest subnormal number is 0.(0^51)1*2^(-1022)= 2^(-1074).)
+
+ The smallest normal power of two is 1.0*2^(-1022).
+ The largest normal power of two is 2^1023.
+ (We have to add one for mpfr since mantissa are between 1/2 and 1.)
+ */
+
+ set_emin (-1021);
+ set_emax (1024);
+}
+
+static void
+check_flags (void)
+{
+ mpfr_t x;
+ mp_exp_t old_emin, old_emax;
+
+ old_emin = mpfr_get_emin ();
+ old_emax = mpfr_get_emax ();
+ mpfr_init (x);
+
+ /* Check the functions not the macros ! */
+ (mpfr_clear_flags)();
+ mpfr_set_double_range ();
+
+ mpfr_set_ui (x, 1, GMP_RNDN);
+ (mpfr_clear_overflow)();
+ mpfr_mul_2exp (x, x, 1024, GMP_RNDN);
+ if (!(mpfr_overflow_p)())
+ ERROR("ERROR: No overflow detected!\n");
+
+ (mpfr_clear_underflow)();
+ mpfr_set_ui (x, 1, GMP_RNDN);
+ mpfr_div_2exp (x, x, 1025, GMP_RNDN);
+ if (!(mpfr_underflow_p)())
+ ERROR("ERROR: No underflow detected!\n");
+
+ (mpfr_clear_nanflag)();
+ MPFR_SET_NAN(x);
+ mpfr_add (x, x, x, GMP_RNDN);
+ if (!(mpfr_nanflag_p)())
+ ERROR("ERROR: No NaN flag!\n");
+
+ (mpfr_clear_inexflag)();
+ mpfr_set_ui(x, 2, GMP_RNDN);
+ mpfr_cos(x, x, GMP_RNDN);
+ if (!(mpfr_inexflag_p)())
+ ERROR("ERROR: No inexact flag!\n");
+
+ (mpfr_clear_erangeflag) ();
+ mpfr_set_ui (x, 1, GMP_RNDN);
+ mpfr_mul_2exp (x, x, 1024, GMP_RNDN);
+ mpfr_get_ui (x, GMP_RNDN);
+ if (!(mpfr_erangeflag_p)())
+ ERROR ("ERROR: No erange flag!\n");
+
+ mpfr_clear (x);
+ set_emin (old_emin);
+ set_emax (old_emax);
+}
+
+static void
+test_set_underflow (void)
+{
+ /* static to allow non-constant initialiers in r */
+ mpfr_t x, zero, min;
+ mpfr_ptr r[4];
+ int t[4] = { 1, -1, 1, -1 };
+ int i;
+ int s;
+
+ mpfr_inits (x, zero, min, (mpfr_ptr) 0);
+ mpfr_set_ui (zero, 0, GMP_RNDN);
+ mpfr_set_ui (min, 0, GMP_RNDN);
+ mpfr_nextabove (min);
+ r[0] = r[2] = min;
+ r[1] = r[3] = zero;
+ for (s = 1; s > 0; s = -1)
+ {
+ for (i = 0; i < GMP_RND_MAX ; i++)
+ {
+ int j;
+ int inex;
+
+ j = s < 0 && i > 1 ? 5 - i : i;
+ inex = mpfr_underflow (x, (mp_rnd_t) i, s);
+ if (mpfr_cmp (x, r[j]) || inex * t[j] <= 0)
+ {
+ printf ("Error in test_set_underflow, sign = %d,"
+ " rnd_mode = %s\n", s, mpfr_print_rnd_mode ((mp_rnd_t) i));
+ printf ("Got\n");
+ mpfr_out_str (stdout, 2, 0, x, GMP_RNDN);
+ printf (", inex = %d\ninstead of\n", inex);
+ mpfr_out_str (stdout, 2, 0, r[j], GMP_RNDN);
+ printf (", inex = %d\n", t[j]);
+ exit (1);
+ }
+ }
+ mpfr_neg (zero, zero, GMP_RNDN);
+ mpfr_neg (min, min, GMP_RNDN);
+ }
+ mpfr_clears (x, zero, min, (mpfr_ptr) 0);
+}
+
+static void
+test_set_overflow (void)
+{
+ /* static to allow non-constant initialiers in r */
+ mpfr_t x, inf, max;
+ mpfr_ptr r[4];
+ int t[4] = { 1, -1, 1, -1 };
+ int i;
+ int s;
+
+ mpfr_inits2 (32, x, inf, max, (mpfr_ptr) 0);
+ mpfr_set_inf (inf, 1);
+ mpfr_set_inf (max, 1);
+ mpfr_nextbelow (max);
+ r[0] = r[2] = inf;
+ r[1] = r[3] = max;
+ for (s = 1; s > 0; s = -1)
+ {
+ for (i = 0; i < GMP_RND_MAX ; i++)
+ {
+ int j;
+ int inex;
+
+ j = s < 0 && i > 1 ? 5 - i : i;
+ inex = mpfr_overflow (x, (mp_rnd_t) i, s);
+ if (mpfr_cmp (x, r[j]) || inex * t[j] <= 0)
+ {
+ printf ("Error in test_set_overflow, sign = %d,"
+ " rnd_mode = %s\n", s, mpfr_print_rnd_mode ((mp_rnd_t) i));
+ printf ("Got\n");
+ mpfr_out_str (stdout, 2, 0, x, GMP_RNDN);
+ printf (", inex = %d\ninstead of\n", inex);
+ mpfr_out_str (stdout, 2, 0, r[j], GMP_RNDN);
+ printf (", inex = %d\n", t[j]);
+ exit (1);
+ }
+ }
+ mpfr_neg (inf, inf, GMP_RNDN);
+ mpfr_neg (max, max, GMP_RNDN);
+ }
+ mpfr_clears (x, inf, max, (mpfr_ptr) 0);
+}
+
+static void
+check_set (void)
+{
+ mpfr_clear_flags ();
+
+ mpfr_set_overflow ();
+ MPFR_ASSERTN ((mpfr_overflow_p) ());
+ mpfr_set_underflow ();
+ MPFR_ASSERTN ((mpfr_underflow_p) ());
+ mpfr_set_nanflag ();
+ MPFR_ASSERTN ((mpfr_nanflag_p) ());
+ mpfr_set_inexflag ();
+ MPFR_ASSERTN ((mpfr_inexflag_p) ());
+ mpfr_set_erangeflag ();
+ MPFR_ASSERTN ((mpfr_erangeflag_p) ());
+
+ mpfr_clear_flags ();
+}
+
+int
+main (int argc, char *argv[])
+{
+ mpfr_t x, y;
+ mp_exp_t emin, emax;
+
+ tests_start_mpfr ();
+
+ test_set_underflow ();
+ test_set_overflow ();
+ check_default_rnd();
+
+ mpfr_init (x);
+ mpfr_init (y);
+
+ emin = mpfr_get_emin ();
+ emax = mpfr_get_emax ();
+ if (emin >= emax)
+ {
+ printf ("Error: emin >= emax\n");
+ exit (1);
+ }
+
+ mpfr_set_ui (x, 1, GMP_RNDN);
+ mpfr_mul_2exp (x, x, 1024, GMP_RNDN);
+ mpfr_set_double_range ();
+ mpfr_check_range (x, 0, GMP_RNDN);
+ if (!mpfr_inf_p (x) || (mpfr_sgn(x) <= 0))
+ {
+ printf ("Error: 2^1024 rounded to nearest should give +Inf\n");
+ exit (1);
+ }
+
+ set_emax (1025);
+ mpfr_set_ui (x, 1, GMP_RNDN);
+ mpfr_mul_2exp (x, x, 1024, GMP_RNDN);
+ mpfr_set_double_range ();
+ mpfr_check_range (x, 0, GMP_RNDD);
+ if (!mpfr_number_p (x))
+ {
+ printf ("Error: 2^1024 rounded down should give a normal number\n");
+ exit (1);
+ }
+
+ mpfr_set_ui (x, 1, GMP_RNDN);
+ mpfr_mul_2exp (x, x, 1023, GMP_RNDN);
+ mpfr_add (x, x, x, GMP_RNDN);
+ if (!mpfr_inf_p (x) || (mpfr_sgn(x) <= 0))
+ {
+ printf ("Error: x+x rounded to nearest for x=2^1023 should give +Inf\n");
+ printf ("emax = %ld\n", mpfr_get_emax ());
+ printf ("got "); mpfr_print_binary (x); puts ("");
+ exit (1);
+ }
+
+ mpfr_set_ui (x, 1, GMP_RNDN);
+ mpfr_mul_2exp (x, x, 1023, GMP_RNDN);
+ mpfr_add (x, x, x, GMP_RNDD);
+ if (!mpfr_number_p (x))
+ {
+ printf ("Error: x+x rounded down for x=2^1023 should give"
+ " a normal number\n");
+ exit (1);
+ }
+
+ mpfr_set_ui (x, 1, GMP_RNDN);
+ mpfr_div_2exp (x, x, 1022, GMP_RNDN);
+ mpfr_set_str_binary (y, "1.1e-1022"); /* y = 3/2*x */
+ mpfr_sub (y, y, x, GMP_RNDZ);
+ if (mpfr_cmp_ui (y, 0))
+ {
+ printf ("Error: y-x rounded to zero should give 0"
+ " for y=3/2*2^(-1022), x=2^(-1022)\n");
+ printf ("y="); mpfr_print_binary (y); puts ("");
+ exit (1);
+ }
+
+ set_emin (-1026);
+ mpfr_set_ui (x, 1, GMP_RNDN);
+ mpfr_div_2exp (x, x, 1025, GMP_RNDN);
+ mpfr_set_double_range ();
+ mpfr_check_range (x, 0, GMP_RNDN);
+ if (!MPFR_IS_ZERO (x) )
+ {
+ printf ("Error: x rounded to nearest for x=2^-1024 should give Zero\n");
+ printf ("emin = %ld\n", mpfr_get_emin ());
+ printf ("got "); mpfr_dump (x);
+ exit (1);
+ }
+
+ mpfr_clear (x);
+ mpfr_clear (y);
+
+ set_emin (emin);
+ set_emax (emax);
+
+ check_emin_emax();
+ check_flags();
+ check_set_get_prec ();
+ check_powerof2 ();
+ check_set ();
+
+ tests_end_mpfr ();
+ return 0;
+}