Use MAP_ANON instead of MAP_ANONYMOUS.
[kopensolaris-gnu/glibc.git] / misc / efgcvt_r.c
index 2dd0c9f..1c237ee 100644 (file)
@@ -23,6 +23,7 @@
 #include <ctype.h>
 #include <math.h>
 #include <stdlib.h>
+#include <sys/param.h>
 
 #ifndef FLOAT_TYPE
 #define FLOAT_TYPE double
@@ -48,6 +49,7 @@ APPEND (FUNC_PREFIX, fcvt_r) (value, ndigit, decpt, sign, buf, len)
      size_t len;
 {
   int n, i;
+  int left;
 
   if (buf == NULL)
     {
@@ -55,25 +57,37 @@ APPEND (FUNC_PREFIX, fcvt_r) (value, ndigit, decpt, sign, buf, len)
       return -1;
     }
 
+  left = 0;
   if (isfinite (value))
     {
       *sign = signbit (value) != 0;
       if (*sign)
        value = -value;
+
+      if (ndigit < 0)
+       {
+         /* Rounding to the left of the decimal point.  */
+         while (ndigit < 0)
+           {
+             FLOAT_TYPE new_value = value * 0.1;
+
+             if (new_value < 1.0)
+               {
+                 ndigit = 0;
+                 break;
+               }
+
+             value = new_value;
+             ++left;
+             ++ndigit;
+           }
+       }
     }
   else
     /* Value is Inf or NaN.  */
     *sign = 0;
 
-  if (ndigit <= 0)
-    {
-      if (len > 0)
-       buf[0] = '\0';
-      *decpt = 0;
-      return 0;
-    }
-
-  n = snprintf (buf, len, "%.*" FLOAT_FMT_FLAG "f", ndigit, value);
+  n = __snprintf (buf, len, "%.*" FLOAT_FMT_FLAG "f", ndigit, value);
   if (n < 0)
     return -1;
 
@@ -91,17 +105,37 @@ APPEND (FUNC_PREFIX, fcvt_r) (value, ndigit, decpt, sign, buf, len)
       do
        ++i;
       while (i < n && !isdigit (buf[i]));
-      memmove (&buf[*decpt], &buf[i], n - i);
-      buf[n - (i - *decpt)] = '\0';
+
+      if (*decpt == 1 && buf[0] == '0' && value != 0.0)
+       {
+         /* We must not have leading zeroes.  Strip them all out and
+            adjust *DECPT if necessary.  */
+         --*decpt;
+         while (i < n && buf[i] == '0')
+           {
+             --*decpt;
+             ++i;
+           }
+       }
+
+      memmove (&buf[MAX (*decpt, 0)], &buf[i], n - i);
+      buf[n - (i - MAX (*decpt, 0))] = '\0';
+    }
+
+  if (left)
+    {
+      *decpt += left;
+      if (--len > n)
+       {
+         while (left-- > 0 && n < len)
+           buf[n++] = '0';
+         buf[n] = '\0';
+       }
     }
 
   return 0;
 }
 
-#define weak_extern2(name) weak_extern (name)
-weak_extern2 (FLOOR) weak_extern2 (LOG10) weak_extern2 (FABS)
-weak_extern2 (EXP)
-
 int
 APPEND (FUNC_PREFIX, ecvt_r) (value, ndigit, decpt, sign, buf, len)
      FLOAT_TYPE value;
@@ -113,51 +147,52 @@ APPEND (FUNC_PREFIX, ecvt_r) (value, ndigit, decpt, sign, buf, len)
 
   if (isfinite (value) && value != 0.0)
     {
-      FLOAT_TYPE (*log10_function) (FLOAT_TYPE) = &LOG10;
-
-      if (log10_function)
-       {
-         /* Use the reasonable code if -lm is included.  */
-         FLOAT_TYPE dexponent;
-         dexponent = FLOOR (LOG10 (FABS (value)));
-         value *= EXP (dexponent * -M_LN10);
-         exponent = (int) dexponent;
-       }
+      /* Slow code that doesn't require -lm functions.  */
+      FLOAT_TYPE d;
+      FLOAT_TYPE f = 1.0;
+      if (value < 0.0)
+       d = -value;
       else
+       d = value;
+      if (d < 1.0)
        {
-         /* Slow code that doesn't require -lm functions.  */
-         FLOAT_TYPE d;
-         if (value < 0.0)
-           d = -value;
-         else
-           d = value;
-         if (d < 1.0)
+         do
            {
-             do
-               {
-                 d *= 10.0;
-                 exponent--;
-               }
-             while (d < 1.0);
+             f *= 10.0;
+             --exponent;
            }
-         else if (d >= 10.0)
+         while (d * f < 1.0);
+
+         value *= f;
+       }
+      else if (d >= 10.0)
+       {
+         do
            {
-             do
-               {
-                 d *= 0.1;
-                 exponent++;
-               }
-             while (d >= 10.0);
+             f *= 10;
+             ++exponent;
            }
-         if (value < 0.0)
-           value = -d;
-         else
-           value = d;
+         while (d >= f * 10.0);
+
+         value /= f;
        }
     }
+  else if (value == 0.0)
+    /* SUSv2 leaves it unspecified whether *DECPT is 0 or 1 for 0.0.
+       This could be changed to -1 if we want to return 0.  */
+    exponent = 0;
+
+  if (ndigit <= 0 && len > 0)
+    {
+      buf[0] = '\0';
+      *decpt = 1;
+      *sign = isfinite (value) ? signbit (value) != 0 : 0;
+    }
+  else
+    if (APPEND (FUNC_PREFIX, fcvt_r) (value, ndigit - 1, decpt, sign,
+                                     buf, len))
+      return -1;
 
-  if (APPEND (FUNC_PREFIX, fcvt_r) (value, ndigit - 1, decpt, sign, buf, len))
-    return -1;
   *decpt += exponent;
   return 0;
 }