diff -r 000000000000 -r 6474c204b198 build/stlport/src/num_put_float.cpp --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/build/stlport/src/num_put_float.cpp Wed Dec 31 06:09:35 2014 +0100 @@ -0,0 +1,919 @@ +/* + * Copyright (c) 1999 + * Silicon Graphics Computer Systems, Inc. + * + * Copyright (c) 1999 + * Boris Fomitchev + * + * This material is provided "as is", with absolutely no warranty expressed + * or implied. Any use is at your own risk. + * + * Permission to use or copy this software for any purpose is hereby granted + * without fee, provided the above notices are retained on all copies. + * Permission to modify the code and to distribute modified code is granted, + * provided the above notices are retained, and a notice that the code was + * modified is included with the above copyright notice. + * + */ + +#include "stlport_prefix.h" + +#include +#include +#include + +#if defined (__DECCXX) +# define NDIG 400 +#else +# define NDIG 82 +#endif + +#define todigit(x) ((x)+'0') + +#if defined (_STLP_UNIX) + +# if defined (__sun) +# include +# endif + +# if defined (__sun) || defined (__digital__) || defined (__sgi) || defined (_STLP_SCO_OPENSERVER) || defined (__NCR_SVR) +// DEC, SGI & Solaris need this +# include +# include +# endif + +# if defined (__QNXNTO__) || ( defined(__GNUC__) && defined(__APPLE__) ) || defined(_STLP_USE_UCLIBC) /* 0.9.26 */ || \ + defined(__FreeBSD__) +# define USE_SPRINTF_INSTEAD +# endif + +# if defined (_AIX) // JFA 3-Aug-2000 +# include +# include +# endif + +# include +#endif + +#include +#include + +#if defined (_STLP_MSVC_LIB) || defined (__MINGW32__) || defined (__BORLANDC__) || defined (__DJGPP) || \ + defined (_STLP_SCO_OPENSERVER) || defined (__NCR_SVR) +# include +#endif + +#if defined (__MRC__) || defined (__SC__) || defined (_CRAY) //*TY 02/24/2000 - added support for MPW +# include +#endif + +#if defined (__CYGWIN__) +# include +#endif + +#if defined (__MSL__) +# include // for atoi +# include // for snprintf +# include +# include +#endif + +#if defined (__ISCPP__) +# include +#endif + +#include + +#if defined (__DMC__) +# define snprintf _snprintf +#endif + +_STLP_BEGIN_NAMESPACE + +_STLP_MOVE_TO_PRIV_NAMESPACE + +#if defined (__MWERKS__) || defined(__BEOS__) +# define USE_SPRINTF_INSTEAD +#endif + +template +struct _Dig +{ + enum { dig = _Dig::dig + 1 }; +}; + +_STLP_TEMPLATE_NULL +struct _Dig<0> +{ + enum { dig = 0 }; +}; + +#ifdef _STLP_NO_LONG_DOUBLE +# define MAXEDIGITS int(_Dig::dig) +# define MAXFSIG DBL_DIG +# define MAXFCVT (DBL_DIG + 1) +#else +# define MAXEDIGITS int(_Dig::dig) +# define MAXFSIG LDBL_DIG +# define MAXFCVT (LDBL_DIG + 1) +#endif + +// Tests for infinity and NaN differ on different OSs. We encapsulate +// these differences here. +#if !defined (USE_SPRINTF_INSTEAD) +# if defined (__hpux) && defined (__GNUC__) +# define _STLP_USE_SIGN_HELPER +# elif defined (__DJGPP) || (defined (_STLP_USE_GLIBC) && ! defined (__MSL__)) || \ + defined (__CYGWIN__) || \ + defined (__FreeBSD__) || defined (__NetBSD__) || defined (__OpenBSD__) || \ + defined (__HP_aCC) +static inline bool _Stl_is_nan_or_inf(double x) +# if defined (isfinite) +{ return !isfinite(x); } +# else +{ return !finite(x); } +# endif +static inline bool _Stl_is_neg_nan(double x) { return isnan(x) && ( copysign(1., x) < 0 ); } +static inline bool _Stl_is_inf(double x) { return isinf(x); } +// inline bool _Stl_is_neg_inf(double x) { return isinf(x) < 0; } +static inline bool _Stl_is_neg_inf(double x) { return isinf(x) && x < 0; } +# elif (defined (__unix) || defined (__unix__)) && \ + !defined (__APPLE__) && !defined (__DJGPP) && !defined(__osf__) && \ + !defined (_CRAY) && !defined (__ANDROID__) +static inline bool _Stl_is_nan_or_inf(double x) { return IsNANorINF(x); } +static inline bool _Stl_is_inf(double x) { return IsNANorINF(x) && IsINF(x); } +static inline bool _Stl_is_neg_inf(double x) { return (IsINF(x)) && (x < 0.0); } +static inline bool _Stl_is_neg_nan(double x) { return IsNegNAN(x); } +# elif defined (_STLP_MSVC_LIB) || defined (__MINGW32__) || defined (__BORLANDC__) +static inline bool _Stl_is_nan_or_inf(double x) { return !_finite(x); } +# if !defined (__BORLANDC__) +static inline bool _Stl_is_inf(double x) { + int fclass = _fpclass(x); + return fclass == _FPCLASS_NINF || fclass == _FPCLASS_PINF; +} +static inline bool _Stl_is_neg_inf(double x) { return _fpclass(x) == _FPCLASS_NINF; } +# else +static inline bool _Stl_is_inf(double x) { return _Stl_is_nan_or_inf(x) && !_isnan(x);} +static inline bool _Stl_is_neg_inf(double x) { return _Stl_is_inf(x) && x < 0 ; } +# endif +static inline bool _Stl_is_neg_nan(double x) { return _isnan(x) && _copysign(1., x) < 0 ; } +# if defined (__BORLANDC__) +static inline bool _Stl_is_nan_or_inf(long double x) { return !_finitel(x); } +static inline bool _Stl_is_inf(long double x) { return _Stl_is_nan_or_inf(x) && !_isnanl(x);} +static inline bool _Stl_is_neg_inf(long double x) { return _Stl_is_inf(x) && x < 0 ; } +static inline bool _Stl_is_neg_nan(long double x) { return _isnanl(x) && _copysignl(1.l, x) < 0 ; } +# elif !defined (_STLP_NO_LONG_DOUBLE) +// Simply there to avoid warning long double -> double implicit conversion: +static inline bool _Stl_is_nan_or_inf(long double x) { return _Stl_is_nan_or_inf(__STATIC_CAST(double, x)); } +static inline bool _Stl_is_inf(long double x) { return _Stl_is_inf(__STATIC_CAST(double, x));} +static inline bool _Stl_is_neg_inf(long double x) { return _Stl_is_neg_inf(__STATIC_CAST(double, x)); } +static inline bool _Stl_is_neg_nan(long double x) { return _Stl_is_neg_nan(__STATIC_CAST(double, x)); } +# endif +# elif defined (__MRC__) || defined (__SC__) || defined (__DMC__) +static bool _Stl_is_nan_or_inf(double x) { return isnan(x) || !isfinite(x); } +static bool _Stl_is_inf(double x) { return !isfinite(x); } +static bool _Stl_is_neg_inf(double x) { return !isfinite(x) && signbit(x); } +static bool _Stl_is_neg_nan(double x) { return isnan(x) && signbit(x); } +# elif /* defined(__FreeBSD__) || defined(__OpenBSD__) || */ (defined(__GNUC__) && defined(__APPLE__)) +static inline bool _Stl_is_nan_or_inf(double x) { return !finite(x); } +static inline bool _Stl_is_inf(double x) { return _Stl_is_nan_or_inf(x) && ! isnan(x); } +static inline bool _Stl_is_neg_inf(double x) { return _Stl_is_inf(x) && x < 0 ; } +static inline bool _Stl_is_neg_nan(double x) { return isnan(x) && copysign(1., x) < 0 ; } +# elif defined( _AIX ) // JFA 11-Aug-2000 +static bool _Stl_is_nan_or_inf(double x) { return isnan(x) || !finite(x); } +static bool _Stl_is_inf(double x) { return !finite(x); } +// bool _Stl_is_neg_inf(double x) { return _class(x) == FP_MINUS_INF; } +static bool _Stl_is_neg_inf(double x) { return _Stl_is_inf(x) && ( copysign(1., x) < 0 ); } +static bool _Stl_is_neg_nan(double x) { return isnan(x) && ( copysign(1., x) < 0 ); } +# elif defined (__ISCPP__) +static inline bool _Stl_is_nan_or_inf (double x) { return _fp_isINF(x) || _fp_isNAN(x); } +static inline bool _Stl_is_inf (double x) { return _fp_isINF(x); } +static inline bool _Stl_is_neg_inf (double x) { return _fp_isINF(x) && x < 0; } +static inline bool _Stl_is_neg_nan (double x) { return _fp_isNAN(x) && x < 0; } +# elif defined (_CRAY) +# if defined (_CRAYIEEE) +static inline bool _Stl_is_nan_or_inf(double x) { return isnan(x) || isinf(x); } +static inline bool _Stl_is_inf(double x) { return isinf(x); } +static inline bool _Stl_is_neg_inf(double x) { return isinf(x) && signbit(x); } +static inline bool _Stl_is_neg_nan(double x) { return isnan(x) && signbit(x); } +# else +static inline bool _Stl_is_nan_or_inf(double x) { return false; } +static inline bool _Stl_is_inf(double x) { return false; } +static inline bool _Stl_is_neg_inf(double x) { return false; } +static inline bool _Stl_is_neg_nan(double x) { return false; } +# endif +# else // nothing from above +# define USE_SPRINTF_INSTEAD +# endif +#endif // !USE_SPRINTF_INSTEAD + +#if !defined (USE_SPRINTF_INSTEAD) +// Reentrant versions of floating-point conversion functions. The argument +// lists look slightly different on different operating systems, so we're +// encapsulating the differences here. + +# if defined (__CYGWIN__) || defined(__DJGPP) +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf) +{ return ecvtbuf(x, n, pt, sign, buf); } +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf) +{ return fcvtbuf(x, n, pt, sign, buf); } +# if !defined (_STLP_NO_LONG_DOUBLE) +# if defined (__CYGWIN__) +# define _STLP_EMULATE_LONG_DOUBLE_CVT +# else +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf) +{ return ecvtbuf(x, n, pt, sign, buf); } +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf) +{ return fcvtbuf(x, n, pt, sign, buf); } +# endif +# endif +# elif defined (_STLP_USE_GLIBC) +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf, size_t bsize) +{ return ecvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0; } +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf, size_t bsize) +{ return fcvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0; } +# ifndef _STLP_NO_LONG_DOUBLE +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf, size_t bsize) +{ return qecvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0; } +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf, size_t bsize) +{ return qfcvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0; } +# endif +# define _STLP_NEED_CVT_BUFFER_SIZE +# elif defined (__sun) +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf) +{ return econvert(x, n, pt, sign, buf); } +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf) +{ return fconvert(x, n, pt, sign, buf); } +# ifndef _STLP_NO_LONG_DOUBLE +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf) +{ return qeconvert(&x, n, pt, sign, buf); } +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf) +{ return qfconvert(&x, n, pt, sign, buf); } +# endif +# elif defined (__DECCXX) +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf, size_t bsize) +{ return (ecvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0); } +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf, size_t bsize) +{ return (fcvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0); } +# ifndef _STLP_NO_LONG_DOUBLE +// fbp : no "long double" conversions ! +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf, size_t bsize) +{ return (ecvt_r((double)x, n, pt, sign, buf, bsize) == 0 ? buf : 0) ; } +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf, size_t bsize) +{ return (fcvt_r((double)x, n, pt, sign, buf, bsize) == 0 ? buf : 0); } +# endif +# define _STLP_NEED_CVT_BUFFER_SIZE +# elif defined (__hpux) +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign) +{ return ecvt(x, n, pt, sign); } +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign) +{ return fcvt(x, n, pt, sign); } +# if !defined (_STLP_NO_LONG_DOUBLE) +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign) +{ return _ldecvt(*(long_double*)&x, n, pt, sign); } +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign) +{ return _ldfcvt(*(long_double*)&x, n, pt, sign); } +# endif +# define _STLP_CVT_NEED_SYNCHRONIZATION +# elif defined (__unix) && !defined (__APPLE__) && !defined (_CRAY) && \ + !defined (__ANDROID__) +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf) +{ return ecvt_r(x, n, pt, sign, buf); } +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf) +{ return fcvt_r(x, n, pt, sign, buf); } +# if !defined (_STLP_NO_LONG_DOUBLE) +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf) +{ return qecvt_r(x, n, pt, sign, buf); } +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf) +{ return qfcvt_r(x, n, pt, sign, buf); } +# endif +# elif defined (_STLP_MSVC_LIB) || defined (__MINGW32__) || defined (__BORLANDC__) +# if defined (_STLP_USE_SAFE_STRING_FUNCTIONS) +# define _STLP_APPEND(a, b) a##b +# define _STLP_BUF_PARAMS , char* buf, size_t bsize +# define _STLP_SECURE_FUN(F, X, N, PT, SIGN) _STLP_APPEND(F, _s)(buf, bsize, X, N, PT, SIGN); return buf +# else +# define _STLP_BUF_PARAMS +# define _STLP_SECURE_FUN(F, X, N, PT, SIGN) return F(X, N, PT, SIGN) +# define _STLP_CVT_NEED_SYNCHRONIZATION +# endif +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign _STLP_BUF_PARAMS) +{ _STLP_SECURE_FUN(_ecvt, x, n, pt, sign); } +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign _STLP_BUF_PARAMS) +{ _STLP_SECURE_FUN(_fcvt, x, n, pt, sign); } +# if !defined (_STLP_NO_LONG_DOUBLE) +# if defined (_STLP_USE_SAFE_STRING_FUNCTIONS) +# define _STLP_PARAMS , buf, bsize +# else +# define _STLP_PARAMS +# endif +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign _STLP_BUF_PARAMS) +{ return _Stl_ecvtR(__STATIC_CAST(double, x), n, pt, sign _STLP_PARAMS); } +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign _STLP_BUF_PARAMS) +{ return _Stl_fcvtR(__STATIC_CAST(double, x), n, pt, sign _STLP_PARAMS); } +# undef _STLP_PARAMS +# endif +# undef _STLP_SECURE_FUN +# undef _STLP_BUF_PARAMS +# undef _STLP_APPEND +# if defined (__BORLANDC__) /* || defined (__GNUC__) MinGW do not support 'L' modifier so emulation do not work */ +# define _STLP_EMULATE_LONG_DOUBLE_CVT +# endif +# elif defined (__ISCPP__) +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf) +{ return _fp_ecvt( x, n, pt, sign, buf); } +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf) +{ return _fp_fcvt(x, n, pt, sign, buf); } +# if !defined (_STLP_NO_LONG_DOUBLE) +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf) +{ return _fp_ecvt( x, n, pt, sign, buf); } +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf) +{ return _fp_fcvt(x, n, pt, sign, buf); } +# endif +# elif defined (_AIX) || defined (__FreeBSD__) || defined (__NetBSD__) || defined (__OpenBSD__) || \ + defined (__MRC__) || defined (__SC__) || defined (_CRAY) || \ + defined (_STLP_SCO_OPENSERVER) || defined (__NCR_SVR) || \ + defined (__DMC__) +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign) +{ return ecvt(x, n, pt, sign ); } +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign) +{ return fcvt(x, n, pt, sign); } +# if !defined (_STLP_NO_LONG_DOUBLE) +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign) +{ return ecvt(x, n, pt, sign ); } +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign) +{ return fcvt(x, n, pt, sign); } +# endif +# define _STLP_CVT_NEED_SYNCHRONIZATION +# else +# error Missing _Stl_ecvtR and _Stl_fcvtR implementations. +# endif + +#if defined (_STLP_CVT_NEED_SYNCHRONIZATION) +/* STLport synchronize access to *cvt functions but those methods might + * be called from outside, in this case we will still have a race condition. */ +# if defined (_STLP_THREADS) +static _STLP_STATIC_MUTEX& put_float_mutex() { + static _STLP_STATIC_MUTEX __put_float_mutex _STLP_MUTEX_INITIALIZER; + return __put_float_mutex; +} +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf) { + _STLP_auto_lock lock(put_float_mutex()); + strcpy(buf, _Stl_ecvtR(x, n, pt, sign)); return buf; +} +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf) { + _STLP_auto_lock lock(put_float_mutex()); + strcpy(buf, _Stl_fcvtR(x, n, pt, sign)); return buf; +} +# if !defined (_STLP_NO_LONG_DOUBLE) && !defined (_STLP_EMULATE_LONG_DOUBLE_CVT) +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf) { + _STLP_auto_lock lock(put_float_mutex()); + strcpy(buf, _Stl_ecvtR(x, n, pt, sign)); return buf; +} +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf) { + _STLP_auto_lock lock(put_float_mutex()); + strcpy(buf, _Stl_fcvtR(x, n, pt, sign)); return buf; +} +# endif +# else +static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char*) +{ return _Stl_ecvtR(x, n, pt, sign); } +static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char*) +{ return _Stl_fcvtR(x, n, pt, sign); } +# if !defined (_STLP_NO_LONG_DOUBLE) && !defined (_STLP_EMULATE_LONG_DOUBLE_CVT) +static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char*) +{ return _Stl_ecvtR(x, n, pt, sign); } +static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char*) +{ return _Stl_fcvtR(x, n, pt, sign); } +# endif +# endif +#endif + +# if !defined (_STLP_USE_SAFE_STRING_FUNCTIONS) && !defined (_STLP_NEED_CVT_BUFFER_SIZE) +# define _STLP_CVT_BUFFER(B) B +# else +# define _STLP_CVT_BUFFER(B) _STLP_ARRAY_AND_SIZE(B) +# endif + +# if defined (_STLP_EMULATE_LONG_DOUBLE_CVT) +static void __fill_fmtbuf(char* fmtbuf, ios_base::fmtflags flags, char long_modifier); + +// Emulation of ecvt/fcvt functions using sprintf: +static char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf) { + // If long double value can be safely converted to double without losing precision + // we use the ecvt function for double: + double y = __STATIC_CAST(double, x); + if (x == y) + return _Stl_ecvtR(y, n, pt, sign, buf); + + char fmtbuf[32]; + __fill_fmtbuf(fmtbuf, 0, 'L'); + sprintf(buf, fmtbuf, n, x < 0.0l ? -x : x); + /* We are waiting for something having the form x.xxxe+yyyy */ + *pt = 0; + *sign = 0; + int i = -1; + int offset = 0; + while (buf[++i] != 0 && n != 0) { + if (buf[i] >= '0' && buf[i] <= '9') { + --n; + if (offset != 0) + buf[i - offset] = buf[i]; + } + else { + if (offset != 0) break; + ++offset; + *pt = i; + } + } + if (offset != 0) + buf[i - offset] = 0; + // Extract exponent part in point position: + int e = 0; + while (buf[++i] != 0) { + if (buf[i] >= '0' && buf[i] <= '9') { + e = e * 10 + (buf[i] - '0'); + } + } + *pt += e; + return buf; +} + +static char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf) { + // If long double value can be safely converted to double without losing precision + // we use the fcvt function for double: + double y = __STATIC_CAST(double, x); + if (x == y) + return _Stl_fcvtR(y, n, pt, sign, buf); + + char fmtbuf[32]; + __fill_fmtbuf(fmtbuf, ios_base::fixed, 'L'); + sprintf(buf, fmtbuf, n, x < 0.0l ? -x : x); + *pt = 0; + *sign = 0; + int i = -1; + int offset = 0; + while (buf[++i] != 0 && (offset == 0 || n != 0)) { + if (buf[i] >= '0' && buf[i] <= '9') { + if (offset != 0) { + --n; + buf[i - offset] = buf[i]; + } + } + else { + ++offset; + *pt = i; + } + } + if (offset != 0) + buf[i - offset] = 0; + else + *pt = i; + return buf; +} +#endif + +//---------------------------------------------------------------------- +// num_put + +// __format_float formats a mantissa and exponent as returned by +// one of the conversion functions (ecvt_r, fcvt_r, qecvt_r, qfcvt_r) +// according to the specified precision and format flags. This is +// based on doprnt but is much simpler since it is concerned only +// with floating point input and does not consider all formats. It +// also does not deal with blank padding, which is handled by +// __copy_float_and_fill. + +static size_t __format_float_scientific( __iostring& buf, const char *bp, + int decpt, int sign, bool is_zero, + ios_base::fmtflags flags, + int precision) { + // sign if required + if (sign) + buf += '-'; + else if (flags & ios_base::showpos) + buf += '+'; + + // first digit of mantissa + buf += *bp++; + + // start of grouping position, grouping won't occur in scientific notation + // as it is impossible to have something like 1234.0e04 but we return a correct + // group position for coherency with __format_float_fixed. + size_t __group_pos = buf.size(); + + // decimal point if required + if (precision != 0 || flags & ios_base::showpoint) { + buf += '.'; + } + + // rest of mantissa + while (*bp != 0 && precision--) + buf += *bp++; + + // trailing 0 if needed + if (precision > 0) + buf.append(precision, '0'); + + // exponent size = number of digits + exponent sign + exponent symbol + trailing zero + char expbuf[MAXEDIGITS + 3]; + //We start filling at the buffer end + char *suffix = expbuf + MAXEDIGITS + 2; + *suffix = 0; + if (!is_zero) { + int nn = decpt - 1; + if (nn < 0) + nn = -nn; + for (; nn > 9; nn /= 10) + *--suffix = (char) todigit(nn % 10); + *--suffix = (char) todigit(nn); + } + + // prepend leading zeros to exponent + // C89 Standard says that it should be at least 2 digits, C99 Standard says that + // we stop prepend zeros if more than 3 digits. To repect both STLport prepend zeros + // until it is 2 digits. + while (suffix > &expbuf[MAXEDIGITS]) + *--suffix = '0'; + + // put in the exponent sign + *--suffix = (char) ((decpt > 0 || is_zero ) ? '+' : '-'); + + // put in the e + *--suffix = flags & ios_base::uppercase ? 'E' : 'e'; + + // copy the suffix + buf += suffix; + return __group_pos; +} + +static size_t __format_float_fixed( __iostring &buf, const char *bp, + int decpt, int sign, + ios_base::fmtflags flags, + int precision) { + if ( sign && (decpt > -precision) && (*bp != 0) ) + buf += '-'; + else if ( flags & ios_base::showpos ) + buf += '+'; + + // digits before decimal point + int nnn = decpt; + do { + buf += (nnn <= 0 || *bp == 0) ? '0' : *bp++; + } while ( --nnn > 0 ); + + // start of grouping position + size_t __group_pos = buf.size(); + + // decimal point if needed + if ( flags & ios_base::showpoint || precision > 0 ) { + buf += '.'; + } + + // digits after decimal point if any + while ( *bp != 0 && --precision >= 0 ) { + buf += (++decpt <= 0) ? '0' : *bp++; + } + + // trailing zeros if needed + if (precision > 0) + buf.append(precision, '0'); + + return __group_pos; +} + +#if defined (_STLP_USE_SIGN_HELPER) +template +struct float_sign_helper { + float_sign_helper(_FloatT __x) + { _M_number._num = __x; } + + bool is_negative() const { + const unsigned short sign_mask(1 << (sizeof(unsigned short) * CHAR_BIT - 1)); + return (get_sign_word() & sign_mask) != 0; + } +private: + union { + unsigned short _Words[8]; + _FloatT _num; + } _M_number; + + unsigned short get_word_higher() const _STLP_NOTHROW + { return _M_number._Words[0]; } + unsigned short get_word_lower() const _STLP_NOTHROW + { return _M_number._Words[(sizeof(_FloatT) >= 12 ? 10 : sizeof(_FloatT)) / sizeof(unsigned short) - 1]; } + unsigned short get_sign_word() const _STLP_NOTHROW +# if defined (_STLP_BIG_ENDIAN) + { return get_word_higher(); } +# else /* _STLP_LITTLE_ENDIAN */ + { return get_word_lower(); } +# endif +}; +#endif + +template +static size_t __format_nan_or_inf(__iostring& buf, _FloatT x, ios_base::fmtflags flags) { + static const char* inf[2] = { "inf", "Inf" }; + static const char* nan[2] = { "nan", "NaN" }; + const char** inf_or_nan; +#if !defined (_STLP_USE_SIGN_HELPER) + if (_Stl_is_inf(x)) { // Infinity + inf_or_nan = inf; + if (_Stl_is_neg_inf(x)) + buf += '-'; + else if (flags & ios_base::showpos) + buf += '+'; + } else { // NaN + inf_or_nan = nan; + if (_Stl_is_neg_nan(x)) + buf += '-'; + else if (flags & ios_base::showpos) + buf += '+'; + } +#else + typedef numeric_limits<_FloatT> limits; + if (x == limits::infinity() || x == -limits::infinity()) { + inf_or_nan = inf; + } else { // NaN + inf_or_nan = nan; + } + float_sign_helper<_FloatT> helper(x); + if (helper.is_negative()) + buf += '-'; + else if (flags & ios_base::showpos) + buf += '+'; +#endif + size_t ret = buf.size(); + buf += inf_or_nan[flags & ios_base::uppercase ? 1 : 0]; + return ret; +} + +static inline size_t __format_float(__iostring &buf, const char * bp, + int decpt, int sign, bool is_zero, + ios_base::fmtflags flags, + int precision) { + size_t __group_pos = 0; + switch (flags & ios_base::floatfield) { + case ios_base::scientific: + __group_pos = __format_float_scientific( buf, bp, decpt, sign, is_zero, + flags, precision); + break; + case ios_base::fixed: + __group_pos = __format_float_fixed( buf, bp, decpt, sign, + flags, precision); + break; + default: // g format + // establish default precision + if (flags & ios_base::showpoint || precision > 0) { + if (precision == 0) precision = 1; + } else + precision = 6; + + // reset exponent if value is zero + if (is_zero) + decpt = 1; + + int kk = precision; + if (!(flags & ios_base::showpoint)) { + size_t n = strlen(bp); + if (n < (size_t)kk) + kk = (int)n; + while (kk >= 1 && bp[kk-1] == '0') + --kk; + } + + if (decpt < -3 || decpt > precision) { + precision = kk - 1; + __group_pos = __format_float_scientific( buf, bp, decpt, sign, is_zero, + flags, precision); + } else { + precision = kk - decpt; + __group_pos = __format_float_fixed( buf, bp, decpt, sign, + flags, precision); + } + break; + } /* switch */ + return __group_pos; +} + +#endif + +#if defined (USE_SPRINTF_INSTEAD) || defined (_STLP_EMULATE_LONG_DOUBLE_CVT) +struct GroupPos { + bool operator () (char __c) const { + return __c == '.' || + __c == 'e' || __c == 'E'; + } +}; + +// Creates a format string for sprintf() +static void __fill_fmtbuf(char* fmtbuf, ios_base::fmtflags flags, char long_modifier) { + fmtbuf[0] = '%'; + int i = 1; + + if (flags & ios_base::showpos) + fmtbuf[i++] = '+'; + + if (flags & ios_base::showpoint) + fmtbuf[i++] = '#'; + + fmtbuf[i++] = '.'; + fmtbuf[i++] = '*'; + + if (long_modifier) + fmtbuf[i++] = long_modifier; + + switch (flags & ios_base::floatfield) + { + case ios_base::scientific: + fmtbuf[i++] = (flags & ios_base::uppercase) ? 'E' : 'e'; + break; + case ios_base::fixed: +# if defined (__FreeBSD__) + fmtbuf[i++] = 'f'; +# else + fmtbuf[i++] = (flags & ios_base::uppercase) ? 'F' : 'f'; +# endif + break; + default: + fmtbuf[i++] = (flags & ios_base::uppercase) ? 'G' : 'g'; + break; + } + + fmtbuf[i] = 0; +} + +#endif /* USE_SPRINTF_INSTEAD */ + +template +static size_t __write_floatT(__iostring &buf, ios_base::fmtflags flags, int precision, + _FloatT x +#if defined (USE_SPRINTF_INSTEAD) + , char modifier) { + /* In theory, if we want 'arbitrary' precision, we should use 'arbitrary' + * buffer size below, but really we limited by exponent part in double. + * - ptr + */ + typedef numeric_limits<_FloatT> limits; + char static_buf[limits::max_exponent10 + 6]; // 6: -xxx.yyyE-zzz (sign, dot, E, exp sign, \0) + char fmtbuf[32]; + __fill_fmtbuf(fmtbuf, flags, modifier); + snprintf(_STLP_ARRAY_AND_SIZE(static_buf), fmtbuf, precision, x); + buf = static_buf; + return find_if(buf.begin(), buf.end(), GroupPos()) - buf.begin(); +#else + ) { + typedef numeric_limits<_FloatT> limits; + //If numeric_limits support is correct we use the exposed values to detect NaN and infinity: + if (limits::has_infinity && limits::has_quiet_NaN) { + if (!(x == x) || // NaN check + (x == limits::infinity() || x == -limits::infinity())) { + return __format_nan_or_inf(buf, x, flags); + } + } + // numeric_limits support is not good enough, we rely on platform dependent function + // _Stl_is_nan_or_inf that do not support long double. + else if (_Stl_is_nan_or_inf(x)) { + return __format_nan_or_inf(buf, x, flags); + } +# if defined (__MINGW32__) + //For the moment MinGW is limited to display at most numeric_limits::max() + if (x > numeric_limits::max() || + x < -numeric_limits::max()) { + return __format_nan_or_inf(buf, x, flags); + } +# endif + + /* Buffer size is max number of digits which is the addition of: + * - max_exponent10: max number of digits in fixed mode + * - digits10 + 2: max number of significant digits + * - trailing '\0' + */ + char cvtbuf[limits::max_exponent10 + limits::digits10 + 2 + 1]; + char *bp; + int decpt, sign; + + switch (flags & ios_base::floatfield) { + case ios_base::fixed: + { + /* Here, number of digits represents digits _after_ decimal point. + * In order to limit static buffer size we have to give 2 different values depending on x value. + * For small values (abs(x) < 1) we need as many digits as requested by precision limited by the maximum number of digits + * which is min_exponent10 + digits10 + 2 + * For bigger values we won't have more than limits::digits10 + 2 digits after decimal point. */ + int digits10 = (x > -1.0 && x < 1.0 ? -limits::min_exponent10 + limits::digits10 + 2 + : limits::digits10 + 2); + bp = _Stl_fcvtR(x, (min) (precision, digits10), &decpt, &sign, _STLP_CVT_BUFFER(cvtbuf) ); + } + break; + case ios_base::scientific: + default: + /* Here, number of digits is total number of digits which is limited to digits10 + 2. */ + { + int digits10 = limits::digits10 + 2; + bp = _Stl_ecvtR(x, (min) (precision, digits10), &decpt, &sign, _STLP_CVT_BUFFER(cvtbuf) ); + } + break; + } + return __format_float(buf, bp, decpt, sign, x == 0.0, flags, precision); +#endif +} + +size_t _STLP_CALL +__write_float(__iostring &buf, ios_base::fmtflags flags, int precision, + double x) { + return __write_floatT(buf, flags, precision, x +#if defined (USE_SPRINTF_INSTEAD) + , 0 +#endif + ); +} + +#if !defined (_STLP_NO_LONG_DOUBLE) +size_t _STLP_CALL +__write_float(__iostring &buf, ios_base::fmtflags flags, int precision, + long double x) { + return __write_floatT(buf, flags, precision, x +#if defined (USE_SPRINTF_INSTEAD) + , 'L' +#endif + ); +} +#endif + +void _STLP_CALL __get_floor_digits(__iostring &out, _STLP_LONGEST_FLOAT_TYPE __x) { + typedef numeric_limits<_STLP_LONGEST_FLOAT_TYPE> limits; +#if defined (USE_SPRINTF_INSTEAD) + char cvtbuf[limits::max_exponent10 + 6]; +# if !defined (_STLP_NO_LONG_DOUBLE) + snprintf(_STLP_ARRAY_AND_SIZE(cvtbuf), "%Lf", __x); // check for 1234.56! +# else + snprintf(_STLP_ARRAY_AND_SIZE(cvtbuf), "%f", __x); // check for 1234.56! +# endif + char *p = strchr( cvtbuf, '.' ); + if ( p == 0 ) { + out.append( cvtbuf ); + } else { + out.append( cvtbuf, p ); + } +#else + char cvtbuf[limits::max_exponent10 + 1]; + char * bp; + int decpt, sign; + bp = _Stl_fcvtR(__x, 0, &decpt, &sign, _STLP_CVT_BUFFER(cvtbuf)); + + if (sign) { + out += '-'; + } + out.append(bp, bp + decpt); +#endif +} + + +#if !defined (_STLP_NO_WCHAR_T) +void _STLP_CALL __convert_float_buffer( __iostring const& str, __iowstring &out, + const ctype& ct, wchar_t dot, bool __check_dot) { + string::const_iterator str_ite(str.begin()), str_end(str.end()); + + //First loop, check the dot char + if (__check_dot) { + while (str_ite != str_end) { + if (*str_ite != '.') { + out += ct.widen(*str_ite++); + } else { + out += dot; + break; + } + } + } else { + if (str_ite != str_end) { + out += ct.widen(*str_ite); + } + } + + if (str_ite != str_end) { + //Second loop, dot has been found, no check anymore + while (++str_ite != str_end) { + out += ct.widen(*str_ite); + } + } +} + +#endif + +void _STLP_CALL +__adjust_float_buffer(__iostring &str, char dot) { + if ('.' != dot) { + size_t __dot_pos = str.find('.'); + if (__dot_pos != string::npos) { + str[__dot_pos] = dot; + } + } +} + +_STLP_MOVE_TO_STD_NAMESPACE +_STLP_END_NAMESPACE + +// Local Variables: +// mode:C++ +// End: