diff options
author | Roland McGrath <roland@gnu.org> | 1996-05-08 02:07:47 +0000 |
---|---|---|
committer | Roland McGrath <roland@gnu.org> | 1996-05-08 02:07:47 +0000 |
commit | 0e3426bbcf2ff61d06d580fc9362fde79953a281 (patch) | |
tree | 6bc37583e3ea48a41eb273fbb112613b51cf64e6 /stdlib/strtod.c | |
parent | b928942eaa0401de9189a2709188f584425f5ca6 (diff) | |
download | glibc-0e3426bbcf2ff61d06d580fc9362fde79953a281.tar.gz glibc-0e3426bbcf2ff61d06d580fc9362fde79953a281.tar.xz glibc-0e3426bbcf2ff61d06d580fc9362fde79953a281.zip |
Tue May 7 19:00:01 1996 Roland McGrath <roland@delasyd.gnu.ai.mit.edu>
* string/argz-extract.c: Remove const from decl. * string/argz.h: Here too. * Makeconfig (version.mk): Fix regexp in sed cmd. Depend on $(..)Makeconfig. * GMP code updated from gmp-2.0 release. * stdlib/Makefile (mpn-routines): Removed add_1, added inlines. * sysdeps/generic/add_1.c: File removed. * stdlib/strtod.c: mp_limb is now mp_limb_t. * stdlib/fpioconst.c, stdlib/fpioconst.h: Likewise. * stdio-common/_itoa.c: Likewise. * stdio-common/printf_fp.c: Likewise. Don't include ansidecl.h. * sysdeps/mach/hurd/getcwd.c: Use io_identity instead of io_stat. * shlib-versions: New file. * Makerules (soversions.mk): New target, include file generated from shlib-versions. Moved shared library rules before installation rules. Rewrote shared library installation rules for versioned libraries. * math/Makefile (libm.so-version): Variable removed. * sysdeps/mach/hurd/i386/exc2signal.c: Use struct hurd_signal_detail. * hurd/report-wait.c (_S_msg_describe_ports): New function. * configure.in: Add AC_PROG_LN_S check. * config.make.in (LN_S): New variable. Sun May 5 03:10:44 1996 Ulrich Drepper <drepper@cygnus.com> * misc/efgcvt_r.c (ecvt_r): Work aroung gcc bug. gcc does not know about weak aliases now and optimizes necessary `if' statement away. * posix/unistd.h: Add swapoff prototype. * sysdeps/generic/confname.h: Add even more POSIX.4 symbols. * sysdeps/posix/fpathconf.c (__fpathconf): Get information for _PC_PATH_MAX from fstatfs function if available. * sysdeps/posix/sysconf.c: Add code to handle _SC_AIO_LISTIO_MAX, _SC_AIO_MAX, _SC_AIO_PRIO_DELTA_MAX, _SC_DELAYTIMER_MAX, _SC_MQ_OPEN_MAX, _SC_MQ_PRIO_MAX, _SC_RTSIG_MAX, _SC_SEM_NSEMS_MAX, _SC_SEM_VALUE_MAX, _SC_SIGQUEUE_MAX, and _SC_TIMER_MAX. * sysdeps/unix/sysv/sysv4/sysconf.c: Ditto. * sysdeps/stub/swapoff.c: New file. Stub version for swapoff function. * sysdeps/unix/syscalls.list: Add swapoff. * sysdeps/unix/sysv/linux/Dist: Add sys/acct.h. * sysdeps/unix/sysv/linux/Makefile [$(subdir) == misc] (sysdep_routines): Add mount, umount, llseek, setfsgid, setfsuid, sysinfo, and uselib. (headers): Add sys/sysinfo.h. * sysdeps/unix/sysv/linux/gethostid.c: Prevent warning. * sysdeps/unix/sysv/linux/i386/Makefile [$(subdir) == misc] (sysdep_routines): Add ioperm, iopl, and vm86. (headers): Add sys/perm.h and sys/vm86.h. * sysdeps/unix/sysv/linux/i386/sys/perm.h: New file. Contains prototypes for iopl and ioperm. * sysdeps/unix/sysv/linux/i386/sys/vm86.h: New file. Contains prototype for vm86. * sysdeps/unix/sysv/linux/i386/syscalls.list: New file. Add vm86 system call. * sysdeps/unix/sysv/linux/sys/acct.h: New file. Contains prototypes for acct function. * sysdeps/unix/sysv/linux/sys/socket.h: Provide real header file with prototypes. * sysdeps/unix/sysv/linux/sys/sysinfo.h: New file. Contains prototype for sysinfo function. * sysdeps/unix/sysv/linux/syscalls.list: Add flock, ioperm, iopl, llseek, setfsgid, setfsuid, sysinfo, and uselib. * sysdeps/unix/sysv/linux/sysconf.c: Instead of duplicating posix/sysconf.c now only handle cases different to that implementation. Tue May 7 15:08:19 1996 Miles Bader <miles@gnu.ai.mit.edu> * stdio/linewrap.c (__line_wrap_output): Renamed from lwoutput (all references changed). Now exported. * stdio/linewrap.c (struct data): Type deleted (moved to linewrap.h). (wrap_stream, unwrap_stream, lwclose, lwfileno, lwoutput, line_wrap_stream, line_unwrap_stream): Use struct line_wrap_data instead of struct data. (lwoutput, line_wrap_stream, line_unwrap_stream): Rename various occurences of `wrap' and `wrapmargin' to `wmargin'. (line_wrapped, line_wrap_lmargin, line_wrap_set_lmargin, line_wrap_rmargin, line_wrap_set_rmargin, line_wrap_wmargin, line_wrap_set_wmargin, line_wrap_point): New functions. * stdio/linewrap.h: New file. * stdio/Makefile (headers): Add linewrap.h. Tue May 7 14:19:12 1996 Roland McGrath <roland@delasyd.gnu.ai.mit.edu> * sysdeps/unix/sysv/linux/i386/Makefile: File removed. * stdio/stdio.h: Remove line_wrap_stream, line_unwap_stream decls. * sysdeps/unix/sysv/linux/schedbits.h: New file. Tue May 7 13:47:02 1996 Miles Bader <miles@gnu.ai.mit.edu> * stdio/linewrap.c (struct data): Make margin fields not-pointers. (lwoutput): Adjust uses acordingly. * sysdeps/mach/hurd/fdatasync.c: New file. * sysdeps/mach/hurd/fsync.c: Pass new flag to file_sync. * sysdeps/mach/hurd/xmknod.c: Pass new flag to dir_link. * sysdeps/mach/hurd/symlink.c: Likewise. * sysdeps/mach/hurd/link.c: Likewise. * sysdeps/mach/hurd/bind.c: Likewise. * hurd/hurdsig.c (write_corefile): Likewise. * hurd/hurdsig.c (write_corefile): Pass cttyid port to crash server. * sysdeps/mach/hurd/fpathconf.c: RPC takes int pointer, not long int. * sysdeps/mach/hurd/_exit.c (_hurd_exit): Pass sigcode arg to proc_mark_exit. * sysdeps/mach/hurd/dl-sysdep.c (_exit): Likewise. * sysdeps/mach/hurd/wait4.c: Pass sigcode arg to proc_wait. * sysdeps/mach/hurd/rename.c: Pass new flag to dir_rename. * hurd/hurdfault.c (_hurdsig_fault_catch_exception_raise): Use struct hurd_signal_detail. * hurd/catch-exc.c (_S_catch_exception_raise): Likewise. * hurd/hurd-raise.c (_hurd_raise_signal): Likewise. * sysdeps/mach/hurd/i386/trampoline.c (_hurd_setup_sighandler): Likewise. * sysdeps/mach/hurd/setitimer.c (restart_itimer): Likewise. * hurd/hurd/signal.h: Fix _hurd_exception2signal prototype. * hurd/hurdsig.c (write_corefile): Take const struct hurd_signal_detail * arg. Pass all details to crash_dump_task. (_hurd_internal_post_signal): Pass DETAIL to write_corefile. (_hurd_internal_post_signal: suspend): Pass code and error to proc_mark_stop. * hurd/hurdprio.c (_hurd_priority_which_map): Pass flags arg to proc_getprocinfo by reference.
Diffstat (limited to 'stdlib/strtod.c')
-rw-r--r-- | stdlib/strtod.c | 90 |
1 files changed, 46 insertions, 44 deletions
diff --git a/stdlib/strtod.c b/stdlib/strtod.c index e262f3de83..525a6f98a9 100644 --- a/stdlib/strtod.c +++ b/stdlib/strtod.c @@ -95,12 +95,12 @@ extern FLOAT MPN2FLOAT (mp_srcptr mpn, int exponent, int negative); # define MAX_DIG_PER_LIMB 19 # define MAX_FAC_PER_LIMB 10000000000000000000UL #else -# error "mp_limb size " BITS_PER_MP_LIMB "not accounted for" +# error "mp_limb_t size " BITS_PER_MP_LIMB "not accounted for" #endif /* Local data structure. */ -static const mp_limb _tens_in_limb[MAX_DIG_PER_LIMB + 1] = +static const mp_limb_t _tens_in_limb[MAX_DIG_PER_LIMB + 1] = { 0, 10, 100, 1000, 10000, 100000, 1000000, 10000000, 100000000, @@ -132,17 +132,17 @@ static const mp_limb _tens_in_limb[MAX_DIG_PER_LIMB + 1] = #define MPNSIZE (howmany (MAX_EXP + 2 * MANT_DIG, BITS_PER_MP_LIMB) \ + 2) /* Declare an mpn integer variable that big. */ -#define MPN_VAR(name) mp_limb name[MPNSIZE]; mp_size_t name##size +#define MPN_VAR(name) mp_limb_t name[MPNSIZE]; mp_size_t name##size /* Copy an mpn integer value. */ #define MPN_ASSIGN(dst, src) \ - memcpy (dst, src, (dst##size = src##size) * sizeof (mp_limb)) + memcpy (dst, src, (dst##size = src##size) * sizeof (mp_limb_t)) /* Return a floating point number of the needed type according to the given multi-precision number after possible rounding. */ static inline FLOAT -round_and_return (mp_limb *retval, int exponent, int negative, - mp_limb round_limb, mp_size_t round_bit, int more_bits) +round_and_return (mp_limb_t *retval, int exponent, int negative, + mp_limb_t round_limb, mp_size_t round_bit, int more_bits) { if (exponent < MIN_EXP - 1) { @@ -154,7 +154,7 @@ round_and_return (mp_limb *retval, int exponent, int negative, return 0.0; } - more_bits |= (round_limb & ((((mp_limb) 1) << round_bit) - 1)) != 0; + more_bits |= (round_limb & ((((mp_limb_t) 1) << round_bit) - 1)) != 0; if (shift == MANT_DIG) /* This is a special case to handle the very seldom case where the mantissa will be empty after the shift. */ @@ -175,7 +175,8 @@ round_and_return (mp_limb *retval, int exponent, int negative, round_bit = (shift - 1) % BITS_PER_MP_LIMB; for (i = 0; i < (shift - 1) / BITS_PER_MP_LIMB; ++i) more_bits |= retval[i] != 0; - more_bits |= (round_limb & ((((mp_limb) 1) << round_bit) - 1)) != 0; + more_bits |= ((round_limb & ((((mp_limb_t) 1) << round_bit) - 1)) + != 0); (void) __mpn_rshift (retval, &retval[shift / BITS_PER_MP_LIMB], RETURN_LIMB_SIZE - (shift / BITS_PER_MP_LIMB), @@ -192,25 +193,25 @@ round_and_return (mp_limb *retval, int exponent, int negative, exponent = MIN_EXP - 2; } - if ((round_limb & (((mp_limb) 1) << round_bit)) != 0 + if ((round_limb & (((mp_limb_t) 1) << round_bit)) != 0 && (more_bits || (retval[0] & 1) != 0 - || (round_limb & ((((mp_limb) 1) << round_bit) - 1)) != 0)) + || (round_limb & ((((mp_limb_t) 1) << round_bit) - 1)) != 0)) { - mp_limb cy = __mpn_add_1 (retval, retval, RETURN_LIMB_SIZE, 1); + mp_limb_t cy = __mpn_add_1 (retval, retval, RETURN_LIMB_SIZE, 1); if (((MANT_DIG % BITS_PER_MP_LIMB) == 0 && cy) || ((MANT_DIG % BITS_PER_MP_LIMB) != 0 && (retval[RETURN_LIMB_SIZE - 1] - & (((mp_limb) 1) << (MANT_DIG % BITS_PER_MP_LIMB))) != 0)) + & (((mp_limb_t) 1) << (MANT_DIG % BITS_PER_MP_LIMB))) != 0)) { ++exponent; (void) __mpn_rshift (retval, retval, RETURN_LIMB_SIZE, 1); retval[RETURN_LIMB_SIZE - 1] - |= ((mp_limb) 1) << ((MANT_DIG - 1) % BITS_PER_MP_LIMB); + |= ((mp_limb_t) 1) << ((MANT_DIG - 1) % BITS_PER_MP_LIMB); } else if (exponent == MIN_EXP - 2 && (retval[RETURN_LIMB_SIZE - 1] - & (((mp_limb) 1) << ((MANT_DIG - 1) % BITS_PER_MP_LIMB))) + & (((mp_limb_t) 1) << ((MANT_DIG - 1) % BITS_PER_MP_LIMB))) != 0) /* The number was denormalized but now normalized. */ exponent = MIN_EXP - 1; @@ -229,13 +230,13 @@ round_and_return (mp_limb *retval, int exponent, int negative, value. If the EXPONENT is small enough to be taken as an additional factor for the resulting number (see code) multiply by it. */ static inline const STRING_TYPE * -str_to_mpn (const STRING_TYPE *str, int digcnt, mp_limb *n, mp_size_t *nsize, +str_to_mpn (const STRING_TYPE *str, int digcnt, mp_limb_t *n, mp_size_t *nsize, int *exponent) { /* Number of digits for actual limb. */ int cnt = 0; - mp_limb low = 0; - mp_limb base; + mp_limb_t low = 0; + mp_limb_t base; *nsize = 0; assert (digcnt > 0); @@ -247,7 +248,7 @@ str_to_mpn (const STRING_TYPE *str, int digcnt, mp_limb *n, mp_size_t *nsize, n[0] = low; else { - mp_limb cy; + mp_limb_t cy; cy = __mpn_mul_1 (n, n, *nsize, MAX_FAC_PER_LIMB); cy += __mpn_add_1 (n, n, *nsize, low); if (cy != 0) @@ -284,7 +285,7 @@ str_to_mpn (const STRING_TYPE *str, int digcnt, mp_limb *n, mp_size_t *nsize, } else { - mp_limb cy; + mp_limb_t cy; cy = __mpn_mul_1 (n, n, *nsize, base); cy += __mpn_add_1 (n, n, *nsize, low); if (cy != 0) @@ -300,7 +301,8 @@ str_to_mpn (const STRING_TYPE *str, int digcnt, mp_limb *n, mp_size_t *nsize, Tege doesn't like this function so I have to write it here myself. :) --drepper */ static inline void -__mpn_lshift_1 (mp_limb *ptr, mp_size_t size, unsigned int count, mp_limb limb) +__mpn_lshift_1 (mp_limb_t *ptr, mp_size_t size, unsigned int count, + mp_limb_t limb) { if (count == BITS_PER_MP_LIMB) { @@ -347,7 +349,7 @@ INTERNAL (STRTOF) (nptr, endptr, group) MPN_VAR (den); /* Representation for the return value. */ - mp_limb retval[RETURN_LIMB_SIZE]; + mp_limb_t retval[RETURN_LIMB_SIZE]; /* Number of bits currently in result value. */ int bits; @@ -638,8 +640,8 @@ INTERNAL (STRTOF) (nptr, endptr, group) if (exponent > 0) { /* We now multiply the gained number by the given power of ten. */ - mp_limb *psrc = num; - mp_limb *pdest = den; + mp_limb_t *psrc = num; + mp_limb_t *pdest = den; int expbit = 1; const struct mp_power *ttab = &_fpioconst_pow10[0]; @@ -647,7 +649,7 @@ INTERNAL (STRTOF) (nptr, endptr, group) { if ((exponent & expbit) != 0) { - mp_limb cy; + mp_limb_t cy; exponent ^= expbit; /* FIXME: not the whole multiplication has to be done. @@ -672,7 +674,7 @@ INTERNAL (STRTOF) (nptr, endptr, group) while (exponent != 0); if (psrc == den) - memcpy (num, den, numsize * sizeof (mp_limb)); + memcpy (num, den, numsize * sizeof (mp_limb_t)); } /* Determine how many bits of the result we already have. */ @@ -702,7 +704,7 @@ INTERNAL (STRTOF) (nptr, endptr, group) if (least_bit == 0) memcpy (retval, &num[least_idx], - RETURN_LIMB_SIZE * sizeof (mp_limb)); + RETURN_LIMB_SIZE * sizeof (mp_limb_t)); else { for (i = least_idx; i < numsize - 1; ++i) @@ -730,7 +732,7 @@ INTERNAL (STRTOF) (nptr, endptr, group) if (target_bit == is_bit) { memcpy (&retval[RETURN_LIMB_SIZE - numsize], num, - numsize * sizeof (mp_limb)); + numsize * sizeof (mp_limb_t)); /* FIXME: the following loop can be avoided if we assume a maximal MANT_DIG value. */ MPN_ZERO (retval, RETURN_LIMB_SIZE - numsize); @@ -745,7 +747,7 @@ INTERNAL (STRTOF) (nptr, endptr, group) } else { - mp_limb cy; + mp_limb_t cy; assert (numsize < RETURN_LIMB_SIZE); cy = __mpn_rshift (&retval[RETURN_LIMB_SIZE - numsize], @@ -761,7 +763,7 @@ INTERNAL (STRTOF) (nptr, endptr, group) } /* Store the bits we already have. */ - memcpy (retval, num, numsize * sizeof (mp_limb)); + memcpy (retval, num, numsize * sizeof (mp_limb_t)); #if RETURN_LIMB_SIZE > 1 if (numsize < RETURN_LIMB_SIZE) retval[numsize] = 0; @@ -779,9 +781,9 @@ INTERNAL (STRTOF) (nptr, endptr, group) int cnt; int neg_exp; int more_bits; - mp_limb cy; - mp_limb *psrc = den; - mp_limb *pdest = num; + mp_limb_t cy; + mp_limb_t *psrc = den; + mp_limb_t *pdest = num; const struct mp_power *ttab = &_fpioconst_pow10[0]; assert (dig_no > int_no && exponent <= 0); @@ -810,14 +812,14 @@ INTERNAL (STRTOF) (nptr, endptr, group) { if ((neg_exp & expbit) != 0) { - mp_limb cy; + mp_limb_t cy; neg_exp ^= expbit; if (densize == 0) { densize = ttab->arraysize - _FPIO_CONST_OFFSET; memcpy (psrc, &ttab->array[_FPIO_CONST_OFFSET], - densize * sizeof (mp_limb)); + densize * sizeof (mp_limb_t)); } else { @@ -836,7 +838,7 @@ INTERNAL (STRTOF) (nptr, endptr, group) while (neg_exp != 0); if (psrc == num) - memcpy (den, num, densize * sizeof (mp_limb)); + memcpy (den, num, densize * sizeof (mp_limb_t)); /* Read the fractional digits from the string. */ (void) str_to_mpn (startp, dig_no - int_no, num, &numsize, &exponent); @@ -873,7 +875,7 @@ INTERNAL (STRTOF) (nptr, endptr, group) { case 1: { - mp_limb d, n, quot; + mp_limb_t d, n, quot; int used = 0; n = num[0]; @@ -930,8 +932,8 @@ INTERNAL (STRTOF) (nptr, endptr, group) } case 2: { - mp_limb d0, d1, n0, n1; - mp_limb quot = 0; + mp_limb_t d0, d1, n0, n1; + mp_limb_t quot = 0; int used = 0; d0 = den[0]; @@ -976,14 +978,14 @@ INTERNAL (STRTOF) (nptr, endptr, group) while (bits <= MANT_DIG) { - mp_limb r; + mp_limb_t r; if (n1 == d1) { /* QUOT should be either 111..111 or 111..110. We need special treatment of this rare case as normal division would give overflow. */ - quot = ~(mp_limb) 0; + quot = ~(mp_limb_t) 0; r = n0 + d1; if (r < d1) /* Carry in the addition? */ @@ -1024,8 +1026,8 @@ INTERNAL (STRTOF) (nptr, endptr, group) default: { int i; - mp_limb cy, dX, d1, n0, n1; - mp_limb quot = 0; + mp_limb_t cy, dX, d1, n0, n1; + mp_limb_t quot = 0; int used = 0; dX = den[densize - 1]; @@ -1105,10 +1107,10 @@ INTERNAL (STRTOF) (nptr, endptr, group) if (n0 == dX) /* This might over-estimate QUOT, but it's probably not worth the extra code here to find out. */ - quot = ~(mp_limb) 0; + quot = ~(mp_limb_t) 0; else { - mp_limb r; + mp_limb_t r; udiv_qrnnd (quot, r, n0, num[densize - 1], dX); umul_ppmm (n1, n0, d1, quot); |