1e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott/* 2e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Copyright (c) 1999 3e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Silicon Graphics Computer Systems, Inc. 4e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 5e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Copyright (c) 1999 6e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Boris Fomitchev 7e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 8e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * This material is provided "as is", with absolutely no warranty expressed 9e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * or implied. Any use is at your own risk. 10e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 11e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Permission to use or copy this software for any purpose is hereby granted 12e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * without fee, provided the above notices are retained on all copies. 13e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Permission to modify the code and to distribute modified code is granted, 14e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * provided the above notices are retained, and a notice that the code was 15e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * modified is included with the above copyright notice. 16e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 17e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 18e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 19e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#include "stlport_prefix.h" 20e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 21e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#include <limits> 22e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#include <locale> 23e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#include <istream> 24e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 25e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#if (defined (__GNUC__) && !defined (__sun) && !defined (__hpux)) || \ 26e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott defined (__DMC__) 27e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# include <stdint.h> 28e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif 29e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 30e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#if defined (__linux__) || defined (__MINGW32__) || defined (__CYGWIN__) || \ 31e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott defined (__BORLANDC__) || defined (__DMC__) || defined (__HP_aCC) 32e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 33e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# if defined (__BORLANDC__) 34e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scotttypedef unsigned int uint32_t; 35e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scotttypedef unsigned __int64 uint64_t; 36e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# endif 37e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 38e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottunion _ll { 39e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64_t i64; 40e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott struct { 41e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# if defined (_STLP_BIG_ENDIAN) 42e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint32_t hi; 43e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint32_t lo; 44e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# elif defined (_STLP_LITTLE_ENDIAN) 45e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint32_t lo; 46e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint32_t hi; 47e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# else 48e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# error Unknown endianess 49e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# endif 50e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } i32; 51e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott}; 52e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 537924895a54b8f0b67514224a7ad9d67c79b3eeafAndrew Hsieh# if defined (__linux__) && !defined (__ANDROID__) 54e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# include <ieee754.h> 55e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# else 56e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottunion ieee854_long_double { 57e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott long double d; 58e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 59e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* This is the IEEE 854 double-extended-precision format. */ 60e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott struct { 61e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott unsigned int mantissa1:32; 62e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott unsigned int mantissa0:32; 63e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott unsigned int exponent:15; 64e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott unsigned int negative:1; 65e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott unsigned int empty:16; 66e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } ieee; 67e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott}; 68e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 69e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# define IEEE854_LONG_DOUBLE_BIAS 0x3fff 70e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# endif 71e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif 72e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 73e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott_STLP_BEGIN_NAMESPACE 74e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott_STLP_MOVE_TO_PRIV_NAMESPACE 75e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 76e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott//---------------------------------------------------------------------- 77e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// num_get 78e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 79e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// Helper functions for _M_do_get_float. 80e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 81e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#if !defined (_STLP_NO_WCHAR_T) 82e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottvoid _STLP_CALL 83e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott_Initialize_get_float( const ctype<wchar_t>& ct, 84e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott wchar_t& Plus, wchar_t& Minus, 85e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott wchar_t& pow_e, wchar_t& pow_E, 86e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott wchar_t* digits) { 87e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott char ndigits[11] = "0123456789"; 88e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott Plus = ct.widen('+'); 89e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott Minus = ct.widen('-'); 90e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott pow_e = ct.widen('e'); 91e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott pow_E = ct.widen('E'); 92e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott ct.widen(ndigits + 0, ndigits + 10, digits); 93e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 94e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif /* WCHAR_T */ 95e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 96e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott/* 97e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * __string_to_double is just lifted from atof, the difference being 98e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * that we just use '.' for the decimal point, rather than let it 99e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * be taken from the current C locale, which of course is not accessible 100e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * to us. 101e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 102e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#if defined (_STLP_MSVC) || defined (__BORLANDC__) || defined (__ICL) 103e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scotttypedef unsigned long uint32; 104e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scotttypedef unsigned __int64 uint64; 105e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# define ULL(x) x##Ui64 106e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#elif defined (__unix) || defined (__MINGW32__) || \ 107e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott (defined (__DMC__) && (__LONGLONG)) || defined (__WATCOMC__) || \ 1087924895a54b8f0b67514224a7ad9d67c79b3eeafAndrew Hsieh defined (__ANDROID__) 109e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scotttypedef uint32_t uint32; 110e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scotttypedef uint64_t uint64; 111e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# define ULL(x) x##ULL 112e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#else 113e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# error There should be some unsigned 64-bit integer on the system! 114e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif 115e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 116e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// Multiplication of two 64-bit integers, giving a 128-bit result. 117e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// Taken from Algorithm M in Knuth section 4.3.1, with the loop 118e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// hand-unrolled. 119e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottstatic void _Stl_mult64(const uint64 u, const uint64 v, 120e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64& high, uint64& low) { 121e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott const uint64 low_mask = ULL(0xffffffff); 122e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott const uint64 u0 = u & low_mask; 123e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott const uint64 u1 = u >> 32; 124e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott const uint64 v0 = v & low_mask; 125e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott const uint64 v1 = v >> 32; 126e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 127e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64 t = u0 * v0; 128e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott low = t & low_mask; 129e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 130e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott t = u1 * v0 + (t >> 32); 131e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64 w1 = t & low_mask; 132e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64 w2 = t >> 32; 133e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 134e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64 x = u0 * v1 + w1; 135e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott low += (x & low_mask) << 32; 136e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott high = u1 * v1 + w2 + (x >> 32); 137e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 138e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 1397924895a54b8f0b67514224a7ad9d67c79b3eeafAndrew Hsieh#if !defined (__linux__) || defined (__ANDROID__) 140e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 141e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# define bit11 ULL(0x7ff) 142e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# define exponent_mask (bit11 << 52) 143e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 144e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# if !defined (__GNUC__) || (__GNUC__ != 3) || (__GNUC_MINOR__ != 4) || \ 145e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott (!defined (__CYGWIN__) && !defined (__MINGW32__)) 146e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott//Generate bad code when compiled with -O2 option. 147e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottinline 148e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott# endif 149e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottvoid _Stl_set_exponent(uint64 &val, uint64 exp) 150e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott{ val = (val & ~exponent_mask) | ((exp & bit11) << 52); } 151e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 152e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif // __linux__ 153e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 154e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott/* Power of ten fractions for tenscale*/ 155e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott/* The constants are factored so that at most two constants 156e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * and two multiplies are needed. Furthermore, one of the constants 157e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * is represented exactly - 10**n where 1<= n <= 27. 158e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 159e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 160e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottstatic const uint64 _Stl_tenpow[80] = { 161e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xa000000000000000), /* _Stl_tenpow[0]=(10**1)/(2**4) */ 162e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xc800000000000000), /* _Stl_tenpow[1]=(10**2)/(2**7) */ 163e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xfa00000000000000), /* _Stl_tenpow[2]=(10**3)/(2**10) */ 164e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0x9c40000000000000), /* _Stl_tenpow[3]=(10**4)/(2**14) */ 165e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xc350000000000000), /* _Stl_tenpow[4]=(10**5)/(2**17) */ 166e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xf424000000000000), /* _Stl_tenpow[5]=(10**6)/(2**20) */ 167e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0x9896800000000000), /* _Stl_tenpow[6]=(10**7)/(2**24) */ 168e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xbebc200000000000), /* _Stl_tenpow[7]=(10**8)/(2**27) */ 169e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xee6b280000000000), /* _Stl_tenpow[8]=(10**9)/(2**30) */ 170e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0x9502f90000000000), /* _Stl_tenpow[9]=(10**10)/(2**34) */ 171e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xba43b74000000000), /* _Stl_tenpow[10]=(10**11)/(2**37) */ 172e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xe8d4a51000000000), /* _Stl_tenpow[11]=(10**12)/(2**40) */ 173e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0x9184e72a00000000), /* _Stl_tenpow[12]=(10**13)/(2**44) */ 174e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xb5e620f480000000), /* _Stl_tenpow[13]=(10**14)/(2**47) */ 175e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xe35fa931a0000000), /* _Stl_tenpow[14]=(10**15)/(2**50) */ 176e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0x8e1bc9bf04000000), /* _Stl_tenpow[15]=(10**16)/(2**54) */ 177e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xb1a2bc2ec5000000), /* _Stl_tenpow[16]=(10**17)/(2**57) */ 178e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xde0b6b3a76400000), /* _Stl_tenpow[17]=(10**18)/(2**60) */ 179e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0x8ac7230489e80000), /* _Stl_tenpow[18]=(10**19)/(2**64) */ 180e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xad78ebc5ac620000), /* _Stl_tenpow[19]=(10**20)/(2**67) */ 181e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xd8d726b7177a8000), /* _Stl_tenpow[20]=(10**21)/(2**70) */ 182e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0x878678326eac9000), /* _Stl_tenpow[21]=(10**22)/(2**74) */ 183e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xa968163f0a57b400), /* _Stl_tenpow[22]=(10**23)/(2**77) */ 184e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xd3c21bcecceda100), /* _Stl_tenpow[23]=(10**24)/(2**80) */ 185e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0x84595161401484a0), /* _Stl_tenpow[24]=(10**25)/(2**84) */ 186e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xa56fa5b99019a5c8), /* _Stl_tenpow[25]=(10**26)/(2**87) */ 187e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xcecb8f27f4200f3a), /* _Stl_tenpow[26]=(10**27)/(2**90) */ 188e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 189e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xd0cf4b50cfe20766), /* _Stl_tenpow[27]=(10**55)/(2**183) */ 190e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xd2d80db02aabd62c), /* _Stl_tenpow[28]=(10**83)/(2**276) */ 191e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xd4e5e2cdc1d1ea96), /* _Stl_tenpow[29]=(10**111)/(2**369) */ 192e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xd6f8d7509292d603), /* _Stl_tenpow[30]=(10**139)/(2**462) */ 193e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xd910f7ff28069da4), /* _Stl_tenpow[31]=(10**167)/(2**555) */ 194e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xdb2e51bfe9d0696a), /* _Stl_tenpow[32]=(10**195)/(2**648) */ 195e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xdd50f1996b947519), /* _Stl_tenpow[33]=(10**223)/(2**741) */ 196e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xdf78e4b2bd342cf7), /* _Stl_tenpow[34]=(10**251)/(2**834) */ 197e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xe1a63853bbd26451), /* _Stl_tenpow[35]=(10**279)/(2**927) */ 198e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xe3d8f9e563a198e5), /* _Stl_tenpow[36]=(10**307)/(2**1020) */ 199e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 200e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// /* _Stl_tenpow[36]=(10**335)/(2**) */ 201e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// /* _Stl_tenpow[36]=(10**335)/(2**) */ 202e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 203e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xfd87b5f28300ca0e), /* _Stl_tenpow[37]=(10**-28)/(2**-93) */ 204e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xfb158592be068d2f), /* _Stl_tenpow[38]=(10**-56)/(2**-186) */ 205e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xf8a95fcf88747d94), /* _Stl_tenpow[39]=(10**-84)/(2**-279) */ 206e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xf64335bcf065d37d), /* _Stl_tenpow[40]=(10**-112)/(2**-372) */ 207e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xf3e2f893dec3f126), /* _Stl_tenpow[41]=(10**-140)/(2**-465) */ 208e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xf18899b1bc3f8ca2), /* _Stl_tenpow[42]=(10**-168)/(2**-558) */ 209e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xef340a98172aace5), /* _Stl_tenpow[43]=(10**-196)/(2**-651) */ 210e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xece53cec4a314ebe), /* _Stl_tenpow[44]=(10**-224)/(2**-744) */ 211e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xea9c227723ee8bcb), /* _Stl_tenpow[45]=(10**-252)/(2**-837) */ 212e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xe858ad248f5c22ca), /* _Stl_tenpow[46]=(10**-280)/(2**-930) */ 213e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xe61acf033d1a45df), /* _Stl_tenpow[47]=(10**-308)/(2**-1023) */ 214e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xe3e27a444d8d98b8), /* _Stl_tenpow[48]=(10**-336)/(2**-1116) */ 215e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottULL(0xe1afa13afbd14d6e) /* _Stl_tenpow[49]=(10**-364)/(2**-1209) */ 216e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott}; 217e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 218e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottstatic const short _Stl_twoexp[80] = { 219e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott4,7,10,14,17,20,24,27,30,34,37,40,44,47,50,54,57,60,64,67,70,74,77,80,84,87,90, 220e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott183,276,369,462,555,648,741,834,927,1020, 221e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott-93,-186,-279,-372,-465,-558,-651,-744,-837,-930,-1023,-1116,-1209 222e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott}; 223e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 224e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#define TEN_1 0 /* offset to 10 ** 1 */ 225e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#define TEN_27 26 /* offset to 10 ** 27 */ 226e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#define TEN_M28 37 /* offset to 10 ** -28 */ 227e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#define NUM_HI_P 11 228e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#define NUM_HI_N 13 229e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 230e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#define _Stl_HIBITULL (ULL(1) << 63) 231e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 232e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottstatic void _Stl_norm_and_round(uint64& p, int& norm, uint64 prodhi, uint64 prodlo) { 233e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott norm = 0; 234e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((prodhi & _Stl_HIBITULL) == 0) { 235e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* leading bit is a zero 236e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * may have to normalize 237e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 238e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((prodhi == ~_Stl_HIBITULL) && 239e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott ((prodlo >> 62) == 0x3)) { /* normalization followed by round 240e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * would cause carry to create 241e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * extra bit, so don't normalize 242e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 243e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott p = _Stl_HIBITULL; 244e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return; 245e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 246e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott p = (prodhi << 1) | (prodlo >> 63); /* normalize */ 247e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott norm = 1; 248e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott prodlo <<= 1; 249e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 250e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott else { 251e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott p = prodhi; 252e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 253e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 254e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((prodlo & _Stl_HIBITULL) != 0) { /* first guard bit a one */ 255e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (((p & 0x1) != 0) || 256e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott prodlo != _Stl_HIBITULL ) { /* not borderline for round to even */ 257e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* round */ 258e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott ++p; 259e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (p == 0) 260e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott ++p; 261e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 262e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 263e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 264e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 265e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// Convert a 64-bitb fraction * 10^exp to a 64-bit fraction * 2^bexp. 266e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// p: 64-bit fraction 267e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// exp: base-10 exponent 268e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// bexp: base-2 exponent (output parameter) 269e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottstatic void _Stl_tenscale(uint64& p, int exp, int& bexp) { 270e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bexp = 0; 271e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 272e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ( exp == 0 ) { /* no scaling needed */ 273e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return; 274e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 275e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 276e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int exp_hi = 0, exp_lo = exp; /* exp = exp_hi*32 + exp_lo */ 277e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int tlo = TEN_1, thi; /* offsets in power of ten table */ 278e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int num_hi; /* number of high exponent powers */ 279e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 280e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (exp > 0) { /* split exponent */ 281e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (exp_lo > 27) { 282e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp_lo++; 283e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott while (exp_lo > 27) { 284e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp_hi++; 285e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp_lo -= 28; 286e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 287e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 288e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott thi = TEN_27; 289e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott num_hi = NUM_HI_P; 290e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else { // exp < 0 291e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott while (exp_lo < 0) { 292e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp_hi++; 293e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp_lo += 28; 294e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 295e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott thi = TEN_M28; 296e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott num_hi = NUM_HI_N; 297e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 298e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 299e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64 prodhi, prodlo; /* 128b product */ 300e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int norm; /* number of bits of normalization */ 301e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 302e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int hi, lo; /* offsets in power of ten table */ 303e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott while (exp_hi) { /* scale */ 304e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott hi = (min) (exp_hi, num_hi); /* only a few large powers of 10 */ 305e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp_hi -= hi; /* could iterate in extreme case */ 306e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott hi += thi-1; 307e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _Stl_mult64(p, _Stl_tenpow[hi], prodhi, prodlo); 308e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _Stl_norm_and_round(p, norm, prodhi, prodlo); 309e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bexp += _Stl_twoexp[hi] - norm; 310e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 311e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 312e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (exp_lo) { 313e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott lo = tlo + exp_lo -1; 314e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _Stl_mult64(p, _Stl_tenpow[lo], prodhi, prodlo); 315e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _Stl_norm_and_round(p, norm, prodhi, prodlo); 316e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bexp += _Stl_twoexp[lo] - norm; 317e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 318e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 319e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return; 320e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 321e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 322e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// First argument is a buffer of values from 0 to 9, NOT ascii. 323e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// Second argument is number of digits in buffer, 1 <= digits <= 17. 324e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// Third argument is base-10 exponent. 325e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 326e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott/* IEEE representation */ 3277924895a54b8f0b67514224a7ad9d67c79b3eeafAndrew Hsieh#if !defined (__linux__) || defined (__ANDROID__) 328e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 329e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottunion _Double_rep { 330e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64 ival; 331e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott double val; 332e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott}; 333e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 334e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottstatic double _Stl_atod(char *buffer, ptrdiff_t ndigit, int dexp) { 335e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott typedef numeric_limits<double> limits; 336e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _Double_rep drep; 337e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64 &value = drep.ival; /* Value develops as follows: 338e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 1) decimal digits as an integer 339e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 2) left adjusted fraction 340e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 3) right adjusted fraction 341e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 4) exponent and fraction 342e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 343e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 344e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint32 guard; /* First guard bit */ 345e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64 rest; /* Remaining guard bits */ 346e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 347e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int bexp; /* binary exponent */ 348e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int nzero; /* number of non-zero bits */ 349e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int sexp; /* scaling exponent */ 350e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 351e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott char *bufferend; /* pointer to char after last digit */ 352e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 353e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Convert the decimal digits to a binary integer. */ 354e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bufferend = buffer + ndigit; 355e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value = 0; 356e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 357e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott while (buffer < bufferend) { 358e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value *= 10; 359e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value += *buffer++; 360e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 361e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 362e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Check for zero and treat it as a special case */ 363e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (value == 0) { 364e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return 0.0; 365e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 366e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 367e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Normalize value */ 368e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bexp = 64; /* convert from 64b int to fraction */ 369e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 370e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Count number of non-zeroes in value */ 371e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott nzero = 0; 372e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((value >> 32) != 0) { nzero = 32; } //*TY 03/25/2000 - added explicit comparison to zero to avoid uint64 to bool conversion operator 373e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((value >> (16 + nzero)) != 0) { nzero += 16; } 374e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((value >> ( 8 + nzero)) != 0) { nzero += 8; } 375e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((value >> ( 4 + nzero)) != 0) { nzero += 4; } 376e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((value >> ( 2 + nzero)) != 0) { nzero += 2; } 377e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((value >> ( 1 + nzero)) != 0) { nzero += 1; } 378e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((value >> ( nzero)) != 0) { nzero += 1; } 379e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 380e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Normalize */ 381e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value <<= /*(uint64)*/ (64 - nzero); //*TY 03/25/2000 - removed extraneous cast to uint64 382e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bexp -= 64 - nzero; 383e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 384e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* At this point we have a 64b fraction and a binary exponent 385e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * but have yet to incorporate the decimal exponent. 386e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 387e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 388e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* multiply by 10^dexp */ 389e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _Stl_tenscale(value, dexp, sexp); 390e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bexp += sexp; 391e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 392e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (bexp <= -1022) { /* HI denorm or underflow */ 393e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bexp += 1022; 394e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (bexp < -53) { /* guaranteed underflow */ 395e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value = 0; 396e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 397e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott else { /* denorm or possible underflow */ 398e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int lead0 = 12 - bexp; /* 12 sign and exponent bits */ 399e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 400e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* we must special case right shifts of more than 63 */ 401e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (lead0 > 64) { 402e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott rest = value; 403e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott guard = 0; 404e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value = 0; 405e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 406e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott else if (lead0 == 64) { 407e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott rest = value & ((ULL(1)<< 63)-1); 408e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott guard = (uint32) ((value>> 63) & 1 ); 409e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value = 0; 410e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 411e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott else { 412e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott rest = value & (((ULL(1) << lead0)-1)-1); 413e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott guard = (uint32) (((value>> lead0)-1) & 1); 414e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value >>= /*(uint64)*/ lead0; /* exponent is zero */ 415e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 416e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 417e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Round */ 418e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (guard && ((value & 1) || rest) ) { 419e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott ++value; 420e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (value == (ULL(1) << (limits::digits - 1))) { /* carry created normal number */ 421e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value = 0; 422e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _Stl_set_exponent(value, 1); 423e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 424e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 425e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 426e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 427e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott else { /* not zero or denorm */ 428e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Round to 53 bits */ 429e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott rest = value & ((1 << 10) - 1); 430e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value >>= 10; 431e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott guard = (uint32) value & 1; 432e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value >>= 1; 433e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 434e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* value&1 guard rest Action 435e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 436e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * dc 0 dc none 437e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 1 1 dc round 438e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 0 1 0 none 439e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 0 1 !=0 round 440e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 441e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (guard) { 442e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (((value&1)!=0) || (rest!=0)) { 443e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott ++value; /* round */ 444e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((value >> 53) != 0) { /* carry all the way across */ 445e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value >>= 1; /* renormalize */ 446e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott ++bexp; 447e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 448e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 449e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 450e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* 451e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Check for overflow 452e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * IEEE Double Precision Format 453e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * (From Table 7-8 of Kane and Heinrich) 454e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 455e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Fraction bits 52 456e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Emax +1023 457e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Emin -1022 458e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Exponent bias +1023 459e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Exponent bits 11 460e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Integer bit hidden 461e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Total width in bits 64 462e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 463e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 464e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (bexp > limits::max_exponent) { /* overflow */ 465e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return limits::infinity(); 466e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 467e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott else { /* value is normal */ 468e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott value &= ~(ULL(1) << (limits::digits - 1)); /* hide hidden bit */ 469e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _Stl_set_exponent(value, bexp + 1022); /* add bias */ 470e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 471e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 472e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 473e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _STLP_STATIC_ASSERT(sizeof(uint64) >= sizeof(double)) 474e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return drep.val; 475e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 476e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 477e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif 478e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 479e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#if defined (__linux__) || defined (__MINGW32__) || defined (__CYGWIN__) || \ 480e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott defined (__BORLANDC__) || defined (__DMC__) || defined (__HP_aCC) 481e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 482e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scotttemplate <class D, class IEEE, int M, int BIAS> 483e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottD _Stl_atodT(char *buffer, ptrdiff_t ndigit, int dexp) 484e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott{ 485e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott typedef numeric_limits<D> limits; 486e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 487e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Convert the decimal digits to a binary integer. */ 488e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott char *bufferend = buffer + ndigit; /* pointer to char after last digit */ 489e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _ll vv; 490e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 = 0L; 491e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 492e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott while ( buffer < bufferend ) { 493e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 *= 10; 494e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 += *buffer++; 495e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 496e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 497e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ( vv.i64 == ULL(0) ) { /* Check for zero and treat it as a special case */ 498e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return D(0.0); 499e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 500e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 501e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Normalize value */ 502e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 503e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int bexp = 64; /* convert from 64b int to fraction */ 504e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 505e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Count number of non-zeroes in value */ 506e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int nzero = 0; 507e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((vv.i64 >> 32) != 0) { nzero = 32; } 508e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((vv.i64 >> (16 + nzero)) != 0) { nzero += 16; } 509e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((vv.i64 >> ( 8 + nzero)) != 0) { nzero += 8; } 510e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((vv.i64 >> ( 4 + nzero)) != 0) { nzero += 4; } 511e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((vv.i64 >> ( 2 + nzero)) != 0) { nzero += 2; } 512e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((vv.i64 >> ( 1 + nzero)) != 0) { nzero += 1; } 513e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((vv.i64 >> ( nzero)) != 0) { nzero += 1; } 514e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 515e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Normalize */ 516e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott nzero = 64 - nzero; 517e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 <<= nzero; // * TY 03/25/2000 - removed extraneous cast to uint64 518e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bexp -= nzero; 519e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 520e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* At this point we have a 64b fraction and a binary exponent 521e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * but have yet to incorporate the decimal exponent. 522e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 523e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 524e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* multiply by 10^dexp */ 525e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int sexp; 526e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _Stl_tenscale(vv.i64, dexp, sexp); 527e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bexp += sexp; 528e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 529e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ( bexp >= limits::min_exponent ) { /* not zero or denorm */ 530e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ( limits::digits < 64 ) { 531e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Round to (64 - M + 1) bits */ 532e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64_t rest = vv.i64 & ((~ULL(0) / ULL(2)) >> (limits::digits - 1)); 533e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 >>= M - 2; 534e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint32_t guard = (uint32) vv.i64 & 1; 535e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 >>= 1; 536e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 537e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* value&1 guard rest Action 538e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 539e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * dc 0 dc none 540e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 1 1 dc round 541e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 0 1 0 none 542e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 0 1 !=0 round 543e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 544e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 545e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (guard) { 546e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ( ((vv.i64 & 1) != 0) || (rest != 0) ) { 547e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64++; /* round */ 548e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ( (vv.i64 >> (limits::digits < 64 ? limits::digits : 0)) != 0 ) { /* carry all the way across */ 549e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 >>= 1; /* renormalize */ 550e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott ++bexp; 551e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 552e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 553e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 554e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 555e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 &= ~(ULL(1) << (limits::digits - 1)); /* hide hidden bit */ 556e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 557e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* 558e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Check for overflow 559e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * IEEE Double Precision Format 560e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * (From Table 7-8 of Kane and Heinrich) 561e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * 562e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Fraction bits 52 563e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Emax +1023 564e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Emin -1022 565e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Exponent bias +1023 566e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Exponent bits 11 567e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Integer bit hidden 568e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Total width in bits 64 569e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 570e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 571e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (bexp > limits::max_exponent) { /* overflow */ 572e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return limits::infinity(); 573e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 574e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 575e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* value is normal */ 576e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 577e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott IEEE v; 578e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 579e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.mantissa0 = vv.i32.hi; 580e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.mantissa1 = vv.i32.lo; 581e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.negative = 0; 582e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.exponent = bexp + BIAS - 1; 583e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 584e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return v.d; 585e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 586e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 587e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* HI denorm or underflow */ 588e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bexp += BIAS - 1; 589e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (bexp < -limits::digits) { /* guaranteed underflow */ 590e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 = 0; 591e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else { /* denorm or possible underflow */ 592e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 593e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* 594e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * Problem point for long double: looks like this code reflect shareing of mantissa 595e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * and exponent in 64b int; not so for long double 596e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 597e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 598e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int lead0 = M - bexp; /* M = 12 sign and exponent bits */ 599e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint64_t rest; 600e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott uint32_t guard; 601e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 602e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* we must special case right shifts of more than 63 */ 603e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 604e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (lead0 > 64) { 605e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott rest = vv.i64; 606e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott guard = 0; 607e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 = 0; 608e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else if (lead0 == 64) { 609e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott rest = vv.i64 & ((ULL(1) << 63)-1); 610e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott guard = (uint32) ((vv.i64 >> 63) & 1 ); 611e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 = 0; 612e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else { 613e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott rest = vv.i64 & (((ULL(1) << lead0)-1)-1); 614e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott guard = (uint32) (((vv.i64 >> lead0)-1) & 1); 615e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64 >>= /*(uint64)*/ lead0; /* exponent is zero */ 616e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 617e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 618e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Round */ 619e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (guard && ( (vv.i64 & 1) || rest)) { 620e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott vv.i64++; 621e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (vv.i64 == (ULL(1) << (limits::digits - 1))) { /* carry created normal number */ 622e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott IEEE v; 623e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 624e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.mantissa0 = 0; 625e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.mantissa1 = 0; 626e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.negative = 0; 627e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.exponent = 1; 628e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return v.d; 629e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 630e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 631e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 632e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 633e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott IEEE v; 634e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 635e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.mantissa0 = vv.i32.hi; 636e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.mantissa1 = vv.i32.lo; 637e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.negative = 0; 638e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott v.ieee.exponent = 0; 639e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 640e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return v.d; 641e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 642e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif // __linux__ 643e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 6447924895a54b8f0b67514224a7ad9d67c79b3eeafAndrew Hsieh#if !defined (__linux__) || defined (__ANDROID__) 645e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottstatic double _Stl_string_to_double(const char *s) { 646e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott typedef numeric_limits<double> limits; 647e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott const int max_digits = limits::digits10 + 2; 648e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott unsigned c; 649e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott unsigned Negate, decimal_point; 650e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott char *d; 651e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int exp; 652e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int dpchar; 653e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott char digits[max_digits]; 654e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 655e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 656e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 657e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* process sign */ 658e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott Negate = 0; 659e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c == '+') { 660e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 661e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else if (c == '-') { 662e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott Negate = 1; 663e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 664e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 665e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 666e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott d = digits; 667e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott dpchar = '.' - '0'; 668e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott decimal_point = 0; 669e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp = 0; 670e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 671e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott for (;;) { 672e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c -= '0'; 673e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c < 10) { 674e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (d == digits + max_digits) { 675e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* ignore more than max_digits digits, but adjust exponent */ 676e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp += (decimal_point ^ 1); 677e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else { 678e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c == 0 && d == digits) { 679e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* ignore leading zeros */ 680e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else { 681e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott *d++ = (char) c; 682e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 683e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp -= decimal_point; 684e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 685e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else if (c == (unsigned int) dpchar && !decimal_point) { /* INTERNATIONAL */ 686e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott decimal_point = 1; 687e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else { 688e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott break; 689e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 690e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 691e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 692e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 693e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* strtod cant return until it finds the end of the exponent */ 694e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (d == digits) { 695e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return 0.0; 696e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 697e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 698e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c == 'e' - '0' || c == 'E' - '0') { 699e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott register unsigned negate_exp = 0; 700e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott register int e = 0; 701e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 702e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c == '+' || c == ' ') { 703e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 704e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else if (c == '-') { 705e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott negate_exp = 1; 706e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 707e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 708e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c -= '0', c < 10) { 709e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott do { 710e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott e = e * 10 + (int)c; 711e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 712e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } while (c -= '0', c < 10); 713e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 714e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (negate_exp) { 715e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott e = -e; 716e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 717e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp += e; 718e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 719e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 720e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 721e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott double x; 722e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott ptrdiff_t n = d - digits; 723e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((exp + n - 1) < limits::min_exponent10) { 724e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott x = 0; 725e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 726e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott else if ((exp + n - 1) > limits::max_exponent10) { 727e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott x = limits::infinity(); 728e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 729e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott else { 730e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* Let _Stl_atod diagnose under- and over-flows. 731e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * If the input was == 0.0, we have already returned, 732e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott * so retval of +-Inf signals OVERFLOW, 0.0 UNDERFLOW */ 733e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott x = _Stl_atod(digits, n, exp); 734e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 735e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 736e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (Negate) { 737e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott x = -x; 738e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 739e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 740e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return x; 741e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 742e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 743e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif 744e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 745e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#if defined (__linux__) || defined (__MINGW32__) || defined (__CYGWIN__) || \ 746e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott defined (__BORLANDC__) || defined (__DMC__) || defined (__HP_aCC) 747e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 748e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scotttemplate <class D, class IEEE, int M, int BIAS> 749e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick ScottD _Stl_string_to_doubleT(const char *s) 750e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott{ 751e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott typedef numeric_limits<D> limits; 752e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott const int max_digits = limits::digits10; /* + 2 17 */; 753e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott unsigned c; 754e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott unsigned decimal_point; 755e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott char *d; 756e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int exp; 757e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott D x; 758e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott int dpchar; 759e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott char digits[max_digits]; 760e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 761e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 762e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 763e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* process sign */ 764e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bool Negate = false; 765e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c == '+') { 766e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 767e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else if (c == '-') { 768e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott Negate = true; 769e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 770e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 771e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 772e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott d = digits; 773e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott dpchar = '.' - '0'; 774e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott decimal_point = 0; 775e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp = 0; 776e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 777e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott for (;;) { 778e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c -= '0'; 779e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c < 10) { 780e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (d == digits + max_digits) { 781e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* ignore more than max_digits digits, but adjust exponent */ 782e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp += (decimal_point ^ 1); 783e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else { 784e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c == 0 && d == digits) { 785e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* ignore leading zeros */ 786e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else { 787e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott *d++ = (char) c; 788e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 789e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp -= decimal_point; 790e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 791e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else if (c == (unsigned int) dpchar && !decimal_point) { /* INTERNATIONAL */ 792e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott decimal_point = 1; 793e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else { 794e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott break; 795e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 796e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 797e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 798e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* strtod cant return until it finds the end of the exponent */ 799e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (d == digits) { 800e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return D(0.0); 801e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 802e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 803e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c == 'e'-'0' || c == 'E'-'0') { 804e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott bool negate_exp = false; 805e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott register int e = 0; 806e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 807e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c == '+' || c == ' ') { 808e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 809e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else if (c == '-') { 810e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott negate_exp = true; 811e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 812e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 813e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (c -= '0', c < 10) { 814e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott do { 815e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott e = e * 10 + (int)c; 816e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott c = *s++; 817e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } while (c -= '0', c < 10); 818e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 819e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if (negate_exp) { 820e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott e = -e; 821e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 822e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott exp += e; 823e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 824e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 825e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 826e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott ptrdiff_t n = d - digits; 827e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott if ((exp + n - 1) < limits::min_exponent10) { 828e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return D(0.0); // +0.0 is the same as -0.0 829e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else if ((exp + n - 1) > limits::max_exponent10 ) { 830e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott // not good, because of x = -x below; this may lead to portability problems 831e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott x = limits::infinity(); 832e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } else { 833e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* let _Stl_atod diagnose under- and over-flows */ 834e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott /* if the input was == 0.0, we have already returned, 835e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott so retval of +-Inf signals OVERFLOW, 0.0 UNDERFLOW 836e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott */ 837e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott x = _Stl_atodT<D,IEEE,M,BIAS>(digits, n, exp); 838e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott } 839e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 840e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott return Negate ? -x : x; 841e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 842e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 843e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif // __linux__ 844e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 845e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottvoid _STLP_CALL 846e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott__string_to_float(const __iostring& v, float& val) 847e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott{ 8487924895a54b8f0b67514224a7ad9d67c79b3eeafAndrew Hsieh#if !defined (__linux__) || defined (__ANDROID__) 849e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott val = (float)_Stl_string_to_double(v.c_str()); 850e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#else 851e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott val = (float)_Stl_string_to_doubleT<double,ieee754_double,12,IEEE754_DOUBLE_BIAS>(v.c_str()); 852e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif 853e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 854e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 855e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottvoid _STLP_CALL 856e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott__string_to_float(const __iostring& v, double& val) 857e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott{ 8587924895a54b8f0b67514224a7ad9d67c79b3eeafAndrew Hsieh#if !defined (__linux__) || defined (__ANDROID__) 859e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott val = _Stl_string_to_double(v.c_str()); 860e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#else 861e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott val = _Stl_string_to_doubleT<double,ieee754_double,12,IEEE754_DOUBLE_BIAS>(v.c_str()); 862e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif 863e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 864e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 865e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#if !defined (_STLP_NO_LONG_DOUBLE) 866e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scottvoid _STLP_CALL 867e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott__string_to_float(const __iostring& v, long double& val) { 868e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#if !defined (__linux__) && !defined (__MINGW32__) && !defined (__CYGWIN__) && \ 869e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott !defined (__BORLANDC__) && !defined (__DMC__) && !defined (__HP_aCC) 870e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott //The following function is valid only if long double is an alias for double. 871e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott _STLP_STATIC_ASSERT( sizeof(long double) <= sizeof(double) ) 872e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott val = _Stl_string_to_double(v.c_str()); 873e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#else 874e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott val = _Stl_string_to_doubleT<long double,ieee854_long_double,16,IEEE854_LONG_DOUBLE_BIAS>(v.c_str()); 875e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif 876e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott} 877e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott#endif 878e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 879e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott_STLP_MOVE_TO_STD_NAMESPACE 880e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott_STLP_END_NAMESPACE 881e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott 882e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// Local Variables: 883e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// mode:C++ 884e46c9386c4f79aa40185f79a19fc5b2a7ef528b3Patrick Scott// End: 885