From 11da511c784eca003deb90c23570f0873954e0de Mon Sep 17 00:00:00 2001 From: Duncan Wilkie Date: Sat, 18 Nov 2023 06:11:09 -0600 Subject: Initial commit. --- gmp-6.3.0/mpz/inp_str.c | 173 ++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 173 insertions(+) create mode 100644 gmp-6.3.0/mpz/inp_str.c (limited to 'gmp-6.3.0/mpz/inp_str.c') diff --git a/gmp-6.3.0/mpz/inp_str.c b/gmp-6.3.0/mpz/inp_str.c new file mode 100644 index 0000000..0756055 --- /dev/null +++ b/gmp-6.3.0/mpz/inp_str.c @@ -0,0 +1,173 @@ +/* mpz_inp_str(dest_integer, stream, base) -- Input a number in base + BASE from stdio stream STREAM and store the result in DEST_INTEGER. + + OF THE FUNCTIONS IN THIS FILE, ONLY mpz_inp_str IS FOR EXTERNAL USE, THE + REST ARE INTERNALS AND ARE ALMOST CERTAIN TO BE SUBJECT TO INCOMPATIBLE + CHANGES OR DISAPPEAR COMPLETELY IN FUTURE GNU MP RELEASES. + +Copyright 1991, 1993, 1994, 1996, 1998, 2000-2003, 2011-2013 Free Software +Foundation, Inc. + +This file is part of the GNU MP Library. + +The GNU MP Library is free software; you can redistribute it and/or modify +it under the terms of either: + + * the GNU Lesser General Public License as published by the Free + Software Foundation; either version 3 of the License, or (at your + option) any later version. + +or + + * the GNU General Public License as published by the Free Software + Foundation; either version 2 of the License, or (at your option) any + later version. + +or both in parallel, as here. + +The GNU MP Library is distributed in the hope that it will be useful, but +WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY +or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License +for more details. + +You should have received copies of the GNU General Public License and the +GNU Lesser General Public License along with the GNU MP Library. If not, +see https://www.gnu.org/licenses/. */ + +#include +#include +#include "gmp-impl.h" +#include "longlong.h" + +#define digit_value_tab __gmp_digit_value_tab + +size_t +mpz_inp_str (mpz_ptr x, FILE *stream, int base) +{ + int c; + size_t nread; + + if (stream == 0) + stream = stdin; + + nread = 0; + + /* Skip whitespace. */ + do + { + c = getc (stream); + nread++; + } + while (isspace (c)); + + return mpz_inp_str_nowhite (x, stream, base, c, nread); +} + +/* shared by mpq_inp_str */ +size_t +mpz_inp_str_nowhite (mpz_ptr x, FILE *stream, int base, int c, size_t nread) +{ + char *str; + size_t alloc_size, str_size; + int negative; + mp_size_t xsize; + const unsigned char *digit_value; + + ASSERT_ALWAYS (EOF == -1); /* FIXME: handle this by adding explicit */ + /* comparisons of c and EOF before each */ + /* read of digit_value[]. */ + + digit_value = digit_value_tab; + if (base > 36) + { + /* For bases > 36, use the collating sequence + 0123456789ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz. */ + digit_value += 208; + if (UNLIKELY (base > 62)) + return 0; /* too large base */ + } + + negative = 0; + if (c == '-') + { + negative = 1; + c = getc (stream); + nread++; + } + + if (c == EOF || digit_value[c] >= (base == 0 ? 10 : base)) + return 0; /* error if no digits */ + + /* If BASE is 0, try to find out the base by looking at the initial + characters. */ + if (base == 0) + { + base = 10; + if (c == '0') + { + base = 8; + c = getc (stream); + nread++; + if (c == 'x' || c == 'X') + { + base = 16; + c = getc (stream); + nread++; + } + else if (c == 'b' || c == 'B') + { + base = 2; + c = getc (stream); + nread++; + } + } + } + + /* Skip leading zeros. */ + while (c == '0') + { + c = getc (stream); + nread++; + } + + alloc_size = 100; + str = __GMP_ALLOCATE_FUNC_TYPE (alloc_size, char); + str_size = 0; + + while (c != EOF) + { + int dig; + dig = digit_value[c]; + if (dig >= base) + break; + if (str_size >= alloc_size) + { + size_t old_alloc_size = alloc_size; + alloc_size = alloc_size * 3 / 2; + str = __GMP_REALLOCATE_FUNC_TYPE (str, old_alloc_size, alloc_size, char); + } + str[str_size++] = dig; + c = getc (stream); + } + nread += str_size; + + ungetc (c, stream); + nread--; + + /* Make sure the string is not empty, mpn_set_str would fail. */ + if (str_size == 0) + { + SIZ (x) = 0; + } + else + { + LIMBS_PER_DIGIT_IN_BASE (xsize, str_size, base); + MPZ_NEWALLOC (x, xsize); + + /* Convert the byte array in base BASE to our bignum format. */ + xsize = mpn_set_str (PTR (x), (unsigned char *) str, str_size, base); + SIZ (x) = negative ? -xsize : xsize; + } + (*__gmp_free_func) (str, alloc_size); + return nread; +} -- cgit v1.2.3