1/* Test file for mpfr_set_z. 2 3Copyright 1999, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011 Free Software Foundation, Inc. 4Contributed by the Arenaire and Cacao projects, INRIA. 5 6This file is part of the GNU MPFR Library. 7 8The GNU MPFR Library is free software; you can redistribute it and/or modify 9it under the terms of the GNU Lesser General Public License as published by 10the Free Software Foundation; either version 3 of the License, or (at your 11option) any later version. 12 13The GNU MPFR Library is distributed in the hope that it will be useful, but 14WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 15or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public 16License for more details. 17 18You should have received a copy of the GNU Lesser General Public License 19along with the GNU MPFR Library; see the file COPYING.LESSER. If not, see 20http://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., 2151 Franklin St, Fifth Floor, Boston, MA 02110-1301, USA. */ 22 23#include <stdio.h> 24#include <stdlib.h> 25#include <limits.h> 26 27#include "mpfr-test.h" 28 29static void check0(void) 30{ 31 mpz_t y; 32 mpfr_t x; 33 int inexact, r; 34 35 /* Check for +0 */ 36 mpfr_init (x); 37 mpz_init (y); 38 mpz_set_si (y, 0); 39 for(r = 0; r < MPFR_RND_MAX; r++) 40 { 41 inexact = mpfr_set_z (x, y, (mpfr_rnd_t) r); 42 if (!MPFR_IS_ZERO(x) || !MPFR_IS_POS(x) || inexact) 43 { 44 printf("mpfr_set_z(x,0) failed for %s\n", 45 mpfr_print_rnd_mode ((mpfr_rnd_t) r)); 46 exit(1); 47 } 48 } 49 mpfr_clear(x); 50 mpz_clear(y); 51} 52 53/* FIXME: It'd be better to examine the actual data in an mpfr_t to see that 54 it's as expected. Comparing mpfr_set_z with mpfr_cmp or against 55 mpfr_get_si is a rather indirect test of a low level routine. */ 56 57static void 58check (long i, mpfr_rnd_t rnd) 59{ 60 mpfr_t f; 61 mpz_t z; 62 63 mpfr_init2 (f, 8 * sizeof(long)); 64 mpz_init (z); 65 mpz_set_ui (z, i); 66 mpfr_set_z (f, z, rnd); 67 if (mpfr_get_si (f, MPFR_RNDZ) != i) 68 { 69 printf ("Error in mpfr_set_z for i=%ld rnd_mode=%d\n", i, rnd); 70 exit (1); 71 } 72 mpfr_clear (f); 73 mpz_clear (z); 74} 75 76static void 77check_large (void) 78{ 79 mpz_t z; 80 mpfr_t x, y; 81 mpfr_exp_t emax, emin; 82 83 mpz_init (z); 84 mpfr_init2 (x, 160); 85 mpfr_init2 (y, 160); 86 87 mpz_set_str (z, "77031627725494291259359895954016675357279104942148788042", 10); 88 mpfr_set_z (x, z, MPFR_RNDN); 89 mpfr_set_str_binary (y, "0.1100100100001111110110101010001000100001011010001100001000110100110001001100011001100010100010111000000011011100000111001101000100101001000000100100111000001001E186"); 90 if (mpfr_cmp (x, y)) 91 { 92 printf ("Error in mpfr_set_z on large input\n"); 93 exit (1); 94 } 95 96 /* check overflow */ 97 emax = mpfr_get_emax (); 98 set_emax (2); 99 mpz_set_str (z, "7", 10); 100 mpfr_set_z (x, z, MPFR_RNDU); 101 MPFR_ASSERTN(mpfr_inf_p (x) && mpfr_sgn (x) > 0); 102 set_emax (3); 103 mpfr_set_prec (x, 2); 104 mpz_set_str (z, "7", 10); 105 mpfr_set_z (x, z, MPFR_RNDU); 106 MPFR_ASSERTN(mpfr_inf_p (x) && mpfr_sgn (x) > 0); 107 set_emax (emax); 108 109 /* check underflow */ 110 emin = mpfr_get_emin (); 111 set_emin (3); 112 mpz_set_str (z, "1", 10); 113 mpfr_set_z (x, z, MPFR_RNDZ); 114 MPFR_ASSERTN(mpfr_cmp_ui (x, 0) == 0 && MPFR_IS_POS(x)); 115 set_emin (2); 116 mpfr_set_z (x, z, MPFR_RNDN); 117 MPFR_ASSERTN(mpfr_cmp_ui (x, 0) == 0 && MPFR_IS_POS(x)); 118 set_emin (emin); 119 120 mpz_clear (z); 121 122 mpfr_clear (x); 123 mpfr_clear (y); 124} 125 126int 127main (int argc, char *argv[]) 128{ 129 long j; 130 131 tests_start_mpfr (); 132 133 check_large (); 134 check (0, MPFR_RNDN); 135 for (j = 0; j < 200000; j++) 136 check (randlimb () & LONG_MAX, RND_RAND ()); 137 check0 (); 138 139 tests_end_mpfr (); 140 141 return 0; 142} 143