/*

 * 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 <cmath>

#include <ios>

#include <locale>

#if defined (__DECCXX)

#  define NDIG 400

#else

#  define NDIG 82

#endif

#if defined (_STLP_NO_LONG_DOUBLE)

#  define MAXECVT 17

#  define MAXFCVT 18

typedef double max_double_type;

#else

#  define MAXECVT 35

#  define MAXFCVT 36

typedef long double max_double_type;

#endif

#define MAXFSIG MAXECVT

#define MAXESIZ 5

#define todigit(x) ((x)+'0')

#if defined (_STLP_UNIX)

#  if defined (__sun)

#    include <floatingpoint.h>

#  endif

#  if defined (__sun) || defined (__digital__) || defined (__sgi) || defined (_STLP_SCO_OPENSERVER) || defined (__NCR_SVR)

// DEC, SGI & Solaris need this

#    include <values.h>

#    include <nan.h>

#  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 <math.h>

#    include <float.h>

#  endif

#endif

#include <cstdio>

#include <cstdlib>

//#if defined(_CRAY)

//# include <stdlib.h>

//#endif

#if defined (_STLP_MSVC_LIB) || defined (__MINGW32__) || defined (__BORLANDC__) || defined (__DJGPP) || \

    defined (_STLP_SCO_OPENSERVER) || defined (__NCR_SVR)

#  include <float.h>

#endif

#if defined(__MRC__) || defined(__SC__)  || defined(_CRAY)  //*TY 02/24/2000 - added support for MPW

#  include <fp.h>

#endif

#if defined (__CYGWIN__)

#  include <ieeefp.h>

#endif

#if defined (__MSL__)

#  include <cstdlib>  // for atoi

#  include <cstdio>  // for snprintf

#  include <algorithm>

#  include <cassert>

#endif

#if defined (__ISCPP__)

#  include <cfloat>

#endif

#include <algorithm>

#if defined (__DMC__)

#  define snprintf _snprintf

#endif

#if defined(__hpux) && (!defined(_INCLUDE_HPUX_SOURCE) || defined(__GNUC__))

extern "C" double erf(double);

extern "C" double erfc(double);

extern "C" double gamma(double);                             /* obsolescent */

extern "C" double hypot(double, double);

extern "C" int    isnan(double);

extern "C" double j0(double);

extern "C" double j1(double);

extern "C" double jn(int, double);

extern "C" double lgamma(double);

extern "C" double y0(double);

extern "C" double y1(double);

extern "C" double yn(int, double);

#  define HUGE_VALF     _SINFINITY

#  define INFINITY      _SINFINITY

#  define NAN           _SQNAN

#  define isnan(x)              _ISNAN(x)

#  define isinf(x)              _ISINF(x)

#  define signbit(x)            _SIGNBIT(x)

#  define isfinite(x)           _ISFINITE(x)

#  define isnormal(x)           _ISNORMAL(x)

#  define fpclassify(x)         _FPCLASSIFY(x)

#  define isunordered(x,y)      _ISUNORDERED(x,y)

#  define isgreater(x,y)        _ISGREATER(x,y)

#  define isgreaterequal(x,y)   _ISGREATEREQUAL(x,y)

#  define isless(x,y)           _ISLESS(x,y)

#  define islessequal(x,y)      _ISLESSEQUAL(x,y)

#  define islessgreater(x,y)    _ISLESSGREATER(x,y)

#  define FP_NORMAL     0

#  define FP_ZERO       1

#  define FP_INFINITE   2

#  define FP_SUBNORMAL  3

#  define FP_NAN        4

#  define DECIMAL_DIG   17

#  define _IS64(x) (sizeof(x) == sizeof(double))

#  define _IS32(x) (sizeof(x) == sizeof(float))

extern "C" {

  extern double copysign(double, double);

  extern const float _SINFINITY;

  extern const float _SQNAN;

#  if defined (_PA_RISC)

#    define _ISNAN(x)          (_IS32(x)?_Isnanf(x):(isnan)(x))

#    define _ISINF(x)          (_IS32(x)?_Isinff(x):_Isinf(x))

#    define _SIGNBIT(x)        (_IS32(x)?_Signbitf(x):_Signbit(x))

#    define _ISFINITE(x)       (_IS32(x)?_Isfinitef(x):_Isfinite(x))

#    define _ISNORMAL(x)       (_IS32(x)?_Isnormalf(x):_Isnormal(x))

#    define _FPCLASSIFY(x)     (_IS32(x)?_Fpclassifyf(x)>>1:_Fpclassify(x)>>1)

#    define _ISUNORDERED(x,y)  (_IS32(x)&&_IS32(y)?_Isunorderedf(x,y):_Isunordered(x,y))

  extern int _Signbit(double);

  extern int _Signbitf(float);

  extern int _Isnanf(float);

  extern int _Isfinite(double);

  extern int _Isfinitef(float);

  extern int _Isinf(double);

  extern int _Isinff(float);

  extern int _Isnormal(double);

  extern int _Isnormalf(float);

  extern int _Isunordered(double, double);

  extern int _Isunorderedf(float, float);

  extern int _Fpclassify(double);

  extern int _Fpclassifyf(float);

#  else

#    include "math_ia64_internal.h"

#    define _FPCLASSIFY(x)     (_IS32(x)?_Fpclassf(x):_Fpclass(x))

  extern int _Fpclass(double);

  extern int _Fpclassf(float);

#  endif

}

#  if !defined (_INCLUDE_XOPEN_SOURCE_EXTENDED)

extern "C" char *fcvt(double, int, int *, int *);

extern "C" char *ecvt(double, int, int *, int *);

#  endif

#  if !defined (_INCLUDE_HPUX_SOURCE)

#    if !defined (_LONG_DOUBLE)

#      define _LONG_DOUBLE

typedef struct {

  uint32_t word1, word2, word3, word4;

} long_double;

#    endif /* _LONG_DOUBLE */

extern "C" char *_ldecvt(long_double, int, int *, int *);

extern "C" char *_ldfcvt(long_double, int, int *, int *);

#  endif

#endif /* __hpux */

_STLP_BEGIN_NAMESPACE

_STLP_MOVE_TO_PRIV_NAMESPACE

#if defined (__MWERKS__) || defined(__BEOS__)

#  define USE_SPRINTF_INSTEAD

#endif

#if defined (_AIX) || defined(__FreeBSD__) || defined(__NetBSD__) || defined(__OpenBSD__)

// Some OS'es only provide non-reentrant primitives, so we have to use additional synchronization here

#  if !defined(_REENTRANT) && !defined(_THREAD_SAFE) && !(defined(_POSIX_THREADS) && defined(__OpenBSD__))

#    define LOCK_CVT

#    define RETURN_CVT(ecvt, x, n, pt, sign, buf) return ecvt(x, n, pt, sign);

#  else

static _STLP_STATIC_MUTEX __put_float_mutex _STLP_MUTEX_INITIALIZER;

#    define LOCK_CVT _STLP_auto_lock lock(__put_float_mutex);

#    define RETURN_CVT(ecvt, x, n, pt, sign, buf) strcpy(buf, ecvt(x, n, pt, sign)); return buf;

#  endif // !_REENTRANT

#endif // _AIX || __FreeBSD__ || __NetBSD__ || __OpenBSD__

// Tests for infinity and NaN differ on different OSs.  We encapsulate

// these differences here.

#if !defined (USE_SPRINTF_INSTEAD)

#  if defined (__hpux) || defined (__DJGPP) || (defined (_STLP_USE_GLIBC) && ! defined (__MSL__)) || \

      defined (__FreeBSD__) || defined (__NetBSD__) || defined (__OpenBSD__)

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)

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 (__BORLANDC__) && ( __BORLANDC__ < 0x540 )

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) && x < 0 ; }

#  elif defined (_STLP_MSVC_LIB) || (defined (__MINGW32__) && !(defined(__SYMBIAN32__) && defined(__GCCXML__)) ) || defined (__BORLANDC__)

static inline bool _Stl_is_nan_or_inf(double x) { return !_finite(x); }

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; }

static inline bool _Stl_is_neg_nan(double x)    { return _isnan(x) && _copysign(1., x) < 0 ; }

#  elif defined (__MRC__) || defined (__SC__)    //*TY 02/24/2000 - added support for MPW

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)

static inline char* _Stl_qecvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return ecvtbuf(x, n, pt, sign, buf); }

static inline char* _Stl_qfcvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return fcvtbuf(x, n, pt, sign, buf); }

#    endif

#  elif defined (_STLP_USE_GLIBC)

static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf)

{ return buf + ecvt_r(x, n, pt, sign, buf, NDIG+2); }

static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf)

{ return buf + fcvt_r(x, n, pt, sign, buf, NDIG+2); }

#    if !defined (_STLP_NO_LONG_DOUBLE)

static inline char* _Stl_qecvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return buf + qecvt_r(x, n, pt, sign, buf, NDIG+2); }

static inline char* _Stl_qfcvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return buf + qfcvt_r(x, n, pt, sign, buf, NDIG+2); }

#    endif

#  elif defined (_STLP_SCO_OPENSERVER) || defined (__NCR_SVR)

static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf)

{ return ecvt(x, n, pt, sign); }

static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf)

{ return fcvt(x, n, pt, sign); }

#    if !defined (_STLP_NO_LONG_DOUBLE)

static inline char* _Stl_qecvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return ecvtl(x, n, pt, sign); }

static inline char* _Stl_qfcvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return fcvtl(x, n, pt, sign); }

#    endif

#  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); }

#    if !defined (_STLP_NO_LONG_DOUBLE)

static inline char* _Stl_qecvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return qeconvert(&x, n, pt, sign, buf); }

static inline char* _Stl_qfcvtR(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)

{ return (ecvt_r(x, n, pt, sign, buf, NDIG)==0 ? buf : 0); }

static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf)

{ return (fcvt_r(x, n, pt, sign, buf, NDIG)==0 ? buf : 0); }

#    if !defined (_STLP_NO_LONG_DOUBLE)

// fbp : no "long double" conversions !

static inline char* _Stl_qecvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return (ecvt_r((double)x, n, pt, sign, buf, NDIG)==0 ? buf : 0) ; }

static inline char* _Stl_qfcvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return (fcvt_r((double)x, n, pt, sign, buf, NDIG)==0 ? buf : 0); }

#    endif

#  elif defined (__hpux)

static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf)

{ return ecvt(x, n, pt, sign); }

static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf)

{ return fcvt(x, n, pt, sign); }

#    if !defined (_STLP_NO_LONG_DOUBLE)

#      if defined( _REENTRANT ) && (defined(_PTHREADS_DRAFT4) || defined(PTHREAD_THREADS_MAX))

static inline char* _Stl_qecvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return (_ldecvt_r(*(long_double*)&x, n, pt, sign, buf, NDIG+2)==0 ? buf : 0); }

static inline char* _Stl_qfcvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return (_ldfcvt_r(*(long_double*)&x, n, pt, sign, buf, NDIG+2)==0 ? buf : 0); }

#      else

static inline char* _Stl_qecvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return _ldecvt(*(long_double*)&x, n, pt, sign); }

static inline char* _Stl_qfcvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return _ldfcvt(*(long_double*)&x, n, pt, sign); }

#      endif

#    endif

#  elif defined (_AIX) || defined (__FreeBSD__) || defined (__NetBSD__) || defined (__OpenBSD__)

static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* buf)

{ LOCK_CVT RETURN_CVT(ecvt, x, n, pt, sign, buf) }

static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* buf)

{ LOCK_CVT RETURN_CVT(fcvt, x, n, pt, sign, buf) }

#    if !defined (_STLP_NO_LONG_DOUBLE)

static inline char* _Stl_qecvtR(long double x, int n, int* pt, int* sign, char* buf)

{ LOCK_CVT RETURN_CVT(ecvt, x, n, pt, sign, buf) }

static inline char* _Stl_qfcvtR(long double x, int n, int* pt, int* sign, char* buf)

{ LOCK_CVT RETURN_CVT(fcvt, x, n, pt, sign, buf) }

#    endif

#  elif defined (__unix) && !defined (__APPLE__) && !defined (_CRAY)

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_qecvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return qecvt_r(x, n, pt, sign, buf); }

static inline char* _Stl_qfcvtR(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__)

// those guys claim _cvt functions being reentrant.

#    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_CVT_DONT_NEED_BUF

#      define _STLP_BUF_PARAMS

#      define _STLP_SECURE_FUN(F, X, N, PT, SIGN) return F(X, N, PT, SIGN)

#    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)

static inline char* _Stl_qecvtR(long double x, int n, int* pt, int* sign _STLP_BUF_PARAMS)

{ _STLP_SECURE_FUN(_ecvt, (double)x, n, pt, sign); }

static inline char* _Stl_qfcvtR(long double x, int n, int* pt, int* sign _STLP_BUF_PARAMS)

{ _STLP_SECURE_FUN(_fcvt, (double)x, n, pt, sign); }

#    endif

#    undef _STLP_SECURE_FUN

#    undef _STLP_BUF_PARAMS

#    undef _STLP_APPEND

#  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_qecvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return _fp_ecvt( x, n, pt, sign, buf); }

static inline char* _Stl_qfcvtR(long double x, int n, int* pt, int* sign, char* buf)

{ return _fp_fcvt(x, n, pt, sign, buf); }

#    endif

#  elif defined (__MRC__) || defined (__SC__) || defined (_CRAY)

static inline char* _Stl_ecvtR(double x, int n, int* pt, int* sign, char* )

{ return ecvt( x, n, pt, sign ); }

static inline char* _Stl_fcvtR(double x, int n, int* pt, int* sign, char* )

{ return fcvt(x, n, pt, sign); }

#    if !defined (_STLP_NO_LONG_DOUBLE)

static inline char* _Stl_qecvtR(long double x, int n, int* pt, int* sign, char* )

{ return ecvt( x, n, pt, sign ); }

static inline char* _Stl_qfcvtR(long double x, int n, int* pt, int* sign, char* )

{ return fcvt(x, n, pt, sign); }

#    endif

#  endif

#  if defined (_STLP_CVT_DONT_NEED_BUF)

#    define _STLP_CVT_BUFFER(B)

#  elif !defined (_STLP_USE_SAFE_STRING_FUNCTIONS)

#    define _STLP_CVT_BUFFER(B) , B

#  else

#    define _STLP_CVT_BUFFER(B) , _STLP_ARRAY_AND_SIZE(B)

#  endif

#  if !defined (_STLP_USE_SAFE_STRING_FUNCTIONS)

#    define _STLP_BUFFER(B) B

#  else

#    define _STLP_BUFFER(B) _STLP_ARRAY_AND_SIZE(B)

#  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, bool /* islong */)

{

  // 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

  int rz = precision;

  while (rz-- > 0 && *bp != 0)

    buf += *bp++;

  // exponent

  char expbuf[MAXESIZ + 2];

  char *suffix = expbuf + MAXESIZ;

  *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

  while (suffix > &expbuf[MAXESIZ - 2])

    *--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, bool /* x */,

                                    ios_base::fmtflags flags,

                                    int precision, bool islong )

{

  if ( sign && (decpt > -precision) && (*bp != 0) )

    buf += '-';

  else if ( flags & ios_base::showpos )

    buf += '+';

  int k       = 0;

  int maxfsig = islong ? 2*MAXFSIG : MAXFSIG;

  // digits before decimal point

  int nnn = decpt;

  do {

    buf += ((nnn <= 0 || *bp == 0 || k >= maxfsig) ? '0' : (++k, *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

  nnn = (min) (precision, MAXFCVT);

  while ( --nnn >= 0 ) {

    buf += (++decpt <= 0 || *bp == 0 || k >= maxfsig) ? '0' : (++k, *bp++);

  }

  // trailing zeros if needed

  if ( precision > MAXFCVT ) {

    buf.append( precision - MAXFCVT, '0' );

  }

  return __group_pos;

}

static void __format_nan_or_inf(__iostring& buf, double 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 (_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 += '+';

  }

  buf += inf_or_nan[flags & ios_base::uppercase ? 1 : 0];

}

template <class max_double_type>

static inline size_t __format_float( __iostring &buf, const char * bp,

                                     int decpt, int sign, max_double_type x,

                                     ios_base::fmtflags flags,

                                     int precision, bool islong)

{

  size_t __group_pos = 0;

  // Output of infinities and NANs does not depend on the format flags

  if (_Stl_is_nan_or_inf((double)x)) {       // Infinity or NaN

    __format_nan_or_inf(buf, (double)x, flags);

  } else {                        // representable number

    switch (flags & ios_base::floatfield) {

      case ios_base::scientific:

        __group_pos = __format_float_scientific( buf, bp, decpt, sign, x == 0.0,

                                                 flags, precision, islong);

        break;

      case ios_base::fixed:

        __group_pos = __format_float_fixed( buf, bp, decpt, sign, true,

                                            flags, precision, islong);

        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 (x == 0)

          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, x == 0,

                                                   flags, precision, islong);

        } else {

          precision = kk - decpt;

          __group_pos = __format_float_fixed( buf, bp, decpt, sign, true,

                                              flags, precision, islong);

        }

        break;

    } /* switch */

  } /* else is_nan_or_inf */

  return __group_pos;

}

#else /* USE_SPRINTF_INSTEAD */

struct GroupPos {

  bool operator () (char __c) const {

    return __c == '.' ||

           __c == 'e' || __c == 'E';

  }

};

// Creates a format string for sprintf()

static int __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;

  return i;

}

#endif  /* USE_SPRINTF_INSTEAD */

_STLP_DECLSPEC size_t  _STLP_CALL

__write_float(__iostring &buf, ios_base::fmtflags flags, int precision,

              double x) {

#if defined (USE_SPRINTF_INSTEAD)

  /* If we want 'abitrary' precision, we should use 'abitrary' buffer size

   * below. - ptr

   */

  char static_buf[128];

  // char *static_buf = new char [128+precision];

  char fmtbuf[32];

  __fill_fmtbuf(fmtbuf, flags, 0);

  // snprintf(static_buf, 128+precision, fmtbuf, precision, x);

#  if !defined (N_PLAT_NLM)

  snprintf(_STLP_ARRAY_AND_SIZE(static_buf), fmtbuf, precision, x);

#  else

  sprintf(static_buf, fmtbuf, precision, x);

#  endif

  buf = static_buf;

  // delete [] static_buf;

  return find_if(buf.begin(), buf.end(), GroupPos()) - buf.begin();

#else

#  if !defined (_STLP_CVT_DONT_NEED_BUF)

  char cvtbuf[NDIG + 2];

#  endif

  char * bp;

  int decpt, sign;

  switch (flags & ios_base::floatfield) {

  case ios_base::fixed:

    bp = _Stl_fcvtR(x, (min) (precision, MAXFCVT), &decpt, &sign _STLP_CVT_BUFFER(cvtbuf));

    break;

  case ios_base::scientific :

    bp = _Stl_ecvtR(x, (min) (precision + 1, MAXECVT), &decpt, &sign _STLP_CVT_BUFFER(cvtbuf));

    break;

  default :

    bp = _Stl_ecvtR(x, (min) (precision, MAXECVT), &decpt, &sign _STLP_CVT_BUFFER(cvtbuf));

    break;

  }

  return __format_float(buf, bp, decpt, sign, x, flags, precision, false);

#endif

}

#if !defined (_STLP_NO_LONG_DOUBLE)

_STLP_DECLSPEC size_t _STLP_CALL

__write_float(__iostring &buf, ios_base::fmtflags flags, int precision,

              long double x) {

#  if defined (USE_SPRINTF_INSTEAD)

  /* If we want 'abitrary' precision, we should use 'abitrary' buffer size

   * below. - ptr

   */

  char static_buf[128];

  // char *static_buf = new char [128+precision];

  char fmtbuf[64];

  int i = __fill_fmtbuf(fmtbuf, flags, 'L');

  // snprintf(static_buf, 128+precision, fmtbuf, precision, x);

#    if !defined (N_PLAT_NLM)

  snprintf(_STLP_ARRAY_AND_SIZE(static_buf), fmtbuf, precision, x);

#    else

  sprintf(static_buf, fmtbuf, precision, x);

#    endif

  // we should be able to return buf + sprintf(), but we do not trust'em...

  buf = static_buf;

  // delete [] static_buf;

  return find_if(buf.begin(), buf.end(), GroupPos()) - buf.begin();

#  else

#    if !defined (_STLP_CVT_DONT_NEED_BUF)

  char cvtbuf[NDIG + 2];

#    endif

  char * bp;

  int decpt, sign;

  switch (flags & ios_base::floatfield) {

  case ios_base::fixed:

    bp = _Stl_qfcvtR(x, (min) (precision, MAXFCVT), &decpt, &sign _STLP_CVT_BUFFER(cvtbuf));

    break;

  case ios_base::scientific:

    bp = _Stl_qecvtR(x, (min) (precision + 1, MAXECVT), &decpt, &sign _STLP_CVT_BUFFER(cvtbuf));

    break;

  default :

    bp = _Stl_qecvtR(x, (min) (precision, MAXECVT), &decpt, &sign _STLP_CVT_BUFFER(cvtbuf));

    break;

  }

  return __format_float(buf, bp, decpt, sign, x, flags, precision, true);

#  endif /* USE_SPRINTF_INSTEAD */

}

#endif /* _STLP_NO_LONG_DOUBLE */

_STLP_DECLSPEC void _STLP_CALL __get_floor_digits(__iostring &out, _STLP_LONGEST_FLOAT_TYPE __x) {

#if defined (USE_SPRINTF_INSTEAD)

  char cvtbuf[128];

#  if !defined (_STLP_NO_LONG_DOUBLE)

#    if !defined (N_PLAT_NLM)

  snprintf(_STLP_ARRAY_AND_SIZE(cvtbuf), "%Lf", __x); // check for 1234.56!

#    else

  sprintf(cvtbuf, "%Lf", __x); // check for 1234.56!

#    endif

#  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

#  if !defined (_STLP_CVT_DONT_NEED_BUF)

  char cvtbuf[NDIG + 2];

#  endif

  char * bp;

  int decpt, sign;

#  if !defined (_STLP_NO_LONG_DOUBLE)

  bp = _Stl_qfcvtR(__x, 0, &decpt, &sign _STLP_CVT_BUFFER(cvtbuf));

#  else

  bp = _Stl_fcvtR(__x, 0, &decpt, &sign _STLP_CVT_BUFFER(cvtbuf));

#  endif

  if (sign) {

    out += '-';

  }

  out.append(bp, bp + decpt);

#endif // USE_PRINTF_INSTEAD

}

#if !defined (_STLP_NO_WCHAR_T)

_STLP_DECLSPEC void _STLP_CALL __convert_float_buffer( __iostring const& str, __iowstring &out,

                                        const ctype<wchar_t>& 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: