=================================================================== RCS file: /home/cvs/OpenXM_contrib/gmp/mpf/Attic/set_str.c,v retrieving revision 1.1.1.1 retrieving revision 1.1.1.3 diff -u -p -r1.1.1.1 -r1.1.1.3 --- OpenXM_contrib/gmp/mpf/Attic/set_str.c 2000/01/10 15:35:22 1.1.1.1 +++ OpenXM_contrib/gmp/mpf/Attic/set_str.c 2003/08/25 16:06:34 1.1.1.3 @@ -2,37 +2,42 @@ in base BASE to a float in dest. If BASE is zero, the leading characters of STRING is used to figure out the base. -Copyright (C) 1993, 1994, 1995, 1996 Free Software Foundation, Inc. +Copyright 1993, 1994, 1995, 1996, 1997, 2000, 2001, 2002 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 Library General Public License as published by -the Free Software Foundation; either version 2 of the License, or (at your +it under the terms of the GNU Lesser General Public License as published by +the Free Software Foundation; either version 2.1 of the License, 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 Library General Public +or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details. -You should have received a copy of the GNU Library General Public License +You should have received a copy of the GNU Lesser General Public License along with the GNU MP Library; see the file COPYING.LIB. If not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ +#include "config.h" + +#include #include #include + +#if HAVE_LOCALE_H +#include /* for localeconv */ +#endif + #include "gmp.h" #include "gmp-impl.h" #include "longlong.h" -long int strtol _PROTO ((const char *, char **ptr, int)); - static int -digit_value_in_base (c, base) - int c; - int base; +digit_value_in_base (int c, int base) { int digit; @@ -50,15 +55,23 @@ digit_value_in_base (c, base) return -1; } -int -#if __STDC__ -mpf_set_str (mpf_ptr x, const char *str, int base) +#if HAVE_LOCALECONV +/* Avoid memcmp for the usual case that point is one character. Don't + bother with str+1,point+1,pointlen-1 since those offsets would add to the + code size. */ +#define POINTCMP_FAST(c, str) \ + ((c) == point0 \ + && (pointlen == 1 || memcmp (str, point, pointlen) == 0)) +#define POINTCMP_SMALL(c, str) (memcmp (str, point, pointlen) == 0) +#define POINTLEN (pointlen) #else -mpf_set_str (x, str, base) - mpf_ptr x; - char *str; - int base; +#define POINTCMP_FAST(c, str) ((c) == '.') +#define POINTCMP_SMALL(c, str) ((c) == '.') +#define POINTLEN 1 #endif + +int +mpf_set_str (mpf_ptr x, const char *str, int base) { size_t str_size; char *s, *begs; @@ -69,6 +82,11 @@ mpf_set_str (x, str, base) char *dotpos = 0; int expflag; int decimal_exponent_flag; +#if HAVE_LOCALECONV + const char *point = localeconv()->decimal_point; + size_t pointlen = strlen (point); + int point0 = point[0]; +#endif TMP_DECL (marker); TMP_MARK (marker); @@ -89,8 +107,14 @@ mpf_set_str (x, str, base) decimal_exponent_flag = base < 0; base = ABS (base); + /* require at least one digit, possibly after an initial decimal point */ if (digit_value_in_base (c, base == 0 ? 10 : base) < 0) - return -1; /* error if no digits */ + { + if (! POINTCMP_SMALL (c, str)) + return -1; + if (digit_value_in_base (str[POINTLEN], base == 0 ? 10 : base) < 0) + return -1; + } /* If BASE is 0, try to find out the base by looking at the initial characters. */ @@ -101,9 +125,12 @@ mpf_set_str (x, str, base) if (c == '0') { base = 8; - c = *++str; - if (c == 'x' || c == 'X') - base = 16; + if (str[1] == 'x' || str[1] == 'X') + { + base = 16; + str += 2; + c = *str; + } } #endif } @@ -131,14 +158,16 @@ mpf_set_str (x, str, base) { int dig; - if (c == '.') - { + if (POINTCMP_FAST (c, str)) + { if (dotpos != 0) { TMP_FREE (marker); return -1; } dotpos = s; + str += POINTLEN - 1; + i += POINTLEN - 1; } else { @@ -159,14 +188,13 @@ mpf_set_str (x, str, base) xsize = str_size / __mp_bases[base].chars_per_limb + 2; { long exp_in_base; - mp_size_t rsize, msize; + mp_size_t ralloc, rsize, msize; int cnt, i; - mp_ptr mp, xp, tp, rp; - mp_limb_t cy; + mp_ptr mp, tp, rp; mp_exp_t exp_in_limbs; - mp_size_t prec = x->_mp_prec; + mp_size_t prec = x->_mp_prec + 1; int divflag; - mp_size_t xxx = 0; + mp_size_t madj, radj; mp = (mp_ptr) TMP_ALLOC (xsize * BYTES_PER_MP_LIMB); msize = mpn_set_str (mp, (unsigned char *) begs, str_size, base); @@ -179,6 +207,15 @@ mpf_set_str (x, str, base) return 0; } + madj = 0; + /* Ignore excess limbs in MP,MSIZE. */ + if (msize > prec) + { + madj = msize - prec; + mp += msize - prec; + msize = prec; + } + if (expflag != 0) exp_in_base = strtol (str + 1, (char **) 0, decimal_exponent_flag ? 10 : base); @@ -193,35 +230,39 @@ mpf_set_str (x, str, base) { MPN_COPY (x->_mp_d, mp, msize); x->_mp_size = negative ? -msize : msize; - x->_mp_exp = msize; + x->_mp_exp = msize + madj; TMP_FREE (marker); return 0; } -#if 1 - rsize = (((mp_size_t) (exp_in_base / __mp_bases[base].chars_per_bit_exactly)) - / BITS_PER_MP_LIMB + 3); -#else - count_leading_zeros (cnt, (mp_limb_t) base); - rsize = exp_in_base - cnt * exp_in_base / BITS_PER_MP_LIMB + 1; -#endif - rp = (mp_ptr) TMP_ALLOC (rsize * BYTES_PER_MP_LIMB); - tp = (mp_ptr) TMP_ALLOC (rsize * BYTES_PER_MP_LIMB); + radj = 0; /* counts number of ignored low limbs in r */ + ralloc = 2 * (prec + 1); + rp = (mp_ptr) TMP_ALLOC (ralloc * BYTES_PER_MP_LIMB); + tp = (mp_ptr) TMP_ALLOC (ralloc * BYTES_PER_MP_LIMB); + rp[0] = base; rsize = 1; - count_leading_zeros (cnt, exp_in_base); - - for (i = BITS_PER_MP_LIMB - cnt - 2; i >= 0; i--) + for (i = GMP_LIMB_BITS - cnt - 2; i >= 0; i--) { - mpn_mul_n (tp, rp, rp, rsize); + mpn_sqr_n (tp, rp, rsize); rsize = 2 * rsize; rsize -= tp[rsize - 1] == 0; - xp = tp; tp = rp; rp = xp; + radj <<= 1; + if (rsize > prec) + { + radj += rsize - prec; + MPN_COPY_INCR (rp, tp + rsize - prec, prec); + rsize = prec; + } + else + MP_PTR_SWAP (rp, tp); + if (((exp_in_base >> i) & 1) != 0) { + mp_limb_t cy; cy = mpn_mul_1 (rp, rp, rsize, (mp_limb_t) base); rp[rsize] = cy; rsize += cy != 0; @@ -230,48 +271,48 @@ mpf_set_str (x, str, base) if (rsize > prec) { - xxx += rsize - prec; + radj += rsize - prec; rp += rsize - prec; rsize = prec; } -#if 0 - if (msize > prec) - { - xxx -= msize - prec; - mp += msize - prec; - msize = prec; - } -#endif + if (divflag) { mp_ptr qp; - mp_limb_t qflag; - mp_size_t xtra; - if (msize <= rsize) + mp_limb_t qlimb; + if (msize < rsize) { - /* Allocate extra limb for current divrem sematics. */ + /* Pad out MP,MSIZE for current divrem semantics. */ mp_ptr tmp = (mp_ptr) TMP_ALLOC ((rsize + 1) * BYTES_PER_MP_LIMB); MPN_ZERO (tmp, rsize - msize); MPN_COPY (tmp + rsize - msize, mp, msize); mp = tmp; - xxx += rsize - msize; + madj -= rsize - msize; msize = rsize; } - count_leading_zeros (cnt, rp[rsize - 1]); - if (cnt != 0) + if ((rp[rsize - 1] & GMP_NUMB_HIGHBIT) == 0) { + mp_limb_t cy; + count_leading_zeros (cnt, rp[rsize - 1]); + cnt -= GMP_NAIL_BITS; mpn_lshift (rp, rp, rsize, cnt); cy = mpn_lshift (mp, mp, msize, cnt); if (cy) mp[msize++] = cy; } + qp = (mp_ptr) TMP_ALLOC ((prec + 1) * BYTES_PER_MP_LIMB); - xtra = prec - (msize - rsize); - qflag = mpn_divrem (qp, xtra, mp, msize, rp, rsize); - qp[prec] = qflag; + qlimb = mpn_divrem (qp, prec - (msize - rsize), mp, msize, rp, rsize); tp = qp; - rsize = prec + qflag; - exp_in_limbs = rsize - xtra - xxx; + exp_in_limbs = qlimb + (msize - rsize) + (madj - radj); + rsize = prec; + if (qlimb != 0) + { + tp[prec] = qlimb; + /* Skip the least significant limb not to overrun the destination + variable. */ + tp++; + } } else { @@ -282,11 +323,10 @@ mpf_set_str (x, str, base) mpn_mul (tp, mp, msize, rp, rsize); rsize += msize; rsize -= tp[rsize - 1] == 0; - exp_in_limbs = rsize + xxx; + exp_in_limbs = rsize + madj + radj; if (rsize > prec) { - xxx = rsize - prec; tp += rsize - prec; rsize = prec; exp_in_limbs += 0;