1 /* Test mpz_set_f. 2 3 Copyright 2001 Free Software Foundation, Inc. 4 5 This file is part of the GNU MP Library test suite. 6 7 The GNU MP Library test suite is free software; you can redistribute it 8 and/or modify it under the terms of the GNU General Public License as 9 published by the Free Software Foundation; either version 3 of the License, 10 or (at your option) any later version. 11 12 The GNU MP Library test suite is distributed in the hope that it will be 13 useful, but WITHOUT ANY WARRANTY; without even the implied warranty of 14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General 15 Public License for more details. 16 17 You should have received a copy of the GNU General Public License along with 18 the GNU MP Library test suite. If not, see https://www.gnu.org/licenses/. */ 19 20 #include <stdio.h> 21 #include <stdlib.h> 22 #include "gmp-impl.h" 23 #include "tests.h" 24 25 26 void 27 check_one (mpz_srcptr z) 28 { 29 static const int shift[] = { 30 0, 1, GMP_LIMB_BITS, 2*GMP_LIMB_BITS, 5*GMP_LIMB_BITS 31 }; 32 33 int sh, shneg, neg; 34 mpf_t f; 35 mpz_t got, want; 36 37 mpf_init2 (f, mpz_sizeinbase(z,2)); 38 mpz_init (got); 39 mpz_init (want); 40 41 for (sh = 0; sh < numberof(shift); sh++) 42 { 43 for (shneg = 0; shneg <= 1; shneg++) 44 { 45 for (neg = 0; neg <= 1; neg++) 46 { 47 mpf_set_z (f, z); 48 mpz_set (want, z); 49 50 if (neg) 51 { 52 mpf_neg (f, f); 53 mpz_neg (want, want); 54 } 55 56 if (shneg) 57 { 58 mpz_tdiv_q_2exp (want, want, shift[sh]); 59 mpf_div_2exp (f, f, shift[sh]); 60 } 61 else 62 { 63 mpz_mul_2exp (want, want, shift[sh]); 64 mpf_mul_2exp (f, f, shift[sh]); 65 } 66 67 mpz_set_f (got, f); 68 MPZ_CHECK_FORMAT (got); 69 70 if (mpz_cmp (got, want) != 0) 71 { 72 printf ("wrong result\n"); 73 printf (" shift %d\n", shneg ? -shift[sh] : shift[sh]); 74 printf (" neg %d\n", neg); 75 mpf_trace (" f", f); 76 mpz_trace (" got", got); 77 mpz_trace (" want", want); 78 abort (); 79 } 80 } 81 } 82 } 83 84 mpf_clear (f); 85 mpz_clear (got); 86 mpz_clear (want); 87 } 88 89 90 void 91 check_various (void) 92 { 93 mpz_t z; 94 95 mpz_init (z); 96 97 mpz_set_ui (z, 0L); 98 check_one (z); 99 100 mpz_set_si (z, 123L); 101 check_one (z); 102 103 mpz_rrandomb (z, RANDS, 2*GMP_LIMB_BITS); 104 check_one (z); 105 106 mpz_rrandomb (z, RANDS, 5*GMP_LIMB_BITS); 107 check_one (z); 108 109 mpz_clear (z); 110 } 111 112 113 int 114 main (int argc, char *argv[]) 115 { 116 #if GMP_NAIL_BITS == 0 117 tests_start (); 118 mp_trace_base = 16; 119 120 check_various (); 121 122 tests_end (); 123 #endif 124 exit (0); 125 } 126