set_str.c
上传用户:qaz666999
上传日期:2022-08-06
资源大小:2570k
文件大小:3k
源码类别:

数学计算

开发平台:

Unix_Linux

  1. /* mpz_set_str(mp_dest, string, base) -- Convert the -terminated
  2.    string STRING in base BASE to multiple precision integer in
  3.    MP_DEST.  Allow white space in the string.  If BASE == 0 determine
  4.    the base in the C standard way, i.e.  0xhh...h means base 16,
  5.    0oo...o means base 8, otherwise assume base 10.
  6. Copyright 1991, 1993, 1994, 1996, 1997, 1998, 2000, 2001, 2002, 2003, 2005
  7. Free Software Foundation, Inc.
  8. This file is part of the GNU MP Library.
  9. The GNU MP Library is free software; you can redistribute it and/or modify
  10. it under the terms of the GNU Lesser General Public License as published by
  11. the Free Software Foundation; either version 3 of the License, or (at your
  12. option) any later version.
  13. The GNU MP Library is distributed in the hope that it will be useful, but
  14. WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  15. or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  16. License for more details.
  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. #include <string.h>
  20. #include <ctype.h>
  21. #include "gmp.h"
  22. #include "gmp-impl.h"
  23. extern const unsigned char __gmp_digit_value_tab[];
  24. #define digit_value_tab __gmp_digit_value_tab
  25. int
  26. mpz_set_str (mpz_ptr x, const char *str, int base)
  27. {
  28.   size_t str_size;
  29.   char *s, *begs;
  30.   size_t i;
  31.   mp_size_t xsize;
  32.   int c;
  33.   int negative;
  34.   const unsigned char *digit_value;
  35.   TMP_DECL;
  36.   digit_value = digit_value_tab;
  37.   if (base > 36)
  38.     {
  39.       /* For bases > 36, use the collating sequence
  40.  0123456789ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz.  */
  41.       digit_value += 224;
  42.       if (base > 62)
  43. return -1; /* too large base */
  44.     }
  45.   /* Skip whitespace.  */
  46.   do
  47.     c = (unsigned char) *str++;
  48.   while (isspace (c));
  49.   negative = 0;
  50.   if (c == '-')
  51.     {
  52.       negative = 1;
  53.       c = (unsigned char) *str++;
  54.     }
  55.   if (digit_value[c] >= (base == 0 ? 10 : base))
  56.     return -1; /* error if no valid digits */
  57.   /* If BASE is 0, try to find out the base by looking at the initial
  58.      characters.  */
  59.   if (base == 0)
  60.     {
  61.       base = 10;
  62.       if (c == '0')
  63. {
  64.   base = 8;
  65.   c = (unsigned char) *str++;
  66.   if (c == 'x' || c == 'X')
  67.     {
  68.       base = 16;
  69.       c = (unsigned char) *str++;
  70.     }
  71.   else if (c == 'b' || c == 'B')
  72.     {
  73.       base = 2;
  74.       c = (unsigned char) *str++;
  75.     }
  76. }
  77.     }
  78.   /* Skip leading zeros and white space.  */
  79.   while (c == '0' || isspace (c))
  80.     c = (unsigned char) *str++;
  81.   /* Make sure the string does not become empty, mpn_set_str would fail.  */
  82.   if (c == 0)
  83.     {
  84.       x->_mp_size = 0;
  85.       return 0;
  86.     }
  87.   TMP_MARK;
  88.   str_size = strlen (str - 1);
  89.   s = begs = (char *) TMP_ALLOC (str_size + 1);
  90.   /* Remove spaces from the string and convert the result from ASCII to a
  91.      byte array.  */
  92.   for (i = 0; i < str_size; i++)
  93.     {
  94.       if (!isspace (c))
  95. {
  96.   int dig = digit_value[c];
  97.   if (dig >= base)
  98.     {
  99.       TMP_FREE;
  100.       return -1;
  101.     }
  102.   *s++ = dig;
  103. }
  104.       c = (unsigned char) *str++;
  105.     }
  106.   str_size = s - begs;
  107.   xsize = 2 + (mp_size_t)
  108.     (str_size / (GMP_NUMB_BITS * mp_bases[base].chars_per_bit_exactly));
  109.   MPZ_REALLOC (x, xsize);
  110.   /* Convert the byte array in base BASE to our bignum format.  */
  111.   xsize = mpn_set_str (x->_mp_d, (unsigned char *) begs, str_size, base);
  112.   x->_mp_size = negative ? -xsize : xsize;
  113.   TMP_FREE;
  114.   return 0;
  115. }