--- /dev/null
+/* _mpz_set_str(mp_dest, string, base) -- Convert the \0-terminated
+ string STRING in base BASE to multiple precision integer in
+ MP_DEST. Allow white space in the string. If BASE == 0 determine
+ the base in the C standard way, i.e. 0xhh...h means base 16,
+ 0oo...o means base 8, otherwise assume base 10.
+
+Copyright (C) 1991 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 the GNU General Public License as published by
+the Free Software Foundation; either version 2, or (at your option)
+any later version.
+
+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 a copy of the GNU General Public License
+along with the GNU MP Library; see the file COPYING. If not, write to
+the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. */
+
+#include "gmp.h"
+#include "gmp-impl.h"
+#include "longlong.h"
+
+enum char_type
+{
+ XX = -3,
+ SPC = -2,
+ EOF = -1
+};
+
+static signed char ascii_to_num[256] =
+{
+ EOF,XX, XX, XX, XX, XX, XX, XX, XX, SPC,SPC,XX, XX, XX, XX, XX,
+ XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX,
+ SPC,XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX,
+ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, XX, XX, XX, XX, XX, XX,
+ XX, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24,
+ 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, XX, XX, XX, XX, XX,
+ XX, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24,
+ 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, XX, XX, XX, XX, XX,
+ XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX,
+ XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX,
+ XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX,
+ XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX,
+ XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX,
+ XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX,
+ XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX,
+ XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX, XX
+};
+
+int
+#ifdef __STDC__
+_mpz_set_str (MP_INT *x, const char *str, int base)
+#else
+_mpz_set_str (x, str, base)
+ MP_INT *x;
+ const char *str;
+ int base;
+#endif
+{
+ mp_ptr xp;
+ mp_size size;
+ mp_limb big_base;
+ int indigits_per_limb;
+ int negative = 0;
+ int inp_rawchar;
+ mp_limb inp_digit;
+ mp_limb res_digit;
+ size_t str_len;
+ mp_size i;
+
+ if (str[0] == '-')
+ {
+ negative = 1;
+ str++;
+ }
+
+ if (base == 0)
+ {
+ if (str[0] == '0')
+ {
+ if (str[1] == 'x' || str[1] == 'X')
+ base = 16;
+ else
+ base = 8;
+ }
+ else
+ base = 10;
+ }
+
+ big_base = __mp_bases[base].big_base;
+ indigits_per_limb = __mp_bases[base].chars_per_limb;
+
+ str_len = strlen (str);
+
+ size = str_len / indigits_per_limb + 1;
+ if (x->alloc < size)
+ _mpz_realloc (x, size);
+ xp = x->d;
+
+ size = 0;
+
+ if ((base & (base - 1)) == 0)
+ {
+ /* The base is a power of 2. Read the input string from
+ least to most significant character/digit. */
+
+ const char *s;
+ int next_bitpos;
+ int bits_per_indigit = big_base;
+
+ /* Accept and ignore 0x or 0X before hexadecimal numbers. */
+ if (base == 16 && str[0] == '0' && (str[1] == 'x' || str[1] == 'X'))
+ {
+ str += 2;
+ str_len -= 2;
+ }
+
+ res_digit = 0;
+ next_bitpos = 0;
+
+ for (s = str + str_len - 1; s >= str; s--)
+ {
+ inp_rawchar = *s;
+ inp_digit = ascii_to_num[inp_rawchar];
+
+ if (inp_digit >= base)
+ {
+ /* Was it white space? Just ignore it. */
+ if ((char) inp_digit == (char) SPC)
+ continue;
+
+ /* We found rubbish in the string. Return -1 to indicate
+ the error. */
+ return -1;
+ }
+
+ res_digit |= inp_digit << next_bitpos;
+ next_bitpos += bits_per_indigit;
+ if (next_bitpos >= BITS_PER_MP_LIMB)
+ {
+ xp[size] = res_digit;
+ size++;
+ next_bitpos -= BITS_PER_MP_LIMB;
+ res_digit = inp_digit >> (bits_per_indigit - next_bitpos);
+ }
+ }
+
+ xp[size] = res_digit;
+ size++;
+ for (i = size - 1; i >= 0; i--)
+ {
+ if (xp[i] != 0)
+ break;
+ }
+ size = i + 1;
+ }
+ else
+ {
+ /* General case. The base is not a power of 2. */
+
+ mp_size i;
+ int j;
+ mp_limb cy;
+
+ for (;;)
+ {
+ res_digit = 0;
+ for (j = 0; j < indigits_per_limb; )
+ {
+ inp_rawchar = (unsigned char) *str++;
+ inp_digit = ascii_to_num[inp_rawchar];
+
+ /* Negative means that the character was not a proper digit. */
+ if (inp_digit >= base)
+ {
+ /* Was it white space? Just ignore it. */
+ if ((char) inp_digit == (char) SPC)
+ continue;
+
+ goto end_or_error;
+ }
+
+ res_digit = res_digit * base + inp_digit;
+
+ /* Increment the loop counter here, since it mustn't be
+ incremented when we do "continue" above. */
+ j++;
+ }
+
+ cy = res_digit;
+
+ /* Insert RES_DIGIT into the result multi prec integer. */
+ for (i = 0; i < size; i++)
+ {
+ mp_limb p1, p0;
+ umul_ppmm (p1, p0, big_base, xp[i]);
+ p0 += cy;
+ cy = p1 + (p0 < cy);
+ xp[i] = p0;
+ }
+ if (cy != 0)
+ {
+ xp[size] = cy;
+ size++;
+ }
+ }
+
+ end_or_error:
+ /* We probably have some digits in RES_DIGIT (J tells how many). */
+ if ((char) inp_digit != (char) EOF)
+ {
+ /* Error return. */
+ return -1;
+ }
+
+ /* J contains number of digits (in base BASE) remaining in
+ RES_DIGIT. */
+ if (j > 0)
+ {
+ big_base = 1;
+ do
+ {
+ big_base *= base;
+ j--;
+ }
+ while (j > 0);
+
+ cy = res_digit;
+
+ /* Insert ultimate RES_DIGIT into the result multi prec integer. */
+ for (i = 0; i < size; i++)
+ {
+ mp_limb p1, p0;
+ umul_ppmm (p1, p0, big_base, xp[i]);
+ p0 += cy;
+ cy = p1 + (p0 < cy);
+ xp[i] = p0;
+ }
+ if (cy != 0)
+ {
+ xp[size] = cy;
+ size++;
+ }
+ }
+ }
+
+ if (negative)
+ size = -size;
+ x->size = size;
+
+ return 0;
+}
projects / ghc-hetmet.git / blobdiff