int main (void) { mpfr_t x; mp_prec_t i, j; tests_start_mpfr (); /* checks that rounds to nearest sets the last bit to zero in case of equal distance */ mpfr_init2 (x, 59); mpfr_set_str_binary (x, "-0.10010001010111000011110010111010111110000000111101100111111E663"); if (mpfr_can_round (x, 54, GMP_RNDZ, GMP_RNDZ, 53) != 0) { printf ("Error (1) in mpfr_can_round\n"); exit (1); } mpfr_set_str_binary (x, "-Inf"); if (mpfr_can_round (x, 2000, GMP_RNDZ, GMP_RNDZ, 2000) != 0) { printf ("Error (2) in mpfr_can_round\n"); exit (1); } mpfr_set_prec (x, 64); mpfr_set_str_binary (x, "0.1011001000011110000110000110001111101011000010001110011000000000"); if (mpfr_can_round (x, 65, GMP_RNDN, GMP_RNDN, 54)) { printf ("Error (3) in mpfr_can_round\n"); exit (1); } mpfr_set_prec (x, 62); mpfr_set_str (x, "0.ff4ca619c76ba69", 16, GMP_RNDZ); for (i = 30; i < 99; i++) for (j = 30; j < 99; j++) { int r1, r2; for (r1 = 0; r1 < GMP_RND_MAX ; r1++) for (r2 = 0; r2 < GMP_RND_MAX ; r2++) mpfr_can_round (x, i, (mp_rnd_t) r1, (mp_rnd_t) r2, j); /* test for assertions */ } mpfr_clear (x); check_round_p (); tests_end_mpfr (); return 0; }
int main (void) { mpfr_t x; mpfr_prec_t i, j; tests_start_mpfr (); /* checks that rounds to nearest sets the last bit to zero in case of equal distance */ mpfr_init2 (x, 59); mpfr_set_str_binary (x, "-0.10010001010111000011110010111010111110000000111101100111111E663"); if (mpfr_can_round (x, 54, MPFR_RNDZ, MPFR_RNDZ, 53) != 0) { printf ("Error (1) in mpfr_can_round\n"); exit (1); } mpfr_set_str_binary (x, "-Inf"); if (mpfr_can_round (x, 2000, MPFR_RNDZ, MPFR_RNDZ, 2000) != 0) { printf ("Error (2) in mpfr_can_round\n"); exit (1); } mpfr_set_prec (x, 64); mpfr_set_str_binary (x, "0.1011001000011110000110000110001111101011000010001110011000000000"); if (mpfr_can_round (x, 65, MPFR_RNDN, MPFR_RNDN, 54)) { printf ("Error (3) in mpfr_can_round\n"); exit (1); } mpfr_set_prec (x, 137); mpfr_set_str_binary (x, "-0.10111001101001010110011000110100111010011101101010010100101100001110000100111111011101010110001010111100100101110111100001000010000000000E-97"); if (mpfr_can_round (x, 132, MPFR_RNDU, MPFR_RNDZ, 128) == 0) { printf ("Error (4) in mpfr_can_round\n"); exit (1); } /* in the following, we can round but cannot determine the inexact flag */ mpfr_set_prec (x, 86); mpfr_set_str_binary (x, "-0.11100100010011001111011010100111101010011000000000000000000000000000000000000000000000E-80"); if (mpfr_can_round (x, 81, MPFR_RNDU, MPFR_RNDZ, 44) == 0) { printf ("Error (5) in mpfr_can_round\n"); exit (1); } mpfr_set_prec (x, 62); mpfr_set_str (x, "0.ff4ca619c76ba69", 16, MPFR_RNDZ); for (i = 30; i < 99; i++) for (j = 30; j < 99; j++) { int r1, r2; for (r1 = 0; r1 < MPFR_RND_MAX ; r1++) for (r2 = 0; r2 < MPFR_RND_MAX ; r2++) mpfr_can_round (x, i, (mpfr_rnd_t) r1, (mpfr_rnd_t) r2, j); /* test for assertions */ } mpfr_clear (x); check_round_p (); tests_end_mpfr (); return 0; }