aboutsummaryrefslogtreecommitdiff
path: root/newlib/libc
diff options
context:
space:
mode:
Diffstat (limited to 'newlib/libc')
-rw-r--r--newlib/libc/stdlib/Makefile.inc1
-rw-r--r--newlib/libc/stdlib/mprec.h2
-rw-r--r--newlib/libc/stdlib/strtold.c56
-rw-r--r--newlib/libc/stdlib/strtorQ.c129
4 files changed, 159 insertions, 29 deletions
diff --git a/newlib/libc/stdlib/Makefile.inc b/newlib/libc/stdlib/Makefile.inc
index 9812add..6510cdf 100644
--- a/newlib/libc/stdlib/Makefile.inc
+++ b/newlib/libc/stdlib/Makefile.inc
@@ -83,6 +83,7 @@ if HAVE_LONG_DOUBLE
libc_a_SOURCES += \
%D%/strtodg.c \
%D%/strtold.c \
+ %D%/strtorQ.c \
%D%/strtorx.c \
%D%/wcstold.c
endif # HAVE_LONG_DOUBLE
diff --git a/newlib/libc/stdlib/mprec.h b/newlib/libc/stdlib/mprec.h
index 83de932..c2d0a0a 100644
--- a/newlib/libc/stdlib/mprec.h
+++ b/newlib/libc/stdlib/mprec.h
@@ -385,6 +385,8 @@ double _strtod_l (struct _reent *ptr, const char *__restrict s00,
#if defined (_HAVE_LONG_DOUBLE) && !defined (_LDBL_EQ_DBL)
int _strtorx_l (struct _reent *, const char *, char **, int,
void *, locale_t);
+int _strtorQ_l (struct _reent *, const char *, char **, int,
+ void *, locale_t);
int _strtodg_l (struct _reent *p, const char *s00, char **se,
struct FPI *fpi, Long *exp, __ULong *bits,
locale_t);
diff --git a/newlib/libc/stdlib/strtold.c b/newlib/libc/stdlib/strtold.c
index 6bd1c2c..57cc38e 100644
--- a/newlib/libc/stdlib/strtold.c
+++ b/newlib/libc/stdlib/strtold.c
@@ -58,47 +58,45 @@ __flt_rounds(void)
#define FLT_ROUNDS 0
#endif
-long double
-_strtold_r (struct _reent *ptr, const char *__restrict s00,
- char **__restrict se)
-{
-#ifdef _LDBL_EQ_DBL
- /* On platforms where long double is as wide as double. */
- return _strtod_l (ptr, s00, se, __get_current_locale ());
-#else
+/*
+ * The core implementation for the strtold family of functions, handling
+ * different long double formats.
+ */
+static long double
+_strtold_impl(struct _reent *ptr, const char *__restrict s00,
+ char **__restrict se, locale_t loc) {
+#if defined(_LDBL_EQ_DBL)
+ /* On platforms where long double is as wide as double. */
+ return _strtod_l(ptr, s00, se, loc);
+#elif LDBL_MANT_DIG == 64 /* For 80 bit long doubles */
long double result;
-
- _strtorx_l (ptr, s00, se, FLT_ROUNDS, &result, __get_current_locale ());
+ _strtorx_l(ptr, s00, se, FLT_ROUNDS, &result, loc);
return result;
+#elif LDBL_MANT_DIG == 113 /* For 128bit long doubles */
+ long double result;
+ _strtorQ_l(ptr, s00, se, FLT_ROUNDS, &result, loc);
+ return result;
+#else
+#warning "strtold not implemented for this long double format"
+ return (long double)_strtod_l(ptr, s00, se, loc);
#endif
}
long double
-strtold_l (const char *__restrict s00, char **__restrict se, locale_t loc)
-{
-#ifdef _LDBL_EQ_DBL
- /* On platforms where long double is as wide as double. */
- return _strtod_l (_REENT, s00, se, loc);
-#else
- long double result;
+_strtold_r(struct _reent *ptr, const char *__restrict s00,
+ char **__restrict se) {
+ return _strtold_impl(ptr, s00, se, __get_current_locale());
+}
- _strtorx_l (_REENT, s00, se, FLT_ROUNDS, &result, loc);
- return result;
-#endif
+long double
+strtold_l(const char *__restrict s00, char **__restrict se, locale_t loc) {
+ return _strtold_impl(_REENT, s00, se, loc);
}
long double
strtold (const char *__restrict s00, char **__restrict se)
{
-#ifdef _LDBL_EQ_DBL
- /* On platforms where long double is as wide as double. */
- return _strtod_l (_REENT, s00, se, __get_current_locale ());
-#else
- long double result;
-
- _strtorx_l (_REENT, s00, se, FLT_ROUNDS, &result, __get_current_locale ());
- return result;
-#endif
+ return _strtold_impl(_REENT, s00, se, __get_current_locale());
}
#endif /* _HAVE_LONG_DOUBLE */
diff --git a/newlib/libc/stdlib/strtorQ.c b/newlib/libc/stdlib/strtorQ.c
new file mode 100644
index 0000000..e9e0b9f
--- /dev/null
+++ b/newlib/libc/stdlib/strtorQ.c
@@ -0,0 +1,129 @@
+/****************************************************************
+
+The author of this software is David M. Gay.
+
+Copyright (C) 1998, 2000 by Lucent Technologies
+All Rights Reserved
+
+Permission to use, copy, modify, and distribute this software and
+its documentation for any purpose and without fee is hereby
+granted, provided that the above copyright notice appear in all
+copies and that both that the copyright notice and this
+permission notice and warranty disclaimer appear in supporting
+documentation, and that the name of Lucent or any of its entities
+not be used in advertising or publicity pertaining to
+distribution of the software without specific, written prior
+permission.
+
+LUCENT DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE,
+INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS.
+IN NO EVENT SHALL LUCENT OR ANY OF ITS ENTITIES BE LIABLE FOR ANY
+SPECIAL, INDIRECT OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
+WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER
+IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION,
+ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF
+THIS SOFTWARE.
+
+****************************************************************/
+
+/* Please send bug reports to David M. Gay (dmg at acm dot org,
+ * with " at " changed at "@" and " dot " changed to "."). */
+
+#include <_ansi.h>
+#include <errno.h>
+#include <stdlib.h>
+#include <string.h>
+#include "mprec.h"
+#include "gdtoa.h"
+
+#if defined (_HAVE_LONG_DOUBLE) && !defined (_LDBL_EQ_DBL)
+
+/* one or the other of IEEE_MC68k or IEEE_8087 should be #defined */
+
+#ifdef IEEE_MC68k
+#define _0 0
+#define _1 1
+#define _2 2
+#define _3 3
+#endif
+#ifdef IEEE_8087
+#define _0 3
+#define _1 2
+#define _2 1
+#define _3 0
+#endif
+
+ void
+#ifdef KR_headers
+ULtoQ(L, bits, exp, k) ULong *L; ULong *bits; Long exp; int k;
+#else
+ULtoQ(ULong *L, ULong *bits, Long exp, int k)
+#endif
+{
+ switch(k & STRTOG_Retmask) {
+ case STRTOG_NoNumber:
+ case STRTOG_Zero:
+ L[0] = L[1] = L[2] = L[3] = 0;
+ break;
+
+ case STRTOG_Normal:
+ L[_3] = bits[0];
+ L[_2] = bits[1];
+ L[_1] = bits[2];
+ L[_0] = (bits[3] & ~0x10000) | ((exp + 0x3fff + 112) << 16);
+ break;
+
+ case STRTOG_NaNbits:
+ L[_3] = bits[0];
+ L[_2] = bits[1];
+ L[_1] = bits[2];
+ L[_0] = (bits[3] & ~0x10000)
+ | (((exp + 0x3fff + 112) << 16) | (1 << 15));
+ break;
+
+ case STRTOG_Denormal:
+ L[_3] = bits[0];
+ L[_2] = bits[1];
+ L[_1] = bits[2];
+ L[_0] = bits[3];
+ break;
+
+ case STRTOG_Infinite:
+ L[_0] = 0x7fff0000;
+ L[_1] = L[_2] = L[_3] = 0;
+ break;
+
+ case STRTOG_NaN:
+ *((long double*)L) = __builtin_nanl ("");
+ }
+ if (k & STRTOG_Neg)
+ L[_0] |= 0x80000000L;
+ }
+
+ int
+#ifdef KR_headers
+_strtorQ_l(p, s, sp, rounding, L, locale) struct _reent *p; CONST char *s;
+char **sp; int rounding; void *L; locale_t locale;
+#else
+_strtorQ_l(struct _reent *p, CONST char *s, char **sp,
+ int rounding, void *L, locale_t locale)
+#endif
+{
+ static FPI fpi0 = { 113, 1-16383-113+1, 32766-16383-113+1, 1, SI };
+ FPI *fpi, fpi1;
+ ULong bits[4];
+ Long exp;
+ int k;
+
+ fpi = &fpi0;
+ if (rounding != FPI_Round_near) {
+ fpi1 = fpi0;
+ fpi1.rounding = rounding;
+ fpi = &fpi1;
+ }
+ k = _strtodg_l(p, s, sp, fpi, &exp, bits, locale);
+ ULtoQ((ULong*)L, bits, exp, k);
+ return k;
+ }
+
+#endif /* _HAVE_LONG_DOUBLE && !_LDBL_EQ_DBL */
generated by cgit v1.1 at 2025-09-27 11:30:39 +0000