| /* Test file for mpfr_set_f. |
| |
| Copyright 1999, 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 <limits.h> /* for ULONG_MAX */ |
| |
| #include "mpfr-test.h" |
| |
| int |
| main (void) |
| { |
| mpfr_t x, u; |
| mpf_t y, z; |
| mp_exp_t emax; |
| unsigned long k, pr; |
| int r, inexact; |
| |
| tests_start_mpfr (); |
| |
| mpf_init (y); |
| mpf_init (z); |
| |
| mpf_set_d (y, 0.0); |
| |
| /* check prototype of mpfr_init_set_f */ |
| mpfr_init_set_f (x, y, GMP_RNDN); |
| mpfr_set_prec (x, 100); |
| mpfr_set_f (x, y, GMP_RNDN); |
| |
| mpf_urandomb (y, RANDS, 10 * GMP_NUMB_BITS); |
| mpfr_set_f (x, y, RND_RAND ()); |
| |
| /* bug found by Jean-Pierre Merlet */ |
| mpfr_set_prec (x, 256); |
| mpf_set_prec (y, 256); |
| mpfr_init2 (u, 256); |
| mpfr_set_str (u, |
| "7.f10872b020c49ba5e353f7ced916872b020c49ba5e353f7ced916872b020c498@2", |
| 16, GMP_RNDN); |
| mpf_set_str (y, "2033033E-3", 10); /* avoid 2033.033 which is |
| locale-sensitive */ |
| mpfr_set_f (x, y, GMP_RNDN); |
| if (mpfr_cmp (x, u)) |
| { |
| printf ("mpfr_set_f failed for y=2033033E-3\n"); |
| exit (1); |
| } |
| mpf_set_str (y, "-2033033E-3", 10); /* avoid -2033.033 which is |
| locale-sensitive */ |
| mpfr_set_f (x, y, GMP_RNDN); |
| mpfr_neg (u, u, GMP_RNDN); |
| if (mpfr_cmp (x, u)) |
| { |
| printf ("mpfr_set_f failed for y=-2033033E-3\n"); |
| exit (1); |
| } |
| |
| mpf_set_prec (y, 300); |
| mpf_set_str (y, "1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111", -2); |
| mpf_mul_2exp (y, y, 600); |
| mpfr_set_prec (x, 300); |
| mpfr_set_f (x, y, GMP_RNDN); |
| if (mpfr_check (x) == 0) |
| { |
| printf ("Error in mpfr_set_f: corrupted result\n"); |
| mpfr_dump (x); |
| exit (1); |
| } |
| MPFR_ASSERTN(mpfr_cmp_ui_2exp (x, 1, 901) == 0); |
| |
| /* random values */ |
| for (k = 1; k <= 1000; k++) |
| { |
| pr = 2 + (randlimb () & 255); |
| mpf_set_prec (z, pr); |
| mpf_urandomb (z, RANDS, z->_mp_prec); |
| mpfr_set_prec (u, ((pr / BITS_PER_MP_LIMB + 1) * BITS_PER_MP_LIMB)); |
| mpfr_set_f (u, z, GMP_RNDN); |
| if (mpfr_cmp_f (u , z) != 0) |
| { |
| printf ("Error in mpfr_set_f:\n"); |
| printf ("mpf (precision=%lu)=", pr); |
| mpf_out_str (stdout, 16, 0, z); |
| printf ("\nmpfr(precision=%lu)=", |
| ((pr / BITS_PER_MP_LIMB + 1) * BITS_PER_MP_LIMB)); |
| mpfr_out_str (stdout, 16, 0, u, GMP_RNDN); |
| putchar ('\n'); |
| exit (1); |
| } |
| mpfr_set_prec (x, pr); |
| mpfr_set_f (x, z, GMP_RNDN); |
| mpfr_sub (u, u, x, GMP_RNDN); |
| mpfr_abs (u, u, GMP_RNDN); |
| if (mpfr_cmp_ui_2exp (u, 1, -pr - 1) > 0) |
| { |
| printf ("Error in mpfr_set_f: precision=%lu\n", pr); |
| printf ("mpf ="); |
| mpf_out_str (stdout, 16, 0, z); |
| printf ("\nmpfr="); |
| mpfr_out_str (stdout, 16, 0, x, GMP_RNDN); |
| putchar ('\n'); |
| exit (1); |
| } |
| } |
| |
| /* Check for +0 */ |
| mpfr_set_prec (x, 53); |
| mpf_set_prec (y, 53); |
| mpf_set_ui (y, 0); |
| for (r = 0 ; r < GMP_RND_MAX ; r++) |
| { |
| int i; |
| for (i = -1; i <= 1; i++) |
| { |
| if (i) |
| mpfr_set_si (x, i, GMP_RNDN); |
| inexact = mpfr_set_f (x, y, (mp_rnd_t) r); |
| if (!MPFR_IS_ZERO(x) || !MPFR_IS_POS(x) || inexact) |
| { |
| printf ("mpfr_set_f(x,0) failed for %s, i = %d\n", |
| mpfr_print_rnd_mode ((mp_rnd_t) r), i); |
| exit (1); |
| } |
| } |
| } |
| |
| /* coverage test */ |
| mpf_set_prec (y, 2); |
| mpfr_set_prec (x, 3 * mp_bits_per_limb); |
| mpf_set_ui (y, 1); |
| for (r = 0; r < mp_bits_per_limb; r++) |
| { |
| mpfr_urandomb (x, RANDS); /* to fill low limbs with random data */ |
| inexact = mpfr_set_f (x, y, GMP_RNDN); |
| MPFR_ASSERTN(inexact == 0 && mpfr_cmp_ui_2exp (x, 1, r) == 0); |
| mpf_mul_2exp (y, y, 1); |
| } |
| |
| mpf_set_ui (y, 1); |
| mpf_mul_2exp (y, y, ULONG_MAX); |
| mpfr_set_f (x, y, GMP_RNDN); |
| mpfr_set_ui (u, 1, GMP_RNDN); |
| mpfr_mul_2ui (u, u, ULONG_MAX, GMP_RNDN); |
| if (!mpfr_equal_p (x, u)) |
| { |
| printf ("Error: mpfr_set_f (x, y, GMP_RNDN) for y = 2^ULONG_MAX\n"); |
| exit (1); |
| } |
| |
| emax = mpfr_get_emax (); |
| |
| /* For mpf_mul_2exp, emax must fit in an unsigned long! */ |
| if (emax >= 0 && emax <= ULONG_MAX) |
| { |
| mpf_set_ui (y, 1); |
| mpf_mul_2exp (y, y, emax); |
| mpfr_set_f (x, y, GMP_RNDN); |
| mpfr_set_ui_2exp (u, 1, emax, GMP_RNDN); |
| if (!mpfr_equal_p (x, u)) |
| { |
| printf ("Error: mpfr_set_f (x, y, GMP_RNDN) for y = 2^emax\n"); |
| exit (1); |
| } |
| } |
| |
| /* For mpf_mul_2exp, emax - 1 must fit in an unsigned long! */ |
| if (emax >= 1 && emax - 1 <= ULONG_MAX) |
| { |
| mpf_set_ui (y, 1); |
| mpf_mul_2exp (y, y, emax - 1); |
| mpfr_set_f (x, y, GMP_RNDN); |
| mpfr_set_ui_2exp (u, 1, emax - 1, GMP_RNDN); |
| if (!mpfr_equal_p (x, u)) |
| { |
| printf ("Error: mpfr_set_f (x, y, GMP_RNDN) for y = 2^(emax-1)\n"); |
| exit (1); |
| } |
| } |
| |
| mpfr_clear (x); |
| mpfr_clear (u); |
| mpf_clear (y); |
| mpf_clear (z); |
| |
| tests_end_mpfr (); |
| return 0; |
| } |