1 /* Test mpf_mul, mpf_div, mpf_ui_div, and mpf_div_ui. 2 3 Copyright 1996, 2000, 2001, 2003 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 23 #include "gmp-impl.h" 24 #include "tests.h" 25 26 #ifndef SIZE 27 #define SIZE 16 28 #endif 29 30 int 31 main (int argc, char **argv) 32 { 33 mp_size_t size; 34 mp_exp_t exp; 35 int reps = 10000; 36 int i; 37 mpf_t u, v, w, x; 38 mp_size_t bprec = SIZE * GMP_LIMB_BITS; 39 mpf_t rerr, limit_rerr; 40 unsigned long ulimb, vlimb; 41 int single_flag; 42 43 tests_start (); 44 45 if (argc > 1) 46 { 47 reps = strtol (argv[1], 0, 0); 48 if (argc > 2) 49 bprec = strtol (argv[2], 0, 0); 50 } 51 52 mpf_set_default_prec (bprec); 53 54 mpf_init (rerr); 55 mpf_init (limit_rerr); 56 57 mpf_init (u); 58 mpf_init (v); 59 mpf_init (w); 60 mpf_init (x); 61 62 for (i = 0; i < reps; i++) 63 { 64 mp_size_t res_prec; 65 66 res_prec = urandom () % bprec + 1; 67 mpf_set_prec (w, res_prec); 68 mpf_set_prec (x, res_prec); 69 70 mpf_set_ui (limit_rerr, 1); 71 mpf_div_2exp (limit_rerr, limit_rerr, res_prec - 1); 72 73 single_flag = 0; 74 75 if ((urandom () & 1) != 0) 76 { 77 size = urandom () % (2 * SIZE) - SIZE; 78 exp = urandom () % SIZE; 79 mpf_random2 (u, size, exp); 80 } 81 else 82 { 83 ulimb = urandom (); 84 mpf_set_ui (u, ulimb); 85 single_flag = 1; 86 } 87 88 if ((urandom () & 1) != 0) 89 { 90 size = urandom () % (2 * SIZE) - SIZE; 91 exp = urandom () % SIZE; 92 mpf_random2 (v, size, exp); 93 } 94 else 95 { 96 vlimb = urandom (); 97 mpf_set_ui (v, vlimb); 98 single_flag = 2; 99 } 100 101 if (mpf_sgn (v) == 0) 102 continue; 103 104 mpf_div (w, u, v); 105 mpf_mul (x, w, v); 106 mpf_reldiff (rerr, u, x); 107 if (mpf_cmp (rerr, limit_rerr) > 0) 108 { 109 printf ("ERROR in mpf_mul or mpf_div after %d tests\n", i); 110 printf (" u = "); mpf_dump (u); 111 printf (" v = "); mpf_dump (v); 112 printf (" x = "); mpf_dump (x); 113 printf (" w = "); mpf_dump (w); 114 abort (); 115 } 116 117 if (single_flag == 2) 118 { 119 mpf_div_ui (x, u, vlimb); 120 mpf_reldiff (rerr, w, x); 121 if (mpf_cmp (rerr, limit_rerr) > 0) 122 { 123 printf ("ERROR in mpf_div or mpf_div_ui after %d tests\n", i); 124 printf (" u = "); mpf_dump (u); 125 printf (" v = "); mpf_dump (v); 126 printf (" x = "); mpf_dump (x); 127 printf (" w = "); mpf_dump (w); 128 abort (); 129 } 130 } 131 132 if (single_flag == 1) 133 { 134 mpf_ui_div (x, ulimb, v); 135 mpf_reldiff (rerr, w, x); 136 if (mpf_cmp (rerr, limit_rerr) > 0) 137 { 138 printf ("ERROR in mpf_div or mpf_ui_div after %d tests\n", i); 139 printf (" u = "); mpf_dump (u); 140 printf (" v = "); mpf_dump (v); 141 printf (" x = "); mpf_dump (x); 142 printf (" w = "); mpf_dump (w); 143 abort (); 144 } 145 } 146 } 147 148 mpf_clear (rerr); 149 mpf_clear (limit_rerr); 150 151 mpf_clear (u); 152 mpf_clear (v); 153 mpf_clear (w); 154 mpf_clear (x); 155 156 tests_end (); 157 exit (0); 158 } 159