strtod.c 16 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619
  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. #endif
  96. #ifdef __UCLIBC_HAS_XLOCALE__
  97. #include <xlocale.h>
  98. #endif /* __UCLIBC_HAS_XLOCALE__ */
  99. /* Handle _STRTOD_HEXADECIMAL_FLOATS via uClibc config now. */
  100. #undef _STRTOD_HEXADECIMAL_FLOATS
  101. #ifdef __UCLIBC_HAS_HEXADECIMAL_FLOATS__
  102. #define _STRTOD_HEXADECIMAL_FLOATS 1
  103. #endif /* __UCLIBC_HAS_HEXADECIMAL_FLOATS__ */
  104. /**********************************************************************/
  105. #undef _STRTOD_FPMAX
  106. #if FPMAX_TYPE == 3
  107. #define NEED_STRTOLD_WRAPPER
  108. #define NEED_STRTOD_WRAPPER
  109. #define NEED_STRTOF_WRAPPER
  110. #elif FPMAX_TYPE == 2
  111. #define NEED_STRTOD_WRAPPER
  112. #define NEED_STRTOF_WRAPPER
  113. #elif FPMAX_TYPE == 1
  114. #define NEED_STRTOF_WRAPPER
  115. #else
  116. #error unknown FPMAX_TYPE!
  117. #endif
  118. extern void __fp_range_check(__fpmax_t y, __fpmax_t x) attribute_hidden;
  119. /**********************************************************************/
  120. #ifdef _STRTOD_RESTRICT_DIGITS
  121. #define EXP_DENORM_ADJUST DECIMAL_DIG
  122. #define MAX_ALLOWED_EXP (DECIMAL_DIG + EXP_DENORM_ADJUST - FPMAX_MIN_10_EXP)
  123. #if MAX_ALLOWED_EXP > INT_MAX
  124. #error size assumption violated for MAX_ALLOWED_EXP
  125. #endif
  126. #else
  127. /* We want some excess if we're not restricting mantissa digits. */
  128. #define MAX_ALLOWED_EXP ((20 - FPMAX_MIN_10_EXP) * 2)
  129. #endif
  130. #if defined(_STRTOD_RESTRICT_DIGITS) || defined(_STRTOD_ENDPTR) || defined(_STRTOD_HEXADECIMAL_FLOATS)
  131. #undef _STRTOD_NEED_NUM_DIGITS
  132. #define _STRTOD_NEED_NUM_DIGITS 1
  133. #endif
  134. /**********************************************************************/
  135. #if defined(L___strtofpmax) || defined(L___strtofpmax_l) || defined(L___wcstofpmax) || defined(L___wcstofpmax_l)
  136. #if defined(L___wcstofpmax) || defined(L___wcstofpmax_l)
  137. #define __strtofpmax __wcstofpmax
  138. #define __strtofpmax_l __wcstofpmax_l
  139. #define Wchar wchar_t
  140. #ifdef __UCLIBC_DO_XLOCALE
  141. #define ISSPACE(C) iswspace_l((C), locale_arg)
  142. #else
  143. #define ISSPACE(C) iswspace((C))
  144. #endif
  145. #else /* defined(L___wcstofpmax) || defined(L___wcstofpmax_l) */
  146. #define Wchar char
  147. #ifdef __UCLIBC_DO_XLOCALE
  148. #define ISSPACE(C) isspace_l((C), locale_arg)
  149. #else
  150. #define ISSPACE(C) isspace((C))
  151. #endif
  152. #endif /* defined(L___wcstofpmax) || defined(L___wcstofpmax_l) */
  153. #if defined(__UCLIBC_HAS_XLOCALE__) && !defined(__UCLIBC_DO_XLOCALE)
  154. __fpmax_t __strtofpmax(const Wchar *str, Wchar **endptr, int exponent_power)
  155. {
  156. return __strtofpmax_l(str, endptr, exponent_power, __UCLIBC_CURLOCALE);
  157. }
  158. #else /* defined(__UCLIBC_HAS_XLOCALE__) && !defined(__UCLIBC_DO_XLOCALE) */
  159. __fpmax_t __XL_NPP(__strtofpmax)(const Wchar *str, Wchar **endptr, int exponent_power
  160. __LOCALE_PARAM )
  161. {
  162. __fpmax_t number;
  163. __fpmax_t p_base = 10; /* Adjusted to 16 in the hex case. */
  164. Wchar *pos0;
  165. #ifdef _STRTOD_ENDPTR
  166. Wchar *pos1;
  167. #endif
  168. Wchar *pos = (Wchar *) str;
  169. int exponent_temp;
  170. int negative; /* A flag for the number, a multiplier for the exponent. */
  171. #ifdef _STRTOD_NEED_NUM_DIGITS
  172. int num_digits;
  173. #endif
  174. #ifdef __UCLIBC_HAS_LOCALE__
  175. #if defined(L___wcstofpmax) || defined(L___wcstofpmax_l)
  176. wchar_t decpt_wc = __LOCALE_PTR->decimal_point;
  177. #else
  178. const char *decpt = __LOCALE_PTR->decimal_point;
  179. int decpt_len = __LOCALE_PTR->decimal_point_len;
  180. #endif
  181. #endif
  182. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  183. Wchar expchar = 'e';
  184. Wchar *poshex = NULL;
  185. __uint16_t is_mask = _ISdigit;
  186. #define EXPCHAR expchar
  187. #define IS_X_DIGIT(C) __isctype((C), is_mask)
  188. #else /* _STRTOD_HEXADECIMAL_FLOATS */
  189. #define EXPCHAR 'e'
  190. #define IS_X_DIGIT(C) isdigit((C))
  191. #endif /* _STRTOD_HEXADECIMAL_FLOATS */
  192. while (ISSPACE(*pos)) { /* Skip leading whitespace. */
  193. ++pos;
  194. }
  195. negative = 0;
  196. switch(*pos) { /* Handle optional sign. */
  197. case '-': negative = 1; /* Fall through to increment position. */
  198. case '+': ++pos;
  199. }
  200. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  201. if ((*pos == '0') && (((pos[1])|0x20) == 'x')) {
  202. poshex = ++pos; /* Save position of 'x' in case no digits */
  203. ++pos; /* and advance past it. */
  204. is_mask = _ISxdigit; /* Used by IS_X_DIGIT. */
  205. expchar = 'p'; /* Adjust exponent char. */
  206. p_base = 16; /* Adjust base multiplier. */
  207. }
  208. #endif
  209. number = 0.;
  210. #ifdef _STRTOD_NEED_NUM_DIGITS
  211. num_digits = -1;
  212. #endif
  213. /* exponent_power = 0; */
  214. pos0 = NULL;
  215. LOOP:
  216. while (IS_X_DIGIT(*pos)) { /* Process string of (hex) digits. */
  217. #ifdef _STRTOD_RESTRICT_DIGITS
  218. if (num_digits < 0) { /* First time through? */
  219. ++num_digits; /* We've now seen a digit. */
  220. }
  221. if (num_digits || (*pos != '0')) { /* Had/have nonzero. */
  222. ++num_digits;
  223. if (num_digits <= DECIMAL_DIG) { /* Is digit significant? */
  224. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  225. number = number * p_base
  226. + (isdigit(*pos)
  227. ? (*pos - '0')
  228. : (((*pos)|0x20) - ('a' - 10)));
  229. #else /* _STRTOD_HEXADECIMAL_FLOATS */
  230. number = number * p_base + (*pos - '0');
  231. #endif /* _STRTOD_HEXADECIMAL_FLOATS */
  232. }
  233. }
  234. #else /* _STRTOD_RESTRICT_DIGITS */
  235. #ifdef _STRTOD_NEED_NUM_DIGITS
  236. ++num_digits;
  237. #endif
  238. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  239. number = number * p_base
  240. + (isdigit(*pos)
  241. ? (*pos - '0')
  242. : (((*pos)|0x20) - ('a' - 10)));
  243. #else /* _STRTOD_HEXADECIMAL_FLOATS */
  244. number = number * p_base + (*pos - '0');
  245. #endif /* _STRTOD_HEXADECIMAL_FLOATS */
  246. #endif /* _STRTOD_RESTRICT_DIGITS */
  247. ++pos;
  248. }
  249. #ifdef __UCLIBC_HAS_LOCALE__
  250. #if defined(L___wcstofpmax) || defined(L___wcstofpmax_l)
  251. if (!pos0 && (*pos == decpt_wc)) { /* First decimal point? */
  252. pos0 = ++pos;
  253. goto LOOP;
  254. }
  255. #else
  256. if (!pos0 && !memcmp(pos, decpt, decpt_len)) { /* First decimal point? */
  257. pos0 = (pos += decpt_len);
  258. goto LOOP;
  259. }
  260. #endif
  261. #else /* __UCLIBC_HAS_LOCALE__ */
  262. if ((*pos == '.') && !pos0) { /* First decimal point? */
  263. pos0 = ++pos; /* Save position of decimal point */
  264. goto LOOP; /* and process rest of digits. */
  265. }
  266. #endif /* __UCLIBC_HAS_LOCALE__ */
  267. #ifdef _STRTOD_NEED_NUM_DIGITS
  268. if (num_digits<0) { /* Must have at least one digit. */
  269. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  270. if (poshex) { /* Back up to '0' in '0x' prefix. */
  271. pos = poshex;
  272. goto DONE;
  273. }
  274. #endif /* _STRTOD_HEXADECIMAL_FLOATS */
  275. #ifdef _STRTOD_NAN_INF_STRINGS
  276. if (!pos0) { /* No decimal point, so check for inf/nan. */
  277. /* Note: nan is the first string so 'number = i/0.;' works. */
  278. static const char nan_inf_str[] = "\05nan\0\012infinity\0\05inf\0";
  279. int i = 0;
  280. #ifdef __UCLIBC_HAS_LOCALE__
  281. /* Avoid tolower problems for INFINITY in the tr_TR locale. (yuk)*/
  282. #undef _tolower
  283. #define _tolower(C) ((C)|0x20)
  284. #endif /* __UCLIBC_HAS_LOCALE__ */
  285. do {
  286. /* Unfortunately, we have no memcasecmp(). */
  287. int j = 0;
  288. while (_tolower(pos[j]) == nan_inf_str[i+1+j]) {
  289. ++j;
  290. if (!nan_inf_str[i+1+j]) {
  291. number = i / 0.;
  292. if (negative) { /* Correct for sign. */
  293. number = -number;
  294. }
  295. pos += nan_inf_str[i] - 2;
  296. goto DONE;
  297. }
  298. }
  299. i += nan_inf_str[i];
  300. } while (nan_inf_str[i]);
  301. }
  302. #endif /* STRTOD_NAN_INF_STRINGS */
  303. #ifdef _STRTOD_ENDPTR
  304. pos = (Wchar *) str;
  305. #endif
  306. goto DONE;
  307. }
  308. #endif /* _STRTOD_NEED_NUM_DIGITS */
  309. #ifdef _STRTOD_RESTRICT_DIGITS
  310. if (num_digits > DECIMAL_DIG) { /* Adjust exponent for skipped digits. */
  311. exponent_power += num_digits - DECIMAL_DIG;
  312. }
  313. #endif
  314. if (pos0) {
  315. exponent_power += pos0 - pos; /* Adjust exponent for decimal point. */
  316. }
  317. #ifdef _STRTOD_HEXADECIMAL_FLOATS
  318. if (poshex) {
  319. exponent_power *= 4; /* Above is 2**4, but below is 2. */
  320. p_base = 2;
  321. }
  322. #endif /* _STRTOD_HEXADECIMAL_FLOATS */
  323. if (negative) { /* Correct for sign. */
  324. number = -number;
  325. }
  326. /* process an exponent string */
  327. if (((*pos)|0x20) == EXPCHAR) {
  328. #ifdef _STRTOD_ENDPTR
  329. pos1 = pos;
  330. #endif
  331. negative = 1;
  332. switch(*++pos) { /* Handle optional sign. */
  333. case '-': negative = -1; /* Fall through to increment pos. */
  334. case '+': ++pos;
  335. }
  336. pos0 = pos;
  337. exponent_temp = 0;
  338. while (isdigit(*pos)) { /* Process string of digits. */
  339. #ifdef _STRTOD_RESTRICT_EXP
  340. if (exponent_temp < MAX_ALLOWED_EXP) { /* Avoid overflow. */
  341. exponent_temp = exponent_temp * 10 + (*pos - '0');
  342. }
  343. #else
  344. exponent_temp = exponent_temp * 10 + (*pos - '0');
  345. #endif
  346. ++pos;
  347. }
  348. #ifdef _STRTOD_ENDPTR
  349. if (pos == pos0) { /* No digits? */
  350. pos = pos1; /* Back up to {e|E}/{p|P}. */
  351. } /* else */
  352. #endif
  353. exponent_power += negative * exponent_temp;
  354. }
  355. #ifdef _STRTOD_ZERO_CHECK
  356. if (number == 0.) {
  357. goto DONE;
  358. }
  359. #endif
  360. /* scale the result */
  361. #ifdef _STRTOD_LOG_SCALING
  362. exponent_temp = exponent_power;
  363. if (exponent_temp < 0) {
  364. exponent_temp = -exponent_temp;
  365. }
  366. while (exponent_temp) {
  367. if (exponent_temp & 1) {
  368. if (exponent_power < 0) {
  369. /* Warning... caluclating a factor for the exponent and
  370. * then dividing could easily be faster. But doing so
  371. * might cause problems when dealing with denormals. */
  372. number /= p_base;
  373. } else {
  374. number *= p_base;
  375. }
  376. }
  377. exponent_temp >>= 1;
  378. p_base *= p_base;
  379. }
  380. #else /* _STRTOD_LOG_SCALING */
  381. while (exponent_power) {
  382. if (exponent_power < 0) {
  383. number /= p_base;
  384. exponent_power++;
  385. } else {
  386. number *= p_base;
  387. exponent_power--;
  388. }
  389. }
  390. #endif /* _STRTOD_LOG_SCALING */
  391. #ifdef _STRTOD_ERRNO
  392. if (__FPMAX_ZERO_OR_INF_CHECK(number)) {
  393. __set_errno(ERANGE);
  394. }
  395. #endif
  396. DONE:
  397. #ifdef _STRTOD_ENDPTR
  398. if (endptr) {
  399. *endptr = pos;
  400. }
  401. #endif
  402. return number;
  403. }
  404. #endif /* defined(__UCLIBC_HAS_XLOCALE__) && !defined(__UCLIBC_DO_XLOCALE) */
  405. #endif
  406. /**********************************************************************/
  407. #ifdef L___fp_range_check
  408. #if defined(NEED_STRTOF_WRAPPER) || defined(NEED_STRTOD_WRAPPER)
  409. void attribute_hidden __fp_range_check(__fpmax_t y, __fpmax_t x)
  410. {
  411. if (__FPMAX_ZERO_OR_INF_CHECK(y) /* y is 0 or +/- infinity */
  412. && (y != 0) /* y is not 0 (could have x>0, y==0 if underflow) */
  413. && !__FPMAX_ZERO_OR_INF_CHECK(x) /* x is not 0 or +/- infinity */
  414. ) {
  415. __set_errno(ERANGE); /* Then x is not in y's range. */
  416. }
  417. }
  418. #endif
  419. #endif
  420. /**********************************************************************/
  421. #if defined(L_strtof) || defined(L_strtof_l) || defined(L_wcstof) || defined(L_wcstof_l)
  422. #if defined(NEED_STRTOF_WRAPPER)
  423. #if defined(L_wcstof) || defined(L_wcstof_l)
  424. #define strtof wcstof
  425. #define strtof_l wcstof_l
  426. #define __strtof __wcstof
  427. #define __strtof_l __wcstof_l
  428. #define __strtofpmax __wcstofpmax
  429. #define __strtofpmax_l __wcstofpmax_l
  430. #define Wchar wchar_t
  431. #else
  432. #define Wchar char
  433. #endif
  434. float __XL(strtof)(const Wchar *str, Wchar **endptr __LOCALE_PARAM )
  435. {
  436. #if FPMAX_TYPE == 1
  437. return __XL_NPP(__strtofpmax)(str, endptr, 0 __LOCALE_ARG );
  438. #else
  439. __fpmax_t x;
  440. float y;
  441. x = __XL_NPP(__strtofpmax)(str, endptr, 0 __LOCALE_ARG );
  442. y = (float) x;
  443. __fp_range_check(y, x);
  444. return y;
  445. #endif
  446. }
  447. __XL_ALIAS(strtof)
  448. #endif
  449. #endif
  450. /**********************************************************************/
  451. #if defined(L_strtod) || defined(L_strtod_l) || defined(L_wcstod) || defined(L_wcstod_l)
  452. #if defined(NEED_STRTOD_WRAPPER)
  453. #if defined(L_wcstod) || defined(L_wcstod_l)
  454. #define strtod wcstod
  455. #define strtod_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. double __XL(strtod)(const Wchar *__restrict str,
  465. Wchar **__restrict endptr __LOCALE_PARAM )
  466. {
  467. #if FPMAX_TYPE == 2
  468. return __XL_NPP(__strtofpmax)(str, endptr, 0 __LOCALE_ARG );
  469. #else
  470. __fpmax_t x;
  471. double y;
  472. x = __XL_NPP(__strtofpmax)(str, endptr, 0 __LOCALE_ARG );
  473. y = (double) x;
  474. __fp_range_check(y, x);
  475. return y;
  476. #endif
  477. }
  478. __XL_ALIAS(strtod)
  479. #endif
  480. #endif
  481. /**********************************************************************/
  482. #if defined(L_strtold) || defined(L_strtold_l) || defined(L_wcstold) || defined(L_wcstold_l)
  483. #if defined(NEED_STRTOLD_WRAPPER)
  484. #if defined(L_wcstold) || defined(L_wcstold_l)
  485. #define strtold wcstold
  486. #define strtold_l wcstold_l
  487. #define __strtold __wcstold
  488. #define __strtold_l __wcstold_l
  489. #define __strtofpmax __wcstofpmax
  490. #define __strtofpmax_l __wcstofpmax_l
  491. #define Wchar wchar_t
  492. #else
  493. #define Wchar char
  494. #endif
  495. long double __XL(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. __XL_ALIAS(strtold)
  509. #endif
  510. #endif
  511. /**********************************************************************/