1 /*
2 * Copyright (c) 1999
3 * Silicon Graphics Computer Systems, Inc.
4 *
5 * Copyright (c) 1999
6 * Boris Fomitchev
7 *
8 * This material is provided "as is", with absolutely no warranty expressed
9 * or implied. Any use is at your own risk.
10 *
11 * Permission to use or copy this software for any purpose is hereby granted
12 * without fee, provided the above notices are retained on all copies.
13 * Permission to modify the code and to distribute modified code is granted,
14 * provided the above notices are retained, and a notice that the code was
15 * modified is included with the above copyright notice.
16 *
17 */
18
19 #include "stlport_prefix.h"
20
21 #include <cmath>
22 #include <ios>
23 #include <locale>
24
25 #if defined (__DECCXX)
26 # define NDIG 400
27 #else
28 # define NDIG 82
29 #endif
30
31 #define todigit(x) ((x)+'0')
32
33 #if defined (_STLP_UNIX)
34
35 # if defined (__sun)
36 # include <floatingpoint.h>
37 # endif
38
39 # if defined (__sun) || defined (__digital__) || defined (__sgi) || defined (_STLP_SCO_OPENSERVER) || defined (__NCR_SVR)
40 // DEC, SGI & Solaris need this
41 # include <values.h>
42 # include <nan.h>
43 # endif
44
45 # if defined (__QNXNTO__) || ( defined(__GNUC__) && defined(__APPLE__) ) || defined(_STLP_USE_UCLIBC) /* 0.9.26 */ || \
46 defined(__FreeBSD__)
47 # define USE_SPRINTF_INSTEAD
48 # endif
49
50 # if defined (_AIX) // JFA 3-Aug-2000
51 # include <math.h>
52 # include <float.h>
53 # endif
54
55 # include <math.h>
56 #endif
57
58 #include <cstdio>
59 #include <cstdlib>
60
61 #if defined (_STLP_MSVC_LIB) || defined (__MINGW32__) || defined (__BORLANDC__) || defined (__DJGPP) || \
62 defined (_STLP_SCO_OPENSERVER) || defined (__NCR_SVR)
63 # include <float.h>
64 #endif
65
66 #if defined (__MRC__) || defined (__SC__) || defined (_CRAY) //*TY 02/24/2000 - added support for MPW
67 # include <fp.h>
68 #endif
69
70 #if defined (__CYGWIN__)
71 # include <ieeefp.h>
72 #endif
73
74 #if defined (__MSL__)
75 # include <cstdlib> // for atoi
76 # include <cstdio> // for snprintf
77 # include <algorithm>
78 # include <cassert>
79 #endif
80
81 #if defined (__ISCPP__)
82 # include <cfloat>
83 #endif
84
85 #include <algorithm>
86
87 #if defined (__DMC__)
88 # define snprintf _snprintf
89 #endif
90
91 _STLP_BEGIN_NAMESPACE
92
93 _STLP_MOVE_TO_PRIV_NAMESPACE
94
95 #if defined (__MWERKS__) || defined(__BEOS__)
96 # define USE_SPRINTF_INSTEAD
97 #endif
98
99 template <int N>
100 struct _Dig
101 {
102 enum { dig = _Dig<N/10>::dig + 1 };
103 };
104
105 _STLP_TEMPLATE_NULL
106 struct _Dig<0>
107 {
108 enum { dig = 0 };
109 };
110
111 #ifdef _STLP_NO_LONG_DOUBLE
112 # define MAXEDIGITS int(_Dig<DBL_MAX_10_EXP>::dig)
113 # define MAXFSIG DBL_DIG
114 # define MAXFCVT (DBL_DIG + 1)
115 #else
116 # define MAXEDIGITS int(_Dig<LDBL_MAX_10_EXP>::dig)
117 # define MAXFSIG LDBL_DIG
118 # define MAXFCVT (LDBL_DIG + 1)
119 #endif
120
121 // Tests for infinity and NaN differ on different OSs. We encapsulate
122 // these differences here.
123 #if !defined (USE_SPRINTF_INSTEAD)
124 # if defined (__hpux) && defined (__GNUC__)
125 # define _STLP_USE_SIGN_HELPER
126 # elif defined (__DJGPP) || (defined (_STLP_USE_GLIBC) && ! defined (__MSL__)) || \
127 defined (__CYGWIN__) || \
128 defined (__FreeBSD__) || defined (__NetBSD__) || defined (__OpenBSD__) || \
129 defined (__HP_aCC)
_Stl_is_nan_or_inf(double x)130 static inline bool _Stl_is_nan_or_inf(double x)
131 # if defined (isfinite)
132 { return !isfinite(x); }
133 # else
134 { return !finite(x); }
135 # endif
_Stl_is_neg_nan(double x)136 static inline bool _Stl_is_neg_nan(double x) { return isnan(x) && ( copysign(1., x) < 0 ); }
_Stl_is_inf(double x)137 static inline bool _Stl_is_inf(double x) { return isinf(x); }
138 // inline bool _Stl_is_neg_inf(double x) { return isinf(x) < 0; }
_Stl_is_neg_inf(double x)139 static inline bool _Stl_is_neg_inf(double x) { return isinf(x) && x < 0; }
140 # elif (defined (__unix) || defined (__unix__)) && \
141 !defined (__APPLE__) && !defined (__DJGPP) && !defined(__osf__) && \
142 !defined (_CRAY) && !defined (__ANDROID__)
_Stl_is_nan_or_inf(double x)143 static inline bool _Stl_is_nan_or_inf(double x) { return IsNANorINF(x); }
_Stl_is_inf(double x)144 static inline bool _Stl_is_inf(double x) { return IsNANorINF(x) && IsINF(x); }
_Stl_is_neg_inf(double x)145 static inline bool _Stl_is_neg_inf(double x) { return (IsINF(x)) && (x < 0.0); }
_Stl_is_neg_nan(double x)146 static inline bool _Stl_is_neg_nan(double x) { return IsNegNAN(x); }
147 # elif defined (_STLP_MSVC_LIB) || defined (__MINGW32__) || defined (__BORLANDC__)
_Stl_is_nan_or_inf(double x)148 static inline bool _Stl_is_nan_or_inf(double x) { return !_finite(x); }
149 # if !defined (__BORLANDC__)
_Stl_is_inf(double x)150 static inline bool _Stl_is_inf(double x) {
151 int fclass = _fpclass(x);
152 return fclass == _FPCLASS_NINF || fclass == _FPCLASS_PINF;
153 }
_Stl_is_neg_inf(double x)154 static inline bool _Stl_is_neg_inf(double x) { return _fpclass(x) == _FPCLASS_NINF; }
155 # else
_Stl_is_inf(double x)156 static inline bool _Stl_is_inf(double x) { return _Stl_is_nan_or_inf(x) && !_isnan(x);}
_Stl_is_neg_inf(double x)157 static inline bool _Stl_is_neg_inf(double x) { return _Stl_is_inf(x) && x < 0 ; }
158 # endif
_Stl_is_neg_nan(double x)159 static inline bool _Stl_is_neg_nan(double x) { return _isnan(x) && _copysign(1., x) < 0 ; }
160 # if defined (__BORLANDC__)
_Stl_is_nan_or_inf(long double x)161 static inline bool _Stl_is_nan_or_inf(long double x) { return !_finitel(x); }
_Stl_is_inf(long double x)162 static inline bool _Stl_is_inf(long double x) { return _Stl_is_nan_or_inf(x) && !_isnanl(x);}
_Stl_is_neg_inf(long double x)163 static inline bool _Stl_is_neg_inf(long double x) { return _Stl_is_inf(x) && x < 0 ; }
_Stl_is_neg_nan(long double x)164 static inline bool _Stl_is_neg_nan(long double x) { return _isnanl(x) && _copysignl(1.l, x) < 0 ; }
165 # elif !defined (_STLP_NO_LONG_DOUBLE)
166 // Simply there to avoid warning long double -> double implicit conversion:
_Stl_is_nan_or_inf(long double x)167 static inline bool _Stl_is_nan_or_inf(long double x) { return _Stl_is_nan_or_inf(__STATIC_CAST(double, x)); }
_Stl_is_inf(long double x)168 static inline bool _Stl_is_inf(long double x) { return _Stl_is_inf(__STATIC_CAST(double, x));}
_Stl_is_neg_inf(long double x)169 static inline bool _Stl_is_neg_inf(long double x) { return _Stl_is_neg_inf(__STATIC_CAST(double, x)); }
_Stl_is_neg_nan(long double x)170 static inline bool _Stl_is_neg_nan(long double x) { return _Stl_is_neg_nan(__STATIC_CAST(double, x)); }
171 # endif
172 # elif defined (__MRC__) || defined (__SC__) || defined (__DMC__)
_Stl_is_nan_or_inf(double x)173 static bool _Stl_is_nan_or_inf(double x) { return isnan(x) || !isfinite(x); }
_Stl_is_inf(double x)174 static bool _Stl_is_inf(double x) { return !isfinite(x); }
_Stl_is_neg_inf(double x)175 static bool _Stl_is_neg_inf(double x) { return !isfinite(x) && signbit(x); }
_Stl_is_neg_nan(double x)176 static bool _Stl_is_neg_nan(double x) { return isnan(x) && signbit(x); }
177 # elif /* defined(__FreeBSD__) || defined(__OpenBSD__) || */ (defined(__GNUC__) && defined(__APPLE__))
_Stl_is_nan_or_inf(double x)178 static inline bool _Stl_is_nan_or_inf(double x) { return !finite(x); }
_Stl_is_inf(double x)179 static inline bool _Stl_is_inf(double x) { return _Stl_is_nan_or_inf(x) && ! isnan(x); }
_Stl_is_neg_inf(double x)180 static inline bool _Stl_is_neg_inf(double x) { return _Stl_is_inf(x) && x < 0 ; }
_Stl_is_neg_nan(double x)181 static inline bool _Stl_is_neg_nan(double x) { return isnan(x) && copysign(1., x) < 0 ; }
182 # elif defined( _AIX ) // JFA 11-Aug-2000
_Stl_is_nan_or_inf(double x)183 static bool _Stl_is_nan_or_inf(double x) { return isnan(x) || !finite(x); }
_Stl_is_inf(double x)184 static bool _Stl_is_inf(double x) { return !finite(x); }
185 // bool _Stl_is_neg_inf(double x) { return _class(x) == FP_MINUS_INF; }
_Stl_is_neg_inf(double x)186 static bool _Stl_is_neg_inf(double x) { return _Stl_is_inf(x) && ( copysign(1., x) < 0 ); }
_Stl_is_neg_nan(double x)187 static bool _Stl_is_neg_nan(double x) { return isnan(x) && ( copysign(1., x) < 0 ); }
188 # elif defined (__ISCPP__)
_Stl_is_nan_or_inf(double x)189 static inline bool _Stl_is_nan_or_inf (double x) { return _fp_isINF(x) || _fp_isNAN(x); }
_Stl_is_inf(double x)190 static inline bool _Stl_is_inf (double x) { return _fp_isINF(x); }
_Stl_is_neg_inf(double x)191 static inline bool _Stl_is_neg_inf (double x) { return _fp_isINF(x) && x < 0; }
_Stl_is_neg_nan(double x)192 static inline bool _Stl_is_neg_nan (double x) { return _fp_isNAN(x) && x < 0; }
193 # elif defined (_CRAY)
194 # if defined (_CRAYIEEE)
_Stl_is_nan_or_inf(double x)195 static inline bool _Stl_is_nan_or_inf(double x) { return isnan(x) || isinf(x); }
_Stl_is_inf(double x)196 static inline bool _Stl_is_inf(double x) { return isinf(x); }
_Stl_is_neg_inf(double x)197 static inline bool _Stl_is_neg_inf(double x) { return isinf(x) && signbit(x); }
_Stl_is_neg_nan(double x)198 static inline bool _Stl_is_neg_nan(double x) { return isnan(x) && signbit(x); }
199 # else
_Stl_is_nan_or_inf(double x)200 static inline bool _Stl_is_nan_or_inf(double x) { return false; }
_Stl_is_inf(double x)201 static inline bool _Stl_is_inf(double x) { return false; }
_Stl_is_neg_inf(double x)202 static inline bool _Stl_is_neg_inf(double x) { return false; }
_Stl_is_neg_nan(double x)203 static inline bool _Stl_is_neg_nan(double x) { return false; }
204 # endif
205 # else // nothing from above
206 # define USE_SPRINTF_INSTEAD
207 # endif
208 #endif // !USE_SPRINTF_INSTEAD
209
210 #if !defined (USE_SPRINTF_INSTEAD)
211 // Reentrant versions of floating-point conversion functions. The argument
212 // lists look slightly different on different operating systems, so we're
213 // encapsulating the differences here.
214
215 # if defined (__CYGWIN__) || defined(__DJGPP)
_Stl_ecvtR(double x,int n,int * pt,int * sign,char * buf)216 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf)
217 { return ecvtbuf(x, n, pt, sign, buf); }
_Stl_fcvtR(double x,int n,int * pt,int * sign,char * buf)218 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf)
219 { return fcvtbuf(x, n, pt, sign, buf); }
220 # if !defined (_STLP_NO_LONG_DOUBLE)
221 # if defined (__CYGWIN__)
222 # define _STLP_EMULATE_LONG_DOUBLE_CVT
223 # else
_Stl_ecvtR(long double x,int n,int * pt,int * sign,char * buf)224 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf)
225 { return ecvtbuf(x, n, pt, sign, buf); }
_Stl_fcvtR(long double x,int n,int * pt,int * sign,char * buf)226 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf)
227 { return fcvtbuf(x, n, pt, sign, buf); }
228 # endif
229 # endif
230 # elif defined (_STLP_USE_GLIBC)
_Stl_ecvtR(double x,int n,int * pt,int * sign,char * buf,size_t bsize)231 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf, size_t bsize)
232 { return ecvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0; }
_Stl_fcvtR(double x,int n,int * pt,int * sign,char * buf,size_t bsize)233 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf, size_t bsize)
234 { return fcvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0; }
235 # ifndef _STLP_NO_LONG_DOUBLE
_Stl_ecvtR(long double x,int n,int * pt,int * sign,char * buf,size_t bsize)236 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf, size_t bsize)
237 { return qecvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0; }
_Stl_fcvtR(long double x,int n,int * pt,int * sign,char * buf,size_t bsize)238 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf, size_t bsize)
239 { return qfcvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0; }
240 # endif
241 # define _STLP_NEED_CVT_BUFFER_SIZE
242 # elif defined (__sun)
_Stl_ecvtR(double x,int n,int * pt,int * sign,char * buf)243 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf)
244 { return econvert(x, n, pt, sign, buf); }
_Stl_fcvtR(double x,int n,int * pt,int * sign,char * buf)245 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf)
246 { return fconvert(x, n, pt, sign, buf); }
247 # ifndef _STLP_NO_LONG_DOUBLE
_Stl_ecvtR(long double x,int n,int * pt,int * sign,char * buf)248 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf)
249 { return qeconvert(&x, n, pt, sign, buf); }
_Stl_fcvtR(long double x,int n,int * pt,int * sign,char * buf)250 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf)
251 { return qfconvert(&x, n, pt, sign, buf); }
252 # endif
253 # elif defined (__DECCXX)
_Stl_ecvtR(double x,int n,int * pt,int * sign,char * buf,size_t bsize)254 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf, size_t bsize)
255 { return (ecvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0); }
_Stl_fcvtR(double x,int n,int * pt,int * sign,char * buf,size_t bsize)256 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf, size_t bsize)
257 { return (fcvt_r(x, n, pt, sign, buf, bsize) == 0 ? buf : 0); }
258 # ifndef _STLP_NO_LONG_DOUBLE
259 // fbp : no "long double" conversions !
_Stl_ecvtR(long double x,int n,int * pt,int * sign,char * buf,size_t bsize)260 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf, size_t bsize)
261 { return (ecvt_r((double)x, n, pt, sign, buf, bsize) == 0 ? buf : 0) ; }
_Stl_fcvtR(long double x,int n,int * pt,int * sign,char * buf,size_t bsize)262 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf, size_t bsize)
263 { return (fcvt_r((double)x, n, pt, sign, buf, bsize) == 0 ? buf : 0); }
264 # endif
265 # define _STLP_NEED_CVT_BUFFER_SIZE
266 # elif defined (__hpux)
_Stl_ecvtR(double x,int n,int * pt,int * sign)267 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign)
268 { return ecvt(x, n, pt, sign); }
_Stl_fcvtR(double x,int n,int * pt,int * sign)269 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign)
270 { return fcvt(x, n, pt, sign); }
271 # if !defined (_STLP_NO_LONG_DOUBLE)
_Stl_ecvtR(long double x,int n,int * pt,int * sign)272 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign)
273 { return _ldecvt(*(long_double*)&x, n, pt, sign); }
_Stl_fcvtR(long double x,int n,int * pt,int * sign)274 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign)
275 { return _ldfcvt(*(long_double*)&x, n, pt, sign); }
276 # endif
277 # define _STLP_CVT_NEED_SYNCHRONIZATION
278 # elif defined (__unix) && !defined (__APPLE__) && !defined (_CRAY) && \
279 !defined (__ANDROID__)
_Stl_ecvtR(double x,int n,int * pt,int * sign,char * buf)280 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf)
281 { return ecvt_r(x, n, pt, sign, buf); }
_Stl_fcvtR(double x,int n,int * pt,int * sign,char * buf)282 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf)
283 { return fcvt_r(x, n, pt, sign, buf); }
284 # if !defined (_STLP_NO_LONG_DOUBLE)
_Stl_ecvtR(long double x,int n,int * pt,int * sign,char * buf)285 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf)
286 { return qecvt_r(x, n, pt, sign, buf); }
_Stl_fcvtR(long double x,int n,int * pt,int * sign,char * buf)287 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf)
288 { return qfcvt_r(x, n, pt, sign, buf); }
289 # endif
290 # elif defined (_STLP_MSVC_LIB) || defined (__MINGW32__) || defined (__BORLANDC__)
291 # if defined (_STLP_USE_SAFE_STRING_FUNCTIONS)
292 # define _STLP_APPEND(a, b) a##b
293 # define _STLP_BUF_PARAMS , char* buf, size_t bsize
294 # define _STLP_SECURE_FUN(F, X, N, PT, SIGN) _STLP_APPEND(F, _s)(buf, bsize, X, N, PT, SIGN); return buf
295 # else
296 # define _STLP_BUF_PARAMS
297 # define _STLP_SECURE_FUN(F, X, N, PT, SIGN) return F(X, N, PT, SIGN)
298 # define _STLP_CVT_NEED_SYNCHRONIZATION
299 # endif
_Stl_ecvtR(double x,int n,int * pt,int * sign _STLP_BUF_PARAMS)300 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign _STLP_BUF_PARAMS)
301 { _STLP_SECURE_FUN(_ecvt, x, n, pt, sign); }
_Stl_fcvtR(double x,int n,int * pt,int * sign _STLP_BUF_PARAMS)302 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign _STLP_BUF_PARAMS)
303 { _STLP_SECURE_FUN(_fcvt, x, n, pt, sign); }
304 # if !defined (_STLP_NO_LONG_DOUBLE)
305 # if defined (_STLP_USE_SAFE_STRING_FUNCTIONS)
306 # define _STLP_PARAMS , buf, bsize
307 # else
308 # define _STLP_PARAMS
309 # endif
_Stl_ecvtR(long double x,int n,int * pt,int * sign _STLP_BUF_PARAMS)310 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign _STLP_BUF_PARAMS)
311 { return _Stl_ecvtR(__STATIC_CAST(double, x), n, pt, sign _STLP_PARAMS); }
_Stl_fcvtR(long double x,int n,int * pt,int * sign _STLP_BUF_PARAMS)312 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign _STLP_BUF_PARAMS)
313 { return _Stl_fcvtR(__STATIC_CAST(double, x), n, pt, sign _STLP_PARAMS); }
314 # undef _STLP_PARAMS
315 # endif
316 # undef _STLP_SECURE_FUN
317 # undef _STLP_BUF_PARAMS
318 # undef _STLP_APPEND
319 # if defined (__BORLANDC__) /* || defined (__GNUC__) MinGW do not support 'L' modifier so emulation do not work */
320 # define _STLP_EMULATE_LONG_DOUBLE_CVT
321 # endif
322 # elif defined (__ISCPP__)
_Stl_ecvtR(double x,int n,int * pt,int * sign,char * buf)323 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf)
324 { return _fp_ecvt( x, n, pt, sign, buf); }
_Stl_fcvtR(double x,int n,int * pt,int * sign,char * buf)325 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf)
326 { return _fp_fcvt(x, n, pt, sign, buf); }
327 # if !defined (_STLP_NO_LONG_DOUBLE)
_Stl_ecvtR(long double x,int n,int * pt,int * sign,char * buf)328 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf)
329 { return _fp_ecvt( x, n, pt, sign, buf); }
_Stl_fcvtR(long double x,int n,int * pt,int * sign,char * buf)330 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf)
331 { return _fp_fcvt(x, n, pt, sign, buf); }
332 # endif
333 # elif defined (_AIX) || defined (__FreeBSD__) || defined (__NetBSD__) || defined (__OpenBSD__) || \
334 defined (__MRC__) || defined (__SC__) || defined (_CRAY) || \
335 defined (_STLP_SCO_OPENSERVER) || defined (__NCR_SVR) || \
336 defined (__DMC__)
_Stl_ecvtR(double x,int n,int * pt,int * sign)337 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign)
338 { return ecvt(x, n, pt, sign ); }
_Stl_fcvtR(double x,int n,int * pt,int * sign)339 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign)
340 { return fcvt(x, n, pt, sign); }
341 # if !defined (_STLP_NO_LONG_DOUBLE)
_Stl_ecvtR(long double x,int n,int * pt,int * sign)342 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign)
343 { return ecvt(x, n, pt, sign ); }
_Stl_fcvtR(long double x,int n,int * pt,int * sign)344 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign)
345 { return fcvt(x, n, pt, sign); }
346 # endif
347 # define _STLP_CVT_NEED_SYNCHRONIZATION
348 # else
349 # error Missing _Stl_ecvtR and _Stl_fcvtR implementations.
350 # endif
351
352 #if defined (_STLP_CVT_NEED_SYNCHRONIZATION)
353 /* STLport synchronize access to *cvt functions but those methods might
354 * be called from outside, in this case we will still have a race condition. */
355 # if defined (_STLP_THREADS)
put_float_mutex()356 static _STLP_STATIC_MUTEX& put_float_mutex() {
357 static _STLP_STATIC_MUTEX __put_float_mutex _STLP_MUTEX_INITIALIZER;
358 return __put_float_mutex;
359 }
_Stl_ecvtR(double x,int n,int * pt,int * sign,char * buf)360 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf) {
361 _STLP_auto_lock lock(put_float_mutex());
362 strcpy(buf, _Stl_ecvtR(x, n, pt, sign)); return buf;
363 }
_Stl_fcvtR(double x,int n,int * pt,int * sign,char * buf)364 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf) {
365 _STLP_auto_lock lock(put_float_mutex());
366 strcpy(buf, _Stl_fcvtR(x, n, pt, sign)); return buf;
367 }
368 # if !defined (_STLP_NO_LONG_DOUBLE) && !defined (_STLP_EMULATE_LONG_DOUBLE_CVT)
_Stl_ecvtR(long double x,int n,int * pt,int * sign,char * buf)369 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf) {
370 _STLP_auto_lock lock(put_float_mutex());
371 strcpy(buf, _Stl_ecvtR(x, n, pt, sign)); return buf;
372 }
_Stl_fcvtR(long double x,int n,int * pt,int * sign,char * buf)373 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf) {
374 _STLP_auto_lock lock(put_float_mutex());
375 strcpy(buf, _Stl_fcvtR(x, n, pt, sign)); return buf;
376 }
377 # endif
378 # else
_Stl_ecvtR(double x,int n,int * pt,int * sign,char *)379 static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char*)
380 { return _Stl_ecvtR(x, n, pt, sign); }
_Stl_fcvtR(double x,int n,int * pt,int * sign,char *)381 static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char*)
382 { return _Stl_fcvtR(x, n, pt, sign); }
383 # if !defined (_STLP_NO_LONG_DOUBLE) && !defined (_STLP_EMULATE_LONG_DOUBLE_CVT)
_Stl_ecvtR(long double x,int n,int * pt,int * sign,char *)384 static inline char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char*)
385 { return _Stl_ecvtR(x, n, pt, sign); }
_Stl_fcvtR(long double x,int n,int * pt,int * sign,char *)386 static inline char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char*)
387 { return _Stl_fcvtR(x, n, pt, sign); }
388 # endif
389 # endif
390 #endif
391
392 # if !defined (_STLP_USE_SAFE_STRING_FUNCTIONS) && !defined (_STLP_NEED_CVT_BUFFER_SIZE)
393 # define _STLP_CVT_BUFFER(B) B
394 # else
395 # define _STLP_CVT_BUFFER(B) _STLP_ARRAY_AND_SIZE(B)
396 # endif
397
398 # if defined (_STLP_EMULATE_LONG_DOUBLE_CVT)
399 static void __fill_fmtbuf(char* fmtbuf, ios_base::fmtflags flags, char long_modifier);
400
401 // Emulation of ecvt/fcvt functions using sprintf:
_Stl_ecvtR(long double x,int n,int * pt,int * sign,char * buf)402 static char* _Stl_ecvtR(long double x, int n, int* pt, int* sign, char* buf) {
403 // If long double value can be safely converted to double without losing precision
404 // we use the ecvt function for double:
405 double y = __STATIC_CAST(double, x);
406 if (x == y)
407 return _Stl_ecvtR(y, n, pt, sign, buf);
408
409 char fmtbuf[32];
410 __fill_fmtbuf(fmtbuf, 0, 'L');
411 sprintf(buf, fmtbuf, n, x < 0.0l ? -x : x);
412 /* We are waiting for something having the form x.xxxe+yyyy */
413 *pt = 0;
414 *sign = 0;
415 int i = -1;
416 int offset = 0;
417 while (buf[++i] != 0 && n != 0) {
418 if (buf[i] >= '0' && buf[i] <= '9') {
419 --n;
420 if (offset != 0)
421 buf[i - offset] = buf[i];
422 }
423 else {
424 if (offset != 0) break;
425 ++offset;
426 *pt = i;
427 }
428 }
429 if (offset != 0)
430 buf[i - offset] = 0;
431 // Extract exponent part in point position:
432 int e = 0;
433 while (buf[++i] != 0) {
434 if (buf[i] >= '0' && buf[i] <= '9') {
435 e = e * 10 + (buf[i] - '0');
436 }
437 }
438 *pt += e;
439 return buf;
440 }
441
_Stl_fcvtR(long double x,int n,int * pt,int * sign,char * buf)442 static char* _Stl_fcvtR(long double x, int n, int* pt, int* sign, char* buf) {
443 // If long double value can be safely converted to double without losing precision
444 // we use the fcvt function for double:
445 double y = __STATIC_CAST(double, x);
446 if (x == y)
447 return _Stl_fcvtR(y, n, pt, sign, buf);
448
449 char fmtbuf[32];
450 __fill_fmtbuf(fmtbuf, ios_base::fixed, 'L');
451 sprintf(buf, fmtbuf, n, x < 0.0l ? -x : x);
452 *pt = 0;
453 *sign = 0;
454 int i = -1;
455 int offset = 0;
456 while (buf[++i] != 0 && (offset == 0 || n != 0)) {
457 if (buf[i] >= '0' && buf[i] <= '9') {
458 if (offset != 0) {
459 --n;
460 buf[i - offset] = buf[i];
461 }
462 }
463 else {
464 ++offset;
465 *pt = i;
466 }
467 }
468 if (offset != 0)
469 buf[i - offset] = 0;
470 else
471 *pt = i;
472 return buf;
473 }
474 #endif
475
476 //----------------------------------------------------------------------
477 // num_put
478
479 // __format_float formats a mantissa and exponent as returned by
480 // one of the conversion functions (ecvt_r, fcvt_r, qecvt_r, qfcvt_r)
481 // according to the specified precision and format flags. This is
482 // based on doprnt but is much simpler since it is concerned only
483 // with floating point input and does not consider all formats. It
484 // also does not deal with blank padding, which is handled by
485 // __copy_float_and_fill.
486
__format_float_scientific(__iostring & buf,const char * bp,int decpt,int sign,bool is_zero,ios_base::fmtflags flags,int precision)487 static size_t __format_float_scientific( __iostring& buf, const char *bp,
488 int decpt, int sign, bool is_zero,
489 ios_base::fmtflags flags,
490 int precision) {
491 // sign if required
492 if (sign)
493 buf += '-';
494 else if (flags & ios_base::showpos)
495 buf += '+';
496
497 // first digit of mantissa
498 buf += *bp++;
499
500 // start of grouping position, grouping won't occur in scientific notation
501 // as it is impossible to have something like 1234.0e04 but we return a correct
502 // group position for coherency with __format_float_fixed.
503 size_t __group_pos = buf.size();
504
505 // decimal point if required
506 if (precision != 0 || flags & ios_base::showpoint) {
507 buf += '.';
508 }
509
510 // rest of mantissa
511 while (*bp != 0 && precision--)
512 buf += *bp++;
513
514 // trailing 0 if needed
515 if (precision > 0)
516 buf.append(precision, '0');
517
518 // exponent size = number of digits + exponent sign + exponent symbol + trailing zero
519 char expbuf[MAXEDIGITS + 3];
520 //We start filling at the buffer end
521 char *suffix = expbuf + MAXEDIGITS + 2;
522 *suffix = 0;
523 if (!is_zero) {
524 int nn = decpt - 1;
525 if (nn < 0)
526 nn = -nn;
527 for (; nn > 9; nn /= 10)
528 *--suffix = (char) todigit(nn % 10);
529 *--suffix = (char) todigit(nn);
530 }
531
532 // prepend leading zeros to exponent
533 // C89 Standard says that it should be at least 2 digits, C99 Standard says that
534 // we stop prepend zeros if more than 3 digits. To repect both STLport prepend zeros
535 // until it is 2 digits.
536 while (suffix > &expbuf[MAXEDIGITS])
537 *--suffix = '0';
538
539 // put in the exponent sign
540 *--suffix = (char) ((decpt > 0 || is_zero ) ? '+' : '-');
541
542 // put in the e
543 *--suffix = flags & ios_base::uppercase ? 'E' : 'e';
544
545 // copy the suffix
546 buf += suffix;
547 return __group_pos;
548 }
549
__format_float_fixed(__iostring & buf,const char * bp,int decpt,int sign,ios_base::fmtflags flags,int precision)550 static size_t __format_float_fixed( __iostring &buf, const char *bp,
551 int decpt, int sign,
552 ios_base::fmtflags flags,
553 int precision) {
554 if ( sign && (decpt > -precision) && (*bp != 0) )
555 buf += '-';
556 else if ( flags & ios_base::showpos )
557 buf += '+';
558
559 // digits before decimal point
560 int nnn = decpt;
561 do {
562 buf += (nnn <= 0 || *bp == 0) ? '0' : *bp++;
563 } while ( --nnn > 0 );
564
565 // start of grouping position
566 size_t __group_pos = buf.size();
567
568 // decimal point if needed
569 if ( flags & ios_base::showpoint || precision > 0 ) {
570 buf += '.';
571 }
572
573 // digits after decimal point if any
574 while ( *bp != 0 && --precision >= 0 ) {
575 buf += (++decpt <= 0) ? '0' : *bp++;
576 }
577
578 // trailing zeros if needed
579 if (precision > 0)
580 buf.append(precision, '0');
581
582 return __group_pos;
583 }
584
585 #if defined (_STLP_USE_SIGN_HELPER)
586 template<class _FloatT>
587 struct float_sign_helper {
float_sign_helperfloat_sign_helper588 float_sign_helper(_FloatT __x)
589 { _M_number._num = __x; }
590
is_negativefloat_sign_helper591 bool is_negative() const {
592 const unsigned short sign_mask(1 << (sizeof(unsigned short) * CHAR_BIT - 1));
593 return (get_sign_word() & sign_mask) != 0;
594 }
595 private:
596 union {
597 unsigned short _Words[8];
598 _FloatT _num;
599 } _M_number;
600
get_word_higherfloat_sign_helper601 unsigned short get_word_higher() const _STLP_NOTHROW
602 { return _M_number._Words[0]; }
get_word_lowerfloat_sign_helper603 unsigned short get_word_lower() const _STLP_NOTHROW
604 { return _M_number._Words[(sizeof(_FloatT) >= 12 ? 10 : sizeof(_FloatT)) / sizeof(unsigned short) - 1]; }
get_sign_wordfloat_sign_helper605 unsigned short get_sign_word() const _STLP_NOTHROW
606 # if defined (_STLP_BIG_ENDIAN)
607 { return get_word_higher(); }
608 # else /* _STLP_LITTLE_ENDIAN */
609 { return get_word_lower(); }
610 # endif
611 };
612 #endif
613
614 template <class _FloatT>
__format_nan_or_inf(__iostring & buf,_FloatT x,ios_base::fmtflags flags)615 static size_t __format_nan_or_inf(__iostring& buf, _FloatT x, ios_base::fmtflags flags) {
616 static const char* inf[2] = { "inf", "Inf" };
617 static const char* nan[2] = { "nan", "NaN" };
618 const char** inf_or_nan;
619 #if !defined (_STLP_USE_SIGN_HELPER)
620 if (_Stl_is_inf(x)) { // Infinity
621 inf_or_nan = inf;
622 if (_Stl_is_neg_inf(x))
623 buf += '-';
624 else if (flags & ios_base::showpos)
625 buf += '+';
626 } else { // NaN
627 inf_or_nan = nan;
628 if (_Stl_is_neg_nan(x))
629 buf += '-';
630 else if (flags & ios_base::showpos)
631 buf += '+';
632 }
633 #else
634 typedef numeric_limits<_FloatT> limits;
635 if (x == limits::infinity() || x == -limits::infinity()) {
636 inf_or_nan = inf;
637 } else { // NaN
638 inf_or_nan = nan;
639 }
640 float_sign_helper<_FloatT> helper(x);
641 if (helper.is_negative())
642 buf += '-';
643 else if (flags & ios_base::showpos)
644 buf += '+';
645 #endif
646 size_t ret = buf.size();
647 buf += inf_or_nan[flags & ios_base::uppercase ? 1 : 0];
648 return ret;
649 }
650
__format_float(__iostring & buf,const char * bp,int decpt,int sign,bool is_zero,ios_base::fmtflags flags,int precision)651 static inline size_t __format_float(__iostring &buf, const char * bp,
652 int decpt, int sign, bool is_zero,
653 ios_base::fmtflags flags,
654 int precision) {
655 size_t __group_pos = 0;
656 switch (flags & ios_base::floatfield) {
657 case ios_base::scientific:
658 __group_pos = __format_float_scientific( buf, bp, decpt, sign, is_zero,
659 flags, precision);
660 break;
661 case ios_base::fixed:
662 __group_pos = __format_float_fixed( buf, bp, decpt, sign,
663 flags, precision);
664 break;
665 default: // g format
666 // establish default precision
667 if (flags & ios_base::showpoint || precision > 0) {
668 if (precision == 0) precision = 1;
669 } else
670 precision = 6;
671
672 // reset exponent if value is zero
673 if (is_zero)
674 decpt = 1;
675
676 int kk = precision;
677 if (!(flags & ios_base::showpoint)) {
678 size_t n = strlen(bp);
679 if (n < (size_t)kk)
680 kk = (int)n;
681 while (kk >= 1 && bp[kk-1] == '0')
682 --kk;
683 }
684
685 if (decpt < -3 || decpt > precision) {
686 precision = kk - 1;
687 __group_pos = __format_float_scientific( buf, bp, decpt, sign, is_zero,
688 flags, precision);
689 } else {
690 precision = kk - decpt;
691 __group_pos = __format_float_fixed( buf, bp, decpt, sign,
692 flags, precision);
693 }
694 break;
695 } /* switch */
696 return __group_pos;
697 }
698
699 #endif
700
701 #if defined (USE_SPRINTF_INSTEAD) || defined (_STLP_EMULATE_LONG_DOUBLE_CVT)
702 struct GroupPos {
operator ()GroupPos703 bool operator () (char __c) const {
704 return __c == '.' ||
705 __c == 'e' || __c == 'E';
706 }
707 };
708
709 // Creates a format string for sprintf()
__fill_fmtbuf(char * fmtbuf,ios_base::fmtflags flags,char long_modifier)710 static void __fill_fmtbuf(char* fmtbuf, ios_base::fmtflags flags, char long_modifier) {
711 fmtbuf[0] = '%';
712 int i = 1;
713
714 if (flags & ios_base::showpos)
715 fmtbuf[i++] = '+';
716
717 if (flags & ios_base::showpoint)
718 fmtbuf[i++] = '#';
719
720 fmtbuf[i++] = '.';
721 fmtbuf[i++] = '*';
722
723 if (long_modifier)
724 fmtbuf[i++] = long_modifier;
725
726 switch (flags & ios_base::floatfield)
727 {
728 case ios_base::scientific:
729 fmtbuf[i++] = (flags & ios_base::uppercase) ? 'E' : 'e';
730 break;
731 case ios_base::fixed:
732 # if defined (__FreeBSD__)
733 fmtbuf[i++] = 'f';
734 # else
735 fmtbuf[i++] = (flags & ios_base::uppercase) ? 'F' : 'f';
736 # endif
737 break;
738 default:
739 fmtbuf[i++] = (flags & ios_base::uppercase) ? 'G' : 'g';
740 break;
741 }
742
743 fmtbuf[i] = 0;
744 }
745
746 #endif /* USE_SPRINTF_INSTEAD */
747
748 template <class _FloatT>
__write_floatT(__iostring & buf,ios_base::fmtflags flags,int precision,_FloatT x,char modifier)749 static size_t __write_floatT(__iostring &buf, ios_base::fmtflags flags, int precision,
750 _FloatT x
751 #if defined (USE_SPRINTF_INSTEAD)
752 , char modifier) {
753 /* In theory, if we want 'arbitrary' precision, we should use 'arbitrary'
754 * buffer size below, but really we limited by exponent part in double.
755 * - ptr
756 */
757 typedef numeric_limits<_FloatT> limits;
758 char static_buf[limits::max_exponent10 + 6]; // 6: -xxx.yyyE-zzz (sign, dot, E, exp sign, \0)
759 char fmtbuf[32];
760 __fill_fmtbuf(fmtbuf, flags, modifier);
761 snprintf(_STLP_ARRAY_AND_SIZE(static_buf), fmtbuf, precision, x);
762 buf = static_buf;
763 return find_if(buf.begin(), buf.end(), GroupPos()) - buf.begin();
764 #else
765 ) {
766 typedef numeric_limits<_FloatT> limits;
767 //If numeric_limits support is correct we use the exposed values to detect NaN and infinity:
768 if (limits::has_infinity && limits::has_quiet_NaN) {
769 if (!(x == x) || // NaN check
770 (x == limits::infinity() || x == -limits::infinity())) {
771 return __format_nan_or_inf(buf, x, flags);
772 }
773 }
774 // numeric_limits support is not good enough, we rely on platform dependent function
775 // _Stl_is_nan_or_inf that do not support long double.
776 else if (_Stl_is_nan_or_inf(x)) {
777 return __format_nan_or_inf(buf, x, flags);
778 }
779 # if defined (__MINGW32__)
780 //For the moment MinGW is limited to display at most numeric_limits<double>::max()
781 if (x > numeric_limits<double>::max() ||
782 x < -numeric_limits<double>::max()) {
783 return __format_nan_or_inf(buf, x, flags);
784 }
785 # endif
786
787 /* Buffer size is max number of digits which is the addition of:
788 * - max_exponent10: max number of digits in fixed mode
789 * - digits10 + 2: max number of significant digits
790 * - trailing '\0'
791 */
792 char cvtbuf[limits::max_exponent10 + limits::digits10 + 2 + 1];
793 char *bp;
794 int decpt, sign;
795
796 switch (flags & ios_base::floatfield) {
797 case ios_base::fixed:
798 {
799 /* Here, number of digits represents digits _after_ decimal point.
800 * In order to limit static buffer size we have to give 2 different values depending on x value.
801 * For small values (abs(x) < 1) we need as many digits as requested by precision limited by the maximum number of digits
802 * which is min_exponent10 + digits10 + 2
803 * For bigger values we won't have more than limits::digits10 + 2 digits after decimal point. */
804 int digits10 = (x > -1.0 && x < 1.0 ? -limits::min_exponent10 + limits::digits10 + 2
805 : limits::digits10 + 2);
806 bp = _Stl_fcvtR(x, (min) (precision, digits10), &decpt, &sign, _STLP_CVT_BUFFER(cvtbuf) );
807 }
808 break;
809 case ios_base::scientific:
810 default:
811 /* Here, number of digits is total number of digits which is limited to digits10 + 2. */
812 {
813 int digits10 = limits::digits10 + 2;
814 bp = _Stl_ecvtR(x, (min) (precision, digits10), &decpt, &sign, _STLP_CVT_BUFFER(cvtbuf) );
815 }
816 break;
817 }
818 return __format_float(buf, bp, decpt, sign, x == 0.0, flags, precision);
819 #endif
820 }
821
822 size_t _STLP_CALL
823 __write_float(__iostring &buf, ios_base::fmtflags flags, int precision,
824 double x) {
825 return __write_floatT(buf, flags, precision, x
826 #if defined (USE_SPRINTF_INSTEAD)
827 , 0
828 #endif
829 );
830 }
831
832 #if !defined (_STLP_NO_LONG_DOUBLE)
833 size_t _STLP_CALL
834 __write_float(__iostring &buf, ios_base::fmtflags flags, int precision,
835 long double x) {
836 return __write_floatT(buf, flags, precision, x
837 #if defined (USE_SPRINTF_INSTEAD)
838 , 'L'
839 #endif
840 );
841 }
842 #endif
843
844 void _STLP_CALL __get_floor_digits(__iostring &out, _STLP_LONGEST_FLOAT_TYPE __x) {
845 typedef numeric_limits<_STLP_LONGEST_FLOAT_TYPE> limits;
846 #if defined (USE_SPRINTF_INSTEAD)
847 char cvtbuf[limits::max_exponent10 + 6];
848 # if !defined (_STLP_NO_LONG_DOUBLE)
849 snprintf(_STLP_ARRAY_AND_SIZE(cvtbuf), "%Lf", __x); // check for 1234.56!
850 # else
851 snprintf(_STLP_ARRAY_AND_SIZE(cvtbuf), "%f", __x); // check for 1234.56!
852 # endif
853 char *p = strchr( cvtbuf, '.' );
854 if ( p == 0 ) {
855 out.append( cvtbuf );
856 } else {
857 out.append( cvtbuf, p );
858 }
859 #else
860 char cvtbuf[limits::max_exponent10 + 1];
861 char * bp;
862 int decpt, sign;
863 bp = _Stl_fcvtR(__x, 0, &decpt, &sign, _STLP_CVT_BUFFER(cvtbuf));
864
865 if (sign) {
866 out += '-';
867 }
868 out.append(bp, bp + decpt);
869 #endif
870 }
871
872
873 #if !defined (_STLP_NO_WCHAR_T)
874 void _STLP_CALL __convert_float_buffer( __iostring const& str, __iowstring &out,
875 const ctype<wchar_t>& ct, wchar_t dot, bool __check_dot) {
876 string::const_iterator str_ite(str.begin()), str_end(str.end());
877
878 //First loop, check the dot char
879 if (__check_dot) {
880 while (str_ite != str_end) {
881 if (*str_ite != '.') {
882 out += ct.widen(*str_ite++);
883 } else {
884 out += dot;
885 break;
886 }
887 }
888 } else {
889 if (str_ite != str_end) {
890 out += ct.widen(*str_ite);
891 }
892 }
893
894 if (str_ite != str_end) {
895 //Second loop, dot has been found, no check anymore
896 while (++str_ite != str_end) {
897 out += ct.widen(*str_ite);
898 }
899 }
900 }
901
902 #endif
903
904 void _STLP_CALL
905 __adjust_float_buffer(__iostring &str, char dot) {
906 if ('.' != dot) {
907 size_t __dot_pos = str.find('.');
908 if (__dot_pos != string::npos) {
909 str[__dot_pos] = dot;
910 }
911 }
912 }
913
914 _STLP_MOVE_TO_STD_NAMESPACE
915 _STLP_END_NAMESPACE
916
917 // Local Variables:
918 // mode:C++
919 // End:
920