strtod.c 16 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614
  1. /*
  2. * Copyright (C) 2000-2005 Manuel Novoa III
  3. *
  4. * Licensed under the LGPL v2.1, see the file COPYING.LIB in this tarball.
  5. */
  6. /* Notes:
  7. *
  8. * The primary objective of this implementation was minimal size and
  9. * portablility, while providing robustness and resonable accuracy.
  10. *
  11. * This implementation depends on IEEE floating point behavior and expects
  12. * to be able to generate +/- infinity as a result.
  13. *
  14. * There are a number of compile-time options below.
  15. */
  16. /* July 27, 2003
  17. *
  18. * General cleanup and some minor size optimizations.
  19. * Change implementation to support __strtofpmax() rather than strtod().
  20. * Now all the strto{floating pt}() funcs are implemented in terms of
  21. * of the internal __strtofpmax() function.
  22. * Support "nan", "inf", and "infinity" strings (case-insensitive).
  23. * Support hexadecimal floating point notation.
  24. * Support wchar variants.
  25. * Support xlocale variants.
  26. *
  27. * TODO:
  28. *
  29. * Consider accumulating blocks of digits in longs to save floating pt mults.
  30. * This would likely be much better on anything that only supported floats
  31. * where DECIMAL_DIG == 9. Actually, if floats have FLT_MAX_10_EXP == 38,
  32. * we could calculate almost all the exponent multipliers (p_base) in
  33. * long arithmetic as well.
  34. */
  35. /**********************************************************************/
  36. /* OPTIONS */
  37. /**********************************************************************/
  38. /* Defined if we want to recognize "nan", "inf", and "infinity". (C99) */
  39. #define _STRTOD_NAN_INF_STRINGS 1
  40. /* Defined if we want support hexadecimal floating point notation. (C99) */
  41. /* Note! Now controlled by uClibc configuration. See below. */
  42. #define _STRTOD_HEXADECIMAL_FLOATS 1
  43. /* Defined if we want to scale with a O(log2(exp)) multiplications.
  44. * This is generally a good thing to do unless you are really tight
  45. * on space and do not expect to convert values of large magnitude. */
  46. #define _STRTOD_LOG_SCALING 1
  47. /* WARNING!!! WARNING!!! WARNING!!! WARNING!!! WARNING!!!
  48. *
  49. * Clearing any of the options below this point is not advised (or tested).
  50. *
  51. * WARNING!!! WARNING!!! WARNING!!! WARNING!!! WARNING!!! */
  52. /* Defined if we want strtod to set errno appropriately. */
  53. /* NOTE: Implies all options below. */
  54. #define _STRTOD_ERRNO 1
  55. /* Defined if we want support for the endptr arg. */
  56. /* Implied by _STRTOD_ERRNO. */
  57. #define _STRTOD_ENDPTR 1
  58. /* Defined if we want to prevent overflow in accumulating the exponent. */
  59. /* Implied by _STRTOD_ERRNO. */
  60. #define _STRTOD_RESTRICT_EXP 1
  61. /* Defined if we want to process mantissa digits more intelligently. */
  62. /* Implied by _STRTOD_ERRNO. */
  63. #define _STRTOD_RESTRICT_DIGITS 1
  64. /* Defined if we want to skip scaling 0 for the exponent. */
  65. /* Implied by _STRTOD_ERRNO. */
  66. #define _STRTOD_ZERO_CHECK 1
  67. /**********************************************************************/
  68. /* Don't change anything that follows. */
  69. /**********************************************************************/
  70. #ifdef _STRTOD_ERRNO
  71. #undef _STRTOD_ENDPTR
  72. #undef _STRTOD_RESTRICT_EXP
  73. #undef _STRTOD_RESTRICT_DIGITS
  74. #undef _STRTOD_ZERO_CHECK
  75. #define _STRTOD_ENDPTR 1
  76. #define _STRTOD_RESTRICT_EXP 1
  77. #define _STRTOD_RESTRICT_DIGITS 1
  78. #define _STRTOD_ZERO_CHECK 1
  79. #endif
  80. /**********************************************************************/
  81. #define _ISOC99_SOURCE 1
  82. #define _GNU_SOURCE
  83. #include <stdlib.h>
  84. #include <string.h>
  85. #include <ctype.h>
  86. #include <errno.h>
  87. #include <limits.h>
  88. #include <float.h>
  89. #include <bits/uClibc_fpmax.h>
  90. #include <locale.h>
  91. #ifdef __UCLIBC_HAS_WCHAR__
  92. #include <wchar.h>
  93. #include <wctype.h>
  94. #include <bits/uClibc_uwchar.h>
  95. libc_hidden_proto(iswspace)
  96. #endif
  97. #ifdef __UCLIBC_HAS_XLOCALE__
  98. #include <xlocale.h>
  99. libc_hidden_proto(iswspace_l)
  100. #endif /* __UCLIBC_HAS_XLOCALE__ */
  101. /* Handle _STRTOD_HEXADECIMAL_FLOATS via uClibc config now. */
  102. #undef _STRTOD_HEXADECIMAL_FLOATS
  103. #ifdef __UCLIBC_HAS_HEXADECIMAL_FLOATS__
  104. #define _STRTOD_HEXADECIMAL_FLOATS 1
  105. #endif /* __UCLIBC_HAS_HEXADECIMAL_FLOATS__ */
  106. /**********************************************************************/
  107. #undef _STRTOD_FPMAX
  108. #if FPMAX_TYPE == 3
  109. #define NEED_STRTOLD_WRAPPER
  110. #define NEED_STRTOD_WRAPPER
  111. #define NEED_STRTOF_WRAPPER
  112. #elif FPMAX_TYPE == 2
  113. #define NEED_STRTOD_WRAPPER
  114. #define NEED_STRTOF_WRAPPER
  115. #elif FPMAX_TYPE == 1
  116. #define NEED_STRTOF_WRAPPER
  117. #else
  118. #error unknown FPMAX_TYPE!
  119. #endif
  120. extern void __fp_range_check(__fpmax_t y, __fpmax_t x) attribute_hidden;
  121. /**********************************************************************/
  122. #ifdef _STRTOD_RESTRICT_DIGITS
  123. #define EXP_DENORM_ADJUST DECIMAL_DIG
  124. #define MAX_ALLOWED_EXP (DECIMAL_DIG + EXP_DENORM_ADJUST - FPMAX_MIN_10_EXP)
  125. #if MAX_ALLOWED_EXP > INT_MAX
  126. #error size assumption violated for MAX_ALLOWED_EXP
  127. #endif
  128. #else
  129. /* We want some excess if we're not restricting mantissa digits. */
  130. #define MAX_ALLOWED_EXP ((20 - FPMAX_MIN_10_EXP) * 2)
  131. #endif
  132. #if defined(_STRTOD_RESTRICT_DIGITS) || defined(_STRTOD_ENDPTR) || defined(_STRTOD_HEXADECIMAL_FLOATS)
  133. #undef _STRTOD_NEED_NUM_DIGITS
  134. #define _STRTOD_NEED_NUM_DIGITS 1
  135. #endif
  136. /**********************************************************************/
  137. #if defined(L___strtofpmax) || defined(L___strtofpmax_l) || defined(L___wcstofpmax) || defined(L___wcstofpmax_l)
  138. #if defined(L___wcstofpmax) || defined(L___wcstofpmax_l)
  139. #define __strtofpmax __wcstofpmax
  140. #define __strtofpmax_l __wcstofpmax_l
  141. #define Wchar wchar_t
  142. #ifdef __UCLIBC_DO_XLOCALE
  143. #define ISSPACE(C) iswspace_l((C), locale_arg)
  144. #else
  145. #define ISSPACE(C) iswspace((C))
  146. #endif
  147. #else /* defined(L___wcstofpmax) || defined(L___wcstofpmax_l) */
  148. #define Wchar char
  149. #ifdef __UCLIBC_DO_XLOCALE
  150. #define ISSPACE(C) isspace_l((C), locale_arg)
  151. #else
  152. #define ISSPACE(C) isspace((C))
  153. #endif
  154. #endif /* defined(L___wcstofpmax) || defined(L___wcstofpmax_l) */
  155. #if defined(__UCLIBC_HAS_XLOCALE__) && !defined(__UCLIBC_DO_XLOCALE)
  156. __fpmax_t attribute_hidden __strtofpmax(const Wchar *str, Wchar **endptr, int exponent_power)
  157. {
  158. return __strtofpmax_l(str, endptr, exponent_power, __UCLIBC_CURLOCALE);
  159. }
  160. #else /* defined(__UCLIBC_HAS_XLOCALE__) && !defined(__UCLIBC_DO_XLOCALE) */
  161. libc_hidden_proto(memcmp)
  162. __fpmax_t attribute_hidden __XL_NPP(__strtofpmax)(const Wchar *str, Wchar **endptr, int exponent_power
  163. __LOCALE_PARAM )
  164. {
  165. __fpmax_t number;
  166. __fpmax_t p_base = 10; /* Adjusted to 16 in the hex case. */
  167. Wchar *pos0;
  168. #ifdef _STRTOD_ENDPTR
  169. Wchar *pos1;
  170. #endif
  171. Wchar *pos = (Wchar *) str;
  172. int exponent_temp;
  173. int negative; /* A flag for the number, a multiplier for the exponent. */
  174. #ifdef _STRTOD_NEED_NUM_DIGITS
  175. int num_digits;
  176. #endif
  177. #ifdef __UCLIBC_HAS_LOCALE__
  178. #if defined(L___wcstofpmax) || defined(L___wcstofpmax_l)
  179. wchar_t decpt_wc = __LOCALE_PTR->decimal_point;
  180. #else
  181. const char *decpt = __LOCALE_PTR->decimal_point;
  182. int decpt_len = __LOCALE_PTR->decimal_point_len;
  183. #endif
  184. #endif
  185. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  186. Wchar expchar = 'e';
  187. Wchar *poshex = NULL;
  188. __uint16_t is_mask = _ISdigit;
  189. #define EXPCHAR expchar
  190. #define IS_X_DIGIT(C) __isctype((C), is_mask)
  191. #else /* _STRTOD_HEXADECIMAL_FLOATS */
  192. #define EXPCHAR 'e'
  193. #define IS_X_DIGIT(C) isdigit((C))
  194. #endif /* _STRTOD_HEXADECIMAL_FLOATS */
  195. while (ISSPACE(*pos)) { /* Skip leading whitespace. */
  196. ++pos;
  197. }
  198. negative = 0;
  199. switch(*pos) { /* Handle optional sign. */
  200. case '-': negative = 1; /* Fall through to increment position. */
  201. case '+': ++pos;
  202. }
  203. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  204. if ((*pos == '0') && (((pos[1])|0x20) == 'x')) {
  205. poshex = ++pos; /* Save position of 'x' in case no digits */
  206. ++pos; /* and advance past it. */
  207. is_mask = _ISxdigit; /* Used by IS_X_DIGIT. */
  208. expchar = 'p'; /* Adjust exponent char. */
  209. p_base = 16; /* Adjust base multiplier. */
  210. }
  211. #endif
  212. number = 0.;
  213. #ifdef _STRTOD_NEED_NUM_DIGITS
  214. num_digits = -1;
  215. #endif
  216. /* exponent_power = 0; */
  217. pos0 = NULL;
  218. LOOP:
  219. while (IS_X_DIGIT(*pos)) { /* Process string of (hex) digits. */
  220. #ifdef _STRTOD_RESTRICT_DIGITS
  221. if (num_digits < 0) { /* First time through? */
  222. ++num_digits; /* We've now seen a digit. */
  223. }
  224. if (num_digits || (*pos != '0')) { /* Had/have nonzero. */
  225. ++num_digits;
  226. if (num_digits <= DECIMAL_DIG) { /* Is digit significant? */
  227. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  228. number = number * p_base
  229. + (isdigit(*pos)
  230. ? (*pos - '0')
  231. : (((*pos)|0x20) - ('a' - 10)));
  232. #else /* _STRTOD_HEXADECIMAL_FLOATS */
  233. number = number * p_base + (*pos - '0');
  234. #endif /* _STRTOD_HEXADECIMAL_FLOATS */
  235. }
  236. }
  237. #else /* _STRTOD_RESTRICT_DIGITS */
  238. #ifdef _STRTOD_NEED_NUM_DIGITS
  239. ++num_digits;
  240. #endif
  241. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  242. number = number * p_base
  243. + (isdigit(*pos)
  244. ? (*pos - '0')
  245. : (((*pos)|0x20) - ('a' - 10)));
  246. #else /* _STRTOD_HEXADECIMAL_FLOATS */
  247. number = number * p_base + (*pos - '0');
  248. #endif /* _STRTOD_HEXADECIMAL_FLOATS */
  249. #endif /* _STRTOD_RESTRICT_DIGITS */
  250. ++pos;
  251. }
  252. #ifdef __UCLIBC_HAS_LOCALE__
  253. #if defined(L___wcstofpmax) || defined(L___wcstofpmax_l)
  254. if (!pos0 && (*pos == decpt_wc)) { /* First decimal point? */
  255. pos0 = ++pos;
  256. goto LOOP;
  257. }
  258. #else
  259. if (!pos0 && !memcmp(pos, decpt, decpt_len)) { /* First decimal point? */
  260. pos0 = (pos += decpt_len);
  261. goto LOOP;
  262. }
  263. #endif
  264. #else /* __UCLIBC_HAS_LOCALE__ */
  265. if ((*pos == '.') && !pos0) { /* First decimal point? */
  266. pos0 = ++pos; /* Save position of decimal point */
  267. goto LOOP; /* and process rest of digits. */
  268. }
  269. #endif /* __UCLIBC_HAS_LOCALE__ */
  270. #ifdef _STRTOD_NEED_NUM_DIGITS
  271. if (num_digits<0) { /* Must have at least one digit. */
  272. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  273. if (poshex) { /* Back up to '0' in '0x' prefix. */
  274. pos = poshex;
  275. goto DONE;
  276. }
  277. #endif /* _STRTOD_HEXADECIMAL_FLOATS */
  278. #ifdef _STRTOD_NAN_INF_STRINGS
  279. if (!pos0) { /* No decimal point, so check for inf/nan. */
  280. /* Note: nan is the first string so 'number = i/0.;' works. */
  281. static const char nan_inf_str[] = "\05nan\0\012infinity\0\05inf\0";
  282. int i = 0;
  283. #ifdef __UCLIBC_HAS_LOCALE__
  284. /* Avoid tolower problems for INFINITY in the tr_TR locale. (yuk)*/
  285. #undef _tolower
  286. #define _tolower(C) ((C)|0x20)
  287. #endif /* __UCLIBC_HAS_LOCALE__ */
  288. do {
  289. /* Unfortunately, we have no memcasecmp(). */
  290. int j = 0;
  291. while (_tolower(pos[j]) == nan_inf_str[i+1+j]) {
  292. ++j;
  293. if (!nan_inf_str[i+1+j]) {
  294. number = i / 0.;
  295. if (negative) { /* Correct for sign. */
  296. number = -number;
  297. }
  298. pos += nan_inf_str[i] - 2;
  299. goto DONE;
  300. }
  301. }
  302. i += nan_inf_str[i];
  303. } while (nan_inf_str[i]);
  304. }
  305. #endif /* STRTOD_NAN_INF_STRINGS */
  306. #ifdef _STRTOD_ENDPTR
  307. pos = (Wchar *) str;
  308. #endif
  309. goto DONE;
  310. }
  311. #endif /* _STRTOD_NEED_NUM_DIGITS */
  312. #ifdef _STRTOD_RESTRICT_DIGITS
  313. if (num_digits > DECIMAL_DIG) { /* Adjust exponent for skipped digits. */
  314. exponent_power += num_digits - DECIMAL_DIG;
  315. }
  316. #endif
  317. if (pos0) {
  318. exponent_power += pos0 - pos; /* Adjust exponent for decimal point. */
  319. }
  320. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  321. if (poshex) {
  322. exponent_power *= 4; /* Above is 2**4, but below is 2. */
  323. p_base = 2;
  324. }
  325. #endif /* _STRTOD_HEXADECIMAL_FLOATS */
  326. if (negative) { /* Correct for sign. */
  327. number = -number;
  328. }
  329. /* process an exponent string */
  330. if (((*pos)|0x20) == EXPCHAR) {
  331. #ifdef _STRTOD_ENDPTR
  332. pos1 = pos;
  333. #endif
  334. negative = 1;
  335. switch(*++pos) { /* Handle optional sign. */
  336. case '-': negative = -1; /* Fall through to increment pos. */
  337. case '+': ++pos;
  338. }
  339. pos0 = pos;
  340. exponent_temp = 0;
  341. while (isdigit(*pos)) { /* Process string of digits. */
  342. #ifdef _STRTOD_RESTRICT_EXP
  343. if (exponent_temp < MAX_ALLOWED_EXP) { /* Avoid overflow. */
  344. exponent_temp = exponent_temp * 10 + (*pos - '0');
  345. }
  346. #else
  347. exponent_temp = exponent_temp * 10 + (*pos - '0');
  348. #endif
  349. ++pos;
  350. }
  351. #ifdef _STRTOD_ENDPTR
  352. if (pos == pos0) { /* No digits? */
  353. pos = pos1; /* Back up to {e|E}/{p|P}. */
  354. } /* else */
  355. #endif
  356. exponent_power += negative * exponent_temp;
  357. }
  358. #ifdef _STRTOD_ZERO_CHECK
  359. if (number == 0.) {
  360. goto DONE;
  361. }
  362. #endif
  363. /* scale the result */
  364. #ifdef _STRTOD_LOG_SCALING
  365. exponent_temp = exponent_power;
  366. if (exponent_temp < 0) {
  367. exponent_temp = -exponent_temp;
  368. }
  369. while (exponent_temp) {
  370. if (exponent_temp & 1) {
  371. if (exponent_power < 0) {
  372. /* Warning... caluclating a factor for the exponent and
  373. * then dividing could easily be faster. But doing so
  374. * might cause problems when dealing with denormals. */
  375. number /= p_base;
  376. } else {
  377. number *= p_base;
  378. }
  379. }
  380. exponent_temp >>= 1;
  381. p_base *= p_base;
  382. }
  383. #else /* _STRTOD_LOG_SCALING */
  384. while (exponent_power) {
  385. if (exponent_power < 0) {
  386. number /= p_base;
  387. exponent_power++;
  388. } else {
  389. number *= p_base;
  390. exponent_power--;
  391. }
  392. }
  393. #endif /* _STRTOD_LOG_SCALING */
  394. #ifdef _STRTOD_ERRNO
  395. if (__FPMAX_ZERO_OR_INF_CHECK(number)) {
  396. __set_errno(ERANGE);
  397. }
  398. #endif
  399. DONE:
  400. #ifdef _STRTOD_ENDPTR
  401. if (endptr) {
  402. *endptr = pos;
  403. }
  404. #endif
  405. return number;
  406. }
  407. #endif /* defined(__UCLIBC_HAS_XLOCALE__) && !defined(__UCLIBC_DO_XLOCALE) */
  408. #endif
  409. /**********************************************************************/
  410. #ifdef L___fp_range_check
  411. #if defined(NEED_STRTOF_WRAPPER) || defined(NEED_STRTOD_WRAPPER)
  412. void attribute_hidden __fp_range_check(__fpmax_t y, __fpmax_t x)
  413. {
  414. if (__FPMAX_ZERO_OR_INF_CHECK(y) /* y is 0 or +/- infinity */
  415. && (y != 0) /* y is not 0 (could have x>0, y==0 if underflow) */
  416. && !__FPMAX_ZERO_OR_INF_CHECK(x) /* x is not 0 or +/- infinity */
  417. ) {
  418. __set_errno(ERANGE); /* Then x is not in y's range. */
  419. }
  420. }
  421. #endif
  422. #endif
  423. /**********************************************************************/
  424. #if defined(L_strtof) || defined(L_strtof_l) || defined(L_wcstof) || defined(L_wcstof_l)
  425. #if defined(NEED_STRTOF_WRAPPER)
  426. #if defined(L_wcstof) || defined(L_wcstof_l)
  427. #define strtof wcstof
  428. #define strtof_l wcstof_l
  429. #define __strtofpmax __wcstofpmax
  430. #define __strtofpmax_l __wcstofpmax_l
  431. #define Wchar wchar_t
  432. #else
  433. #define Wchar char
  434. #endif
  435. libc_hidden_proto(__XL_NPP(strtof))
  436. float __XL_NPP(strtof)(const Wchar *str, Wchar **endptr __LOCALE_PARAM )
  437. {
  438. #if FPMAX_TYPE == 1
  439. return __XL_NPP(__strtofpmax)(str, endptr, 0 __LOCALE_ARG );
  440. #else
  441. __fpmax_t x;
  442. float y;
  443. x = __XL_NPP(__strtofpmax)(str, endptr, 0 __LOCALE_ARG );
  444. y = (float) x;
  445. __fp_range_check(y, x);
  446. return y;
  447. #endif
  448. }
  449. libc_hidden_def(__XL_NPP(strtof))
  450. #endif
  451. #endif
  452. /**********************************************************************/
  453. #if defined(L_strtod) || defined(L_strtod_l) || defined(L_wcstod) || defined(L_wcstod_l)
  454. #if defined(NEED_STRTOD_WRAPPER)
  455. #if defined(L_wcstod) || defined(L_wcstod_l)
  456. #define strtod wcstod
  457. #define strtod_l wcstod_l
  458. #define __strtofpmax __wcstofpmax
  459. #define __strtofpmax_l __wcstofpmax_l
  460. #define Wchar wchar_t
  461. #else
  462. #define Wchar char
  463. #endif
  464. libc_hidden_proto(__XL_NPP(strtod))
  465. double __XL_NPP(strtod)(const Wchar *__restrict str,
  466. Wchar **__restrict endptr __LOCALE_PARAM )
  467. {
  468. #if FPMAX_TYPE == 2
  469. return __XL_NPP(__strtofpmax)(str, endptr, 0 __LOCALE_ARG );
  470. #else
  471. __fpmax_t x;
  472. double y;
  473. x = __XL_NPP(__strtofpmax)(str, endptr, 0 __LOCALE_ARG );
  474. y = (double) x;
  475. __fp_range_check(y, x);
  476. return y;
  477. #endif
  478. }
  479. libc_hidden_def(__XL_NPP(strtod))
  480. #endif
  481. #endif
  482. /**********************************************************************/
  483. #if defined(L_strtold) || defined(L_strtold_l) || defined(L_wcstold) || defined(L_wcstold_l)
  484. #if defined(NEED_STRTOLD_WRAPPER)
  485. #if defined(L_wcstold) || defined(L_wcstold_l)
  486. #define strtold wcstold
  487. #define strtold_l wcstold_l
  488. #define __strtofpmax __wcstofpmax
  489. #define __strtofpmax_l __wcstofpmax_l
  490. #define Wchar wchar_t
  491. #else
  492. #define Wchar char
  493. #endif
  494. libc_hidden_proto(__XL_NPP(strtold))
  495. long double __XL_NPP(strtold) (const Wchar *str, Wchar **endptr __LOCALE_PARAM )
  496. {
  497. #if FPMAX_TYPE == 3
  498. return __XL_NPP(__strtofpmax)(str, endptr, 0 __LOCALE_ARG );
  499. #else
  500. __fpmax_t x;
  501. long double y;
  502. x = __XL_NPP(__strtofpmax)(str, endptr, 0 __LOCALE_ARG );
  503. y = (long double) x;
  504. __fp_range_check(y, x);
  505. return y;
  506. #endif
  507. }
  508. libc_hidden_def(__XL_NPP(strtold))
  509. #endif
  510. #endif
  511. /**********************************************************************/