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

数学计算

开发平台:

Unix_Linux

  1. /* mpf_set_str (dest, string, base) -- Convert the string STRING
  2.    in base BASE to a float in dest.  If BASE is zero, the leading characters
  3.    of STRING is used to figure out the base.
  4. Copyright 1993, 1994, 1995, 1996, 1997, 2000, 2001, 2002, 2003, 2005, 2007,
  5. 2008 Free Software Foundation, Inc.
  6. This file is part of the GNU MP Library.
  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. The GNU MP Library is distributed in the hope that it will be useful, but
  12. WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  13. or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  14. License for more details.
  15. You should have received a copy of the GNU Lesser General Public License
  16. along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
  17. /*
  18.   This still needs work, as suggested by some FIXME comments.
  19.   1. Don't depend on superfluous mantissa digits.
  20.   2. Allocate temp space more cleverly.
  21.   3. Use mpn_tdiv_qr instead of mpn_lshift+mpn_divrem.
  22. */
  23. #define _GNU_SOURCE    /* for DECIMAL_POINT in langinfo.h */
  24. #include "config.h"
  25. #include <stdlib.h>
  26. #include <string.h>
  27. #include <ctype.h>
  28. #if HAVE_LANGINFO_H
  29. #include <langinfo.h>  /* for nl_langinfo */
  30. #endif
  31. #if HAVE_LOCALE_H
  32. #include <locale.h>    /* for localeconv */
  33. #endif
  34. #include "gmp.h"
  35. #include "gmp-impl.h"
  36. #include "longlong.h"
  37. extern const unsigned char __gmp_digit_value_tab[];
  38. #define digit_value_tab __gmp_digit_value_tab
  39. /* Compute base^exp and return the most significant prec limbs in rp[].
  40.    Put the count of omitted low limbs in *ign.
  41.    Return the actual size (which might be less than prec).  */
  42. static mp_size_t
  43. mpn_pow_1_highpart (mp_ptr rp, mp_size_t *ignp,
  44.     mp_limb_t base, mp_exp_t exp,
  45.     mp_size_t prec, mp_ptr tp)
  46. {
  47.   mp_size_t ign; /* counts number of ignored low limbs in r */
  48.   mp_size_t off; /* keeps track of offset where value starts */
  49.   mp_ptr passed_rp = rp;
  50.   mp_size_t rn;
  51.   int cnt;
  52.   int i;
  53.   rp[0] = base;
  54.   rn = 1;
  55.   off = 0;
  56.   ign = 0;
  57.   count_leading_zeros (cnt, exp);
  58.   for (i = GMP_LIMB_BITS - cnt - 2; i >= 0; i--)
  59.     {
  60.       mpn_sqr (tp, rp + off, rn);
  61.       rn = 2 * rn;
  62.       rn -= tp[rn - 1] == 0;
  63.       ign <<= 1;
  64.       off = 0;
  65.       if (rn > prec)
  66. {
  67.   ign += rn - prec;
  68.   off = rn - prec;
  69.   rn = prec;
  70. }
  71.       MP_PTR_SWAP (rp, tp);
  72.       if (((exp >> i) & 1) != 0)
  73. {
  74.   mp_limb_t cy;
  75.   cy = mpn_mul_1 (rp, rp + off, rn, base);
  76.   rp[rn] = cy;
  77.   rn += cy != 0;
  78.   off = 0;
  79. }
  80.     }
  81.   if (rn > prec)
  82.     {
  83.       ign += rn - prec;
  84.       rp += rn - prec;
  85.       rn = prec;
  86.     }
  87.   MPN_COPY_INCR (passed_rp, rp + off, rn);
  88.   *ignp = ign;
  89.   return rn;
  90. }
  91. int
  92. mpf_set_str (mpf_ptr x, const char *str, int base)
  93. {
  94.   size_t str_size;
  95.   char *s, *begs;
  96.   size_t i, j;
  97.   int c;
  98.   int negative;
  99.   char *dotpos = 0;
  100.   const char *expptr;
  101.   int exp_base;
  102.   const char  *point = GMP_DECIMAL_POINT;
  103.   size_t      pointlen = strlen (point);
  104.   const unsigned char *digit_value;
  105.   TMP_DECL;
  106.   c = (unsigned char) *str;
  107.   /* Skip whitespace.  */
  108.   while (isspace (c))
  109.     c = (unsigned char) *++str;
  110.   negative = 0;
  111.   if (c == '-')
  112.     {
  113.       negative = 1;
  114.       c = (unsigned char) *++str;
  115.     }
  116.   /* Default base to decimal.  */
  117.   if (base == 0)
  118.     base = 10;
  119.   exp_base = base;
  120.   if (base < 0)
  121.     {
  122.       exp_base = 10;
  123.       base = -base;
  124.     }
  125.   digit_value = digit_value_tab;
  126.   if (base > 36)
  127.     {
  128.       /* For bases > 36, use the collating sequence
  129.  0123456789ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz.  */
  130.       digit_value += 224;
  131.       if (base > 62)
  132. return -1; /* too large base */
  133.     }
  134.   /* Require at least one digit, possibly after an initial decimal point.  */
  135.   if (digit_value[c] >= (base == 0 ? 10 : base))
  136.     {
  137.       /* not a digit, must be a decimal point */
  138.       for (i = 0; i < pointlen; i++)
  139.         if (str[i] != point[i])
  140.           return -1;
  141.       if (digit_value[(unsigned char) str[pointlen]] >= (base == 0 ? 10 : base))
  142. return -1;
  143.     }
  144.   /* Locate exponent part of the input.  Look from the right of the string,
  145.      since the exponent is usually a lot shorter than the mantissa.  */
  146.   expptr = NULL;
  147.   str_size = strlen (str);
  148.   for (i = str_size - 1; i > 0; i--)
  149.     {
  150.       c = (unsigned char) str[i];
  151.       if (c == '@' || (base <= 10 && (c == 'e' || c == 'E')))
  152. {
  153.   expptr = str + i + 1;
  154.   str_size = i;
  155.   break;
  156. }
  157.     }
  158.   TMP_MARK;
  159.   s = begs = (char *) TMP_ALLOC (str_size + 1);
  160.   /* Loop through mantissa, converting it from ASCII to raw byte values.  */
  161.   for (i = 0; i < str_size; i++)
  162.     {
  163.       c = (unsigned char) *str;
  164.       if (!isspace (c))
  165. {
  166.   int dig;
  167.           for (j = 0; j < pointlen; j++)
  168.             if (str[j] != point[j])
  169.               goto not_point;
  170.           if (1)
  171.     {
  172.       if (dotpos != 0)
  173. {
  174.   /* already saw a decimal point, another is invalid */
  175.   TMP_FREE;
  176.   return -1;
  177. }
  178.       dotpos = s;
  179.       str += pointlen - 1;
  180.       i += pointlen - 1;
  181.     }
  182.   else
  183.     {
  184.             not_point:
  185.       dig = digit_value[c];
  186.       if (dig >= base)
  187. {
  188.   TMP_FREE;
  189.   return -1;
  190. }
  191.       *s++ = dig;
  192.     }
  193. }
  194.       c = (unsigned char) *++str;
  195.     }
  196.   str_size = s - begs;
  197.   {
  198.     long exp_in_base;
  199.     mp_size_t ra, ma, rn, mn;
  200.     int cnt;
  201.     mp_ptr mp, tp, rp;
  202.     mp_exp_t exp_in_limbs;
  203.     mp_size_t prec = PREC(x) + 1;
  204.     int divflag;
  205.     mp_size_t madj, radj;
  206. #if 0
  207.     size_t n_chars_needed;
  208.     /* This breaks things like 0.000...0001.  To safely ignore superfluous
  209.        digits, we need to skip over leading zeros.  */
  210.     /* Just consider the relevant leading digits of the mantissa.  */
  211.     n_chars_needed = 2 + (size_t)
  212.       (((size_t) prec * GMP_NUMB_BITS) * mp_bases[base].chars_per_bit_exactly);
  213.     if (str_size > n_chars_needed)
  214.       str_size = n_chars_needed;
  215. #endif
  216.     ma = 2 + (mp_size_t)
  217.       (str_size / (GMP_NUMB_BITS * mp_bases[base].chars_per_bit_exactly));
  218.     mp = TMP_ALLOC_LIMBS (ma);
  219.     mn = mpn_set_str (mp, (unsigned char *) begs, str_size, base);
  220.     if (mn == 0)
  221.       {
  222. SIZ(x) = 0;
  223. EXP(x) = 0;
  224. TMP_FREE;
  225. return 0;
  226.       }
  227.     madj = 0;
  228.     /* Ignore excess limbs in MP,MSIZE.  */
  229.     if (mn > prec)
  230.       {
  231. madj = mn - prec;
  232. mp += mn - prec;
  233. mn = prec;
  234.       }
  235.     if (expptr != 0)
  236.       {
  237. /* Scan and convert the exponent, in base exp_base.  */
  238. long dig, minus, plusminus;
  239. c = (unsigned char) *expptr;
  240. minus = -(long) (c == '-');
  241. plusminus = minus | -(long) (c == '+');
  242. expptr -= plusminus; /* conditional increment */
  243. c = (unsigned char) *expptr++;
  244. dig = digit_value[c];
  245. if (dig >= exp_base)
  246.   {
  247.     TMP_FREE;
  248.     return -1;
  249.   }
  250. exp_in_base = dig;
  251. c = (unsigned char) *expptr++;
  252. dig = digit_value[c];
  253. while (dig < exp_base)
  254.   {
  255.     exp_in_base = exp_in_base * exp_base;
  256.     exp_in_base += dig;
  257.     c = (unsigned char) *expptr++;
  258.     dig = digit_value[c];
  259.   }
  260. exp_in_base = (exp_in_base ^ minus) - minus; /* conditional negation */
  261.       }
  262.     else
  263.       exp_in_base = 0;
  264.     if (dotpos != 0)
  265.       exp_in_base -= s - dotpos;
  266.     divflag = exp_in_base < 0;
  267.     exp_in_base = ABS (exp_in_base);
  268.     if (exp_in_base == 0)
  269.       {
  270. MPN_COPY (PTR(x), mp, mn);
  271. SIZ(x) = negative ? -mn : mn;
  272. EXP(x) = mn + madj;
  273. TMP_FREE;
  274. return 0;
  275.       }
  276.     ra = 2 * (prec + 1);
  277.     rp = TMP_ALLOC_LIMBS (ra);
  278.     tp = TMP_ALLOC_LIMBS (ra);
  279.     rn = mpn_pow_1_highpart (rp, &radj, (mp_limb_t) base, exp_in_base, prec, tp);
  280.     if (divflag)
  281.       {
  282. #if 0
  283. /* FIXME: Should use mpn_tdiv here.  */
  284. mpn_tdiv_qr (qp, mp, 0L, mp, mn, rp, rn);
  285. #else
  286. mp_ptr qp;
  287. mp_limb_t qlimb;
  288. if (mn < rn)
  289.   {
  290.     /* Pad out MP,MSIZE for current divrem semantics.  */
  291.     mp_ptr tmp = TMP_ALLOC_LIMBS (rn + 1);
  292.     MPN_ZERO (tmp, rn - mn);
  293.     MPN_COPY (tmp + rn - mn, mp, mn);
  294.     mp = tmp;
  295.     madj -= rn - mn;
  296.     mn = rn;
  297.   }
  298. if ((rp[rn - 1] & GMP_NUMB_HIGHBIT) == 0)
  299.   {
  300.     mp_limb_t cy;
  301.     count_leading_zeros (cnt, rp[rn - 1]);
  302.     cnt -= GMP_NAIL_BITS;
  303.     mpn_lshift (rp, rp, rn, cnt);
  304.     cy = mpn_lshift (mp, mp, mn, cnt);
  305.     if (cy)
  306.       mp[mn++] = cy;
  307.   }
  308. qp = TMP_ALLOC_LIMBS (prec + 1);
  309. qlimb = mpn_divrem (qp, prec - (mn - rn), mp, mn, rp, rn);
  310. tp = qp;
  311. exp_in_limbs = qlimb + (mn - rn) + (madj - radj);
  312. rn = prec;
  313. if (qlimb != 0)
  314.   {
  315.     tp[prec] = qlimb;
  316.     /* Skip the least significant limb not to overrun the destination
  317.        variable.  */
  318.     tp++;
  319.   }
  320. #endif
  321.       }
  322.     else
  323.       {
  324. tp = TMP_ALLOC_LIMBS (rn + mn);
  325. if (rn > mn)
  326.   mpn_mul (tp, rp, rn, mp, mn);
  327. else
  328.   mpn_mul (tp, mp, mn, rp, rn);
  329. rn += mn;
  330. rn -= tp[rn - 1] == 0;
  331. exp_in_limbs = rn + madj + radj;
  332. if (rn > prec)
  333.   {
  334.     tp += rn - prec;
  335.     rn = prec;
  336.     exp_in_limbs += 0;
  337.   }
  338.       }
  339.     MPN_COPY (PTR(x), tp, rn);
  340.     SIZ(x) = negative ? -rn : rn;
  341.     EXP(x) = exp_in_limbs;
  342.     TMP_FREE;
  343.     return 0;
  344.   }
  345. }