Index: gcc/mpfr/tests/tpow_z.c |
diff --git a/gcc/mpfr/tests/tpow_z.c b/gcc/mpfr/tests/tpow_z.c |
deleted file mode 100644 |
index c506782916a8203e2c829b919715553ad0b0b320..0000000000000000000000000000000000000000 |
--- a/gcc/mpfr/tests/tpow_z.c |
+++ /dev/null |
@@ -1,367 +0,0 @@ |
-/* Test file for mpfr_pow_z -- power function x^z with z a MPZ |
- |
-Copyright 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 <stdlib.h> |
-#include <float.h> |
-#include <math.h> |
- |
-#include "mpfr-test.h" |
- |
-#define ERROR(str) do { printf("Error for "str"\n"); exit (1); } while (0) |
- |
-static void |
-check_special (void) |
-{ |
- mpfr_t x, y; |
- mpz_t z; |
- int res; |
- |
- mpfr_init (x); |
- mpfr_init (y); |
- mpz_init (z); |
- |
- /* x^0 = 1 except for NAN */ |
- mpfr_set_ui (x, 23, GMP_RNDN); |
- mpz_set_ui (z, 0); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_cmp_ui (y, 1) != 0) |
- ERROR ("23^0"); |
- mpfr_set_nan (x); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_nan_p (y) || mpfr_cmp_si (y, 1) != 0) |
- ERROR ("NAN^0"); |
- mpfr_set_inf (x, 1); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_cmp_ui (y, 1) != 0) |
- ERROR ("INF^0"); |
- |
- /* sINF^N = INF if s==1 or n even if N > 0*/ |
- mpz_set_ui (z, 42); |
- mpfr_set_inf (x, 1); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_inf_p (y) == 0 || mpfr_sgn (y) <= 0) |
- ERROR ("INF^42"); |
- mpfr_set_inf (x, -1); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_inf_p (y) == 0 || mpfr_sgn (y) <= 0) |
- ERROR ("-INF^42"); |
- mpz_set_ui (z, 17); |
- mpfr_set_inf (x, 1); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_inf_p (y) == 0 || mpfr_sgn (y) <= 0) |
- ERROR ("INF^17"); |
- mpfr_set_inf (x, -1); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_inf_p (y) == 0 || mpfr_sgn (y) >= 0) |
- ERROR ("-INF^17"); |
- |
- mpz_set_si (z, -42); |
- mpfr_set_inf (x, 1); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_zero_p (y) == 0 || MPFR_SIGN (y) <= 0) |
- ERROR ("INF^-42"); |
- mpfr_set_inf (x, -1); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_zero_p (y) == 0 || MPFR_SIGN (y) <= 0) |
- ERROR ("-INF^-42"); |
- mpz_set_si (z, -17); |
- mpfr_set_inf (x, 1); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_zero_p (y) == 0 || MPFR_SIGN (y) <= 0) |
- ERROR ("INF^-17"); |
- mpfr_set_inf (x, -1); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_zero_p (y) == 0 || MPFR_SIGN (y) >= 0) |
- ERROR ("-INF^-17"); |
- |
- /* s0^N = +0 if s==+ or n even if N > 0*/ |
- mpz_set_ui (z, 42); |
- MPFR_SET_ZERO (x); MPFR_SET_POS (x); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_zero_p (y) == 0 || MPFR_SIGN (y) <= 0) |
- ERROR ("+0^42"); |
- MPFR_SET_NEG (x); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_zero_p (y) == 0 || MPFR_SIGN (y) <= 0) |
- ERROR ("-0^42"); |
- mpz_set_ui (z, 17); |
- MPFR_SET_POS (x); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_zero_p (y) == 0 || MPFR_SIGN (y) <= 0) |
- ERROR ("+0^17"); |
- MPFR_SET_NEG (x); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_zero_p (y) == 0 || MPFR_SIGN (y) >= 0) |
- ERROR ("-0^17"); |
- |
- mpz_set_si (z, -42); |
- MPFR_SET_ZERO (x); MPFR_SET_POS (x); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_inf_p (y) == 0 || MPFR_SIGN (y) <= 0) |
- ERROR ("+0^-42"); |
- MPFR_SET_NEG (x); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_inf_p (y) == 0 || MPFR_SIGN (y) <= 0) |
- ERROR ("-0^-42"); |
- mpz_set_si (z, -17); |
- MPFR_SET_POS (x); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_inf_p (y) == 0 || MPFR_SIGN (y) <= 0) |
- ERROR ("+0^-17"); |
- MPFR_SET_NEG (x); |
- res = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (res != 0 || mpfr_inf_p (y) == 0 || MPFR_SIGN (y) >= 0) |
- ERROR ("-0^-17"); |
- |
- mpz_clear (z); |
- mpfr_clear (y); |
- mpfr_clear (x); |
-} |
- |
-static void |
-check_integer (mp_prec_t begin, mp_prec_t end, unsigned long max) |
-{ |
- mpfr_t x, y1, y2; |
- mpz_t z; |
- unsigned long i, n; |
- mp_prec_t p; |
- int res1, res2; |
- mp_rnd_t rnd; |
- |
- mpfr_inits2 (begin, x, y1, y2, (mpfr_ptr) 0); |
- mpz_init (z); |
- for (p = begin ; p < end ; p+=4) |
- { |
- mpfr_set_prec (x, p); |
- mpfr_set_prec (y1, p); |
- mpfr_set_prec (y2, p); |
- for (i = 0 ; i < max ; i++) |
- { |
- mpz_urandomb (z, RANDS, GMP_NUMB_BITS); |
- if ((i & 1) != 0) |
- mpz_neg (z, z); |
- mpfr_urandomb (x, RANDS); |
- mpfr_mul_2ui (x, x, 1, GMP_RNDN); /* 0 <= x < 2 */ |
- rnd = RND_RAND (); |
- if (mpz_fits_slong_p (z)) |
- { |
- n = mpz_get_si (z); |
- /* printf ("New test for x=%ld\nCheck Pow_si\n", n); */ |
- res1 = mpfr_pow_si (y1, x, n, rnd); |
- /* printf ("Check pow_z\n"); */ |
- res2 = mpfr_pow_z (y2, x, z, rnd); |
- if (mpfr_cmp (y1, y2) != 0) |
- { |
- printf ("Error for p = %lu, z = %lu, rnd = %s and x = ", |
- p, n, mpfr_print_rnd_mode (rnd)); |
- mpfr_dump (x); |
- printf ("Ypowsi = "); mpfr_dump (y1); |
- printf ("Ypowz = "); mpfr_dump (y2); |
- exit (1); |
- } |
- if (res1 != res2) |
- { |
- printf ("Wrong inexact flags for p = %lu, z = %lu, rnd = %s" |
- " and x = ", p, n, mpfr_print_rnd_mode (rnd)); |
- mpfr_dump (x); |
- printf ("Ypowsi(inex = %2d) = ", res1); mpfr_dump (y1); |
- printf ("Ypowz (inex = %2d) = ", res2); mpfr_dump (y2); |
- exit (1); |
- } |
- } |
- } /* for i */ |
- } /* for p */ |
- mpfr_clears (x, y1, y2, (mpfr_ptr) 0); |
- mpz_clear (z); |
-} |
- |
-static void |
-check_regression (void) |
-{ |
- mpfr_t x, y; |
- mpz_t z; |
- int res1, res2; |
- |
- mpz_init_set_ui (z, 2026876995); |
- mpfr_init2 (x, 122); |
- mpfr_init2 (y, 122); |
- |
- mpfr_set_str_binary (x, "0.10000010010000111101001110100101101010011110011100001111000001001101000110011001001001001011001011010110110110101000111011E1"); |
- res1 = mpfr_pow_z (y, x, z, GMP_RNDU); |
- res2 = mpfr_pow_ui (x, x, 2026876995UL, GMP_RNDU); |
- if (mpfr_cmp (x, y) || res1 != res2) |
- { |
- printf ("Regression (1) tested failed (%d=?%d)\n",res1, res2); |
- printf ("pow_ui: "); mpfr_dump (x); |
- printf ("pow_z: "); mpfr_dump (y); |
- |
- exit (1); |
- } |
- |
- mpfr_clear (x); |
- mpfr_clear (y); |
- mpz_clear (z); |
-} |
- |
-/* Bug found by Kevin P. Rauch */ |
-static void |
-bug20071104 (void) |
-{ |
- mpfr_t x, y; |
- mpz_t z; |
- int inex; |
- |
- mpz_init_set_si (z, -2); |
- mpfr_inits2 (20, x, y, (mpfr_ptr) 0); |
- mpfr_set_ui (x, 0, GMP_RNDN); |
- mpfr_nextbelow (x); /* x = -2^(emin-1) */ |
- mpfr_clear_flags (); |
- inex = mpfr_pow_z (y, x, z, GMP_RNDN); |
- if (! mpfr_inf_p (y) || MPFR_SIGN (y) < 0) |
- { |
- printf ("Error in bug20071104: expected +Inf, got "); |
- mpfr_dump (y); |
- exit (1); |
- } |
- if (inex <= 0) |
- { |
- printf ("Error in bug20071104: bad ternary value (%d)\n", inex); |
- exit (1); |
- } |
- if (__gmpfr_flags != (MPFR_FLAGS_OVERFLOW | MPFR_FLAGS_INEXACT)) |
- { |
- printf ("Error in bug20071104: bad flags (%u)\n", __gmpfr_flags); |
- exit (1); |
- } |
- mpfr_clears (x, y, (mpfr_ptr) 0); |
- mpz_clear (z); |
-} |
- |
-static void |
-check_overflow (void) |
-{ |
- mpfr_t a; |
- mpz_t z; |
- unsigned long n; |
- int res; |
- |
- mpfr_init2 (a, 53); |
- |
- mpfr_set_str_binary (a, "1E10"); |
- mpz_init_set_ui (z, ULONG_MAX); |
- res = mpfr_pow_z (a, a, z, GMP_RNDN); |
- if (!MPFR_IS_INF (a) || MPFR_SIGN (a) < 0) |
- { |
- printf ("Error for (1e10)^ULONG_MAX\n"); |
- exit (1); |
- } |
- |
- /* Bug in pow_z.c up to r5109: if x = y (same mpfr_t argument), the |
- input argument is negative and not a power of two, z is positive |
- and odd, an overflow or underflow occurs, and the temporary result |
- res is positive, then the result gets a wrong sign (positive |
- instead of negative). */ |
- mpfr_set_str_binary (a, "-1.1E10"); |
- n = (ULONG_MAX ^ (ULONG_MAX >> 1)) + 1; |
- mpz_set_ui (z, n); |
- res = mpfr_pow_z (a, a, z, GMP_RNDN); |
- if (!MPFR_IS_INF (a) || MPFR_SIGN (a) > 0) |
- { |
- printf ("Error for (-1e10)^%lu, expected -Inf,\ngot ", n); |
- mpfr_dump (a); |
- exit (1); |
- } |
- |
- mpfr_clear (a); |
- mpz_clear (z); |
-} |
- |
-/* bug reported by Carl Witty (32-bit architecture) */ |
-static void |
-bug20080223 (void) |
-{ |
- mpfr_t a, exp, answer; |
- |
- mpfr_init2 (a, 53); |
- mpfr_init2 (exp, 53); |
- mpfr_init2 (answer, 53); |
- |
- mpfr_set_si (exp, -1073741824, GMP_RNDN); |
- |
- mpfr_set_str (a, "1.999999999", 10, GMP_RNDN); |
- /* a = 562949953139837/2^48 */ |
- mpfr_pow (answer, a, exp, GMP_RNDN); |
- mpfr_set_str_binary (a, "0.110110101111011001110000111111100011101000111011101E-1073741823"); |
- MPFR_ASSERTN(mpfr_cmp0 (answer, a) == 0); |
- |
- mpfr_clear (a); |
- mpfr_clear (exp); |
- mpfr_clear (answer); |
-} |
- |
-static void |
-bug20080904 (void) |
-{ |
- mpz_t exp; |
- mpfr_t a, answer; |
- mp_exp_t emin_default; |
- |
- mpz_init (exp); |
- mpfr_init2 (a, 70); |
- mpfr_init2 (answer, 70); |
- |
- emin_default = mpfr_get_emin (); |
- mpfr_set_emin (MPFR_EMIN_MIN); |
- |
- mpz_set_str (exp, "-4eb92f8c7b7bf81e", 16); |
- mpfr_set_str_binary (a, "1.110000101110100110100011111000011110111101000011111001111001010011100"); |
- |
- mpfr_pow_z (answer, a, exp, GMP_RNDN); |
- /* The correct result is near 2^(-2^62), so it underflows when |
- MPFR_EMIN_MIN > -2^62 (i.e. with 32 and 64 bits machines). */ |
- mpfr_set_str (a, "AA500C0D7A69275DBp-4632850503556296886", 16, GMP_RNDN); |
- MPFR_ASSERTN(mpfr_cmp0 (answer, a) == 0); |
- |
- mpfr_set_emin (emin_default); |
- |
- mpz_clear (exp); |
- mpfr_clear (a); |
- mpfr_clear (answer); |
-} |
- |
-int |
-main (void) |
-{ |
- tests_start_mpfr (); |
- |
- check_special (); |
- |
- check_integer (2, 163, 100); |
- check_regression (); |
- bug20071104 (); |
- bug20080223 (); |
- bug20080904 (); |
- check_overflow (); |
- |
- tests_end_mpfr (); |
- return 0; |
-} |