Put a @table around @include summary.out.
[kopensolaris-gnu/glibc.git] / stdlib / strtod.c
1 /* Copyright (C) 1991 Free Software Foundation, Inc.
2 This file is part of the GNU C Library.
3
4 The GNU C Library is free software; you can redistribute it and/or
5 modify it under the terms of the GNU Library General Public License as
6 published by the Free Software Foundation; either version 2 of the
7 License, or (at your option) any later version.
8
9 The GNU C Library is distributed in the hope that it will be useful,
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
12 Library General Public License for more details.
13
14 You should have received a copy of the GNU Library General Public
15 License along with the GNU C Library; see the file COPYING.LIB.  If
16 not, write to the Free Software Foundation, Inc., 675 Mass Ave,
17 Cambridge, MA 02139, USA.  */
18
19 #include <ansidecl.h>
20 #include <localeinfo.h>
21 #include <errno.h>
22 #include <float.h>
23 #include <ctype.h>
24 #include <stdlib.h>
25 #include <string.h>
26 #include <math.h>
27
28
29 /* Convert NPTR to a double.  If ENDPTR is not NULL, a pointer to the
30    character after the last one used in the number is put in *ENDPTR.  */
31 double
32 DEFUN(strtod, (nptr, endptr), CONST char *nptr AND char **endptr)
33 {
34   register CONST char *s;
35   short int sign;
36   wchar_t decimal;      /* Decimal point character.  */
37
38   /* The number so far.  */
39   double num;
40
41   int got_dot;          /* Found a decimal point.  */
42   int got_digit;        /* Seen any digits.  */
43   int nonzero;          /* Seen any nonzero digits.  */
44
45   /* The exponent of the number.  */
46   long int exponent;
47
48   if (nptr == NULL)
49     {
50       errno = EINVAL;
51       goto noconv;
52     }
53
54   /* Figure out the decimal point character.  */
55   if (mbtowc(&decimal, _numeric_info->decimal_point, 1) <= 0)
56     decimal = (wchar_t) *_numeric_info->decimal_point;
57
58   s = nptr;
59
60   /* Eat whitespace.  */
61   while (isspace(*s))
62     ++s;
63
64   /* Get the sign.  */
65   sign = *s == '-' ? -1 : 1;
66   if (*s == '-' || *s == '+')
67     ++s;
68
69   num = 0.0;
70   got_dot = 0;
71   got_digit = 0;
72   nonzero = 0;
73   exponent = 0;
74   for (;; ++s)
75     {
76       if (isdigit (*s))
77         {
78           got_digit = 1;
79
80           if (*s != '0')
81             nonzero = 1;
82
83           /* Make sure that multiplication by 10 will not overflow.  */
84           if (num > DBL_MAX * 0.1)
85             /* The value of the digit doesn't matter, since we have already
86                gotten as many digits as can be represented in a `double'.
87                This doesn't necessarily mean the result will overflow.
88                The exponent may reduce it to within range.
89
90                We just need to record that there was another
91                digit so that we can multiply by 10 later.  */
92             ++exponent;
93           else
94             num = (num * 10.0) + (*s - '0');
95
96           /* Keep track of the number of digits after the decimal point.
97              If we just divided by 10 here, we would lose precision.  */
98           if (got_dot)
99             --exponent;
100         }
101       else if (!got_dot && (wchar_t) *s == decimal)
102         /* Record that we have found the decimal point.  */
103         got_dot = 1;
104       else
105         /* Any other character terminates the number.  */
106         break;
107     }
108
109   if (!got_digit)
110     goto noconv;
111
112   if (tolower(*s) == 'e')
113     {
114       /* Get the exponent specified after the `e' or `E'.  */
115       int save = errno;
116       char *end;
117       long int exp;
118
119       errno = 0;
120       ++s;
121       exp = strtol(s, &end, 10);
122       if (end == s)
123         {
124           if (errno == ERANGE)
125             {
126               /* The exponent overflowed a `long int'.  It is probably a safe
127                  assumption that an exponent that cannot be represented by
128                  a `long int' exceeds the limits of a `double'.  */
129               if (exp < 0)
130                 goto underflow;
131               else
132                 goto overflow;
133             }
134           else
135             /* There was no exponent.  Reset END to point to
136                the 'e' or 'E', so *ENDPTR will be set there.  */
137             end = (char *) s - 1;
138         }
139       errno = save;
140       s = end;
141       exponent += exp;
142     }
143
144   if (num == 0.0)
145     return 0.0;
146
147   /* Multiply NUM by 10 to the EXPONENT power,
148      checking for overflow and underflow.  */
149
150   if (exponent < 0)
151     {
152       if (num < DBL_MIN * pow(10.0, (double) -exponent))
153         goto underflow;
154     }
155   else if (exponent > 0)
156     {
157       if (num > DBL_MAX * pow(10.0, (double) exponent))
158         goto overflow;
159     }
160
161   num *= pow(10.0, (double) exponent);
162
163   /* Set *ENDPTR and return the result.  */
164   if (endptr != NULL)
165     *endptr = (char *) s;
166   return num * sign;
167
168 overflow:;
169   /* Return an overflow error.  */
170   if (endptr != NULL)
171     *endptr = (char *) nptr;
172   errno = ERANGE;
173   return HUGE_VAL * sign;
174
175 underflow:;
176   /* Return an underflow error.  */
177   if (endptr != NULL)
178     *endptr = (char *) nptr;
179   errno = ERANGE;
180   return 0.0;
181
182 noconv:;
183   /* There was no number.  */
184   if (endptr != NULL)
185     *endptr = (char *) nptr;
186   return 0.0;
187 }