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. 6 7 The GNU MP Library is free software; you can redistribute it and/or modify 8 it under the terms of the GNU Lesser General Public License as published by 9 the Free Software Foundation; either version 3 of the License, or (at your 10 option) any later version. 11 12 The GNU MP Library is distributed in the hope that it will be useful, but 13 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 14 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public 15 License for more details. 16 17 You should have received a copy of the GNU Lesser General Public License 18 along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */ 19 20 #include <stdio.h> 21 #include <stdlib.h> 22 23 #include "gmp.h" 24 #include "gmp-impl.h" 25 #include "tests.h" 26 27 #ifndef SIZE 28 #define SIZE 16 29 #endif 30 31 int 32 main (int argc, char **argv) 33 { 34 mp_size_t size; 35 mp_exp_t exp; 36 int reps = 10000; 37 int i; 38 mpf_t u, v, w, x; 39 mp_size_t bprec = SIZE * GMP_LIMB_BITS; 40 mpf_t rerr, limit_rerr; 41 unsigned long ulimb, vlimb; 42 int single_flag; 43 44 tests_start (); 45 46 if (argc > 1) 47 { 48 reps = strtol (argv[1], 0, 0); 49 if (argc > 2) 50 bprec = strtol (argv[2], 0, 0); 51 } 52 53 mpf_set_default_prec (bprec); 54 55 mpf_init (rerr); 56 mpf_init (limit_rerr); 57 58 mpf_init (u); 59 mpf_init (v); 60 mpf_init (w); 61 mpf_init (x); 62 63 for (i = 0; i < reps; i++) 64 { 65 mp_size_t res_prec; 66 67 res_prec = urandom () % bprec + 1; 68 mpf_set_prec (w, res_prec); 69 mpf_set_prec (x, res_prec); 70 71 mpf_set_ui (limit_rerr, 1); 72 mpf_div_2exp (limit_rerr, limit_rerr, res_prec - 1); 73 74 single_flag = 0; 75 76 if ((urandom () & 1) != 0) 77 { 78 size = urandom () % (2 * SIZE) - SIZE; 79 exp = urandom () % SIZE; 80 mpf_random2 (u, size, exp); 81 } 82 else 83 { 84 ulimb = urandom (); 85 mpf_set_ui (u, ulimb); 86 single_flag = 1; 87 } 88 89 if ((urandom () & 1) != 0) 90 { 91 size = urandom () % (2 * SIZE) - SIZE; 92 exp = urandom () % SIZE; 93 mpf_random2 (v, size, exp); 94 } 95 else 96 { 97 vlimb = urandom (); 98 mpf_set_ui (v, vlimb); 99 single_flag = 2; 100 } 101 102 if (mpf_sgn (v) == 0) 103 continue; 104 105 mpf_div (w, u, v); 106 mpf_mul (x, w, v); 107 mpf_reldiff (rerr, u, x); 108 if (mpf_cmp (rerr, limit_rerr) > 0) 109 { 110 printf ("ERROR in mpf_mul or mpf_div after %d tests\n", i); 111 printf (" u = "); mpf_dump (u); 112 printf (" v = "); mpf_dump (v); 113 printf (" x = "); mpf_dump (x); 114 printf (" w = "); mpf_dump (w); 115 abort (); 116 } 117 118 if (single_flag == 2) 119 { 120 mpf_div_ui (x, u, vlimb); 121 mpf_reldiff (rerr, w, x); 122 if (mpf_cmp (rerr, limit_rerr) > 0) 123 { 124 printf ("ERROR in mpf_div or mpf_div_ui after %d tests\n", i); 125 printf (" u = "); mpf_dump (u); 126 printf (" v = "); mpf_dump (v); 127 printf (" x = "); mpf_dump (x); 128 printf (" w = "); mpf_dump (w); 129 abort (); 130 } 131 } 132 133 if (single_flag == 1) 134 { 135 mpf_ui_div (x, ulimb, v); 136 mpf_reldiff (rerr, w, x); 137 if (mpf_cmp (rerr, limit_rerr) > 0) 138 { 139 printf ("ERROR in mpf_div or mpf_ui_div after %d tests\n", i); 140 printf (" u = "); mpf_dump (u); 141 printf (" v = "); mpf_dump (v); 142 printf (" x = "); mpf_dump (x); 143 printf (" w = "); mpf_dump (w); 144 abort (); 145 } 146 } 147 } 148 149 mpf_clear (rerr); 150 mpf_clear (limit_rerr); 151 152 mpf_clear (u); 153 mpf_clear (v); 154 mpf_clear (w); 155 mpf_clear (x); 156 157 tests_end (); 158 exit (0); 159 } 160