1 /* Test mpf_get_str and mpf_set_str. 2 3 Copyright 1996, 2000, 2001, 2008 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 #include <string.h> /* for strlen */ 23 24 #include "gmp.h" 25 #include "gmp-impl.h" 26 #include "tests.h" 27 28 #ifndef SIZE 29 #define SIZE 10 30 #endif 31 32 #ifndef EXPO 33 #define EXPO 200 34 #endif 35 36 int 37 main (int argc, char **argv) 38 { 39 mpf_t x, y; 40 int reps = 20000; 41 int i; 42 mp_size_t bprec = 100; 43 mpf_t d, rerr, max_rerr, limit_rerr; 44 char *str; 45 mp_exp_t bexp; 46 long size, exp; 47 int base; 48 char buf[SIZE * GMP_LIMB_BITS + 5]; 49 50 tests_start (); 51 52 if (argc > 1) 53 { 54 reps = strtol (argv[1], 0, 0); 55 if (argc > 2) 56 bprec = strtol (argv[2], 0, 0); 57 } 58 59 mpf_set_default_prec (bprec); 60 61 mpf_init_set_ui (limit_rerr, 1); 62 mpf_div_2exp (limit_rerr, limit_rerr, bprec); 63 #if VERBOSE 64 mpf_dump (limit_rerr); 65 #endif 66 mpf_init (rerr); 67 mpf_init_set_ui (max_rerr, 0); 68 69 mpf_init (x); 70 mpf_init (y); 71 mpf_init (d); 72 73 /* First test some specific values. */ 74 75 mpf_set_str (y, "1.23456e1000", 0); 76 77 mpf_set_str (x, "1.23456e1000", 10); 78 if (mpf_cmp (x, y) != 0) 79 abort (); 80 mpf_set_str (x, "1.23456e+1000", 0); 81 if (mpf_cmp (x, y) != 0) 82 abort (); 83 mpf_set_str (x, "1.23456e+1000", 10); 84 if (mpf_cmp (x, y) != 0) 85 abort (); 86 87 /* Now test random values. */ 88 89 for (i = 0; i < reps; i++) 90 { 91 if (i == 0) 92 { 93 /* exercise the special case in get_str for for x==0 */ 94 mpf_set_ui (x, 0L); 95 base = 10; 96 } 97 else 98 { 99 size = urandom () % (2 * SIZE) - SIZE; 100 exp = urandom () % EXPO; 101 mpf_random2 (x, size, exp); 102 base = urandom () % 61 + 2; 103 } 104 105 str = mpf_get_str (0, &bexp, base, 0, x); 106 107 if (str[0] == '-') 108 sprintf (buf, "-0.%s@%ld", str + 1, bexp); 109 else 110 sprintf (buf, "0.%s@%ld", str, bexp); 111 112 mpf_set_str_or_abort (y, buf, -base); 113 (*__gmp_free_func) (str, strlen (str) + 1); 114 115 mpf_reldiff (rerr, x, y); 116 if (mpf_cmp (rerr, max_rerr) > 0) 117 { 118 mpf_set (max_rerr, rerr); 119 #if VERBOSE 120 mpf_dump (max_rerr); 121 #endif 122 if (mpf_cmp (rerr, limit_rerr) > 0) 123 { 124 printf ("ERROR after %d tests\n", i); 125 printf ("base = %d\n", base); 126 printf (" x = "); mpf_dump (x); 127 printf (" y = "); mpf_dump (y); 128 abort (); 129 } 130 } 131 } 132 133 mpf_clear (limit_rerr); 134 mpf_clear (rerr); 135 mpf_clear (max_rerr); 136 137 mpf_clear (x); 138 mpf_clear (y); 139 mpf_clear (d); 140 141 tests_end (); 142 exit (0); 143 } 144