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