libs/newlib-cygwin
libs/newlib-cygwin
4
0
Fork 0
mirror of git://sourceware.org/git/newlib-cygwin.git synced 2025-02-28 12:05:47 +08:00
Code Issues Projects Releases Wiki Activity

frexpl: Support smaller long double of LDBL_MANT_DIG == 53.

Browse Source 
- Currently, frexpl() supports only the following cases.
    1) LDBL_MANT_DIG == 64 or 113
    2) 'long double' is equivalent to 'double'
  This patch add support for LDBL_MANT_DIG == 53.
This commit is contained in:
Takashi Yano 2021-12-02 18:51:21 +09:00
parent dbc7d28ba9
commit 34876c05a5
1 changed files with 42 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

54
newlib/libm/common/frexpl.c
View File

@ -33,16 +33,42 @@ POSSIBILITY OF SUCH DAMAGE.
#include "local.h" #include "local.h"
/* On platforms where long double is as wide as double. */ /* On platforms where long double is as wide as double. */
#if defined(_LDBL_EQ_DBL) || (LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024) #if defined(_LDBL_EQ_DBL)
long double long double
frexpl (long double x, int *eptr) frexpl (long double x, int *eptr)
{ {
return frexp(x, eptr); return frexp(x, eptr);
} }
#endif #else /* !_DBL_EQ_DBL */
# if (LDBL_MANT_DIG == 53) /* 64-bit long double */
static const double scale = 0x1p54;
union ldbl {
long double x;
struct {
# ifdef __IEEE_LITTLE_ENDIAN /* for Intel CPU */
__uint32_t fracl;
__uint32_t frach:20;
__uint32_t exp:11;
__uint32_t sign:1;
# endif
# ifdef __IEEE_BIG_ENDIAN
__uint32_t sign:1;
__uint32_t exp:11;
__uint32_t frach:20;
# ifndef ___IEEE_BYTES_LITTLE_ENDIAN
# else /* ARMEL without __VFP_FP__ */
__uint32_t frach:20;
__uint32_t exp:11;
__uint32_t sign:1;
# endif
__uint32_t fracl;
# endif
} u32;
};
# elif (LDBL_MANT_DIG == 64) /* 80-bit long double */
static const double scale = 0x1p65;
#if (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
# if (LDBL_MANT_DIG == 64) /* 80-bit long double */
union ldbl { union ldbl {
long double x; long double x;
struct { struct {
@ -68,7 +94,9 @@ union ldbl {
# endif # endif
} u32; } u32;
}; };
# else /* LDBL_MANT_DIG == 113, 128-bit long double */ # elif (LDBL_MANT_DIG == 113) /* 128-bit long double */
static const double scale = 0x1p114;
union ldbl { union ldbl {
long double x; long double x;
struct { struct {
@ -96,9 +124,11 @@ union ldbl {
# endif # endif
} u32; } u32;
}; };
# else
# error Unsupported long double format.
# endif # endif
static const double two114 = 0x1p114; static const int scale_exp = LDBL_MANT_DIG + 1;
long double long double
frexpl (long double x, int *eptr) frexpl (long double x, int *eptr)
@ -107,16 +137,16 @@ frexpl (long double x, int *eptr)
u.x = x; u.x = x;
int e = u.u32.exp; int e = u.u32.exp;
*eptr = 0; *eptr = 0;
if (e == 0x7fff || x == 0) if (e == (LDBL_MAX_EXP*2 - 1) || x == 0)
return x; /* inf,nan,0 */ return x; /* inf,nan,0 */
if (e == 0) /* subnormal */ if (e == 0) /* subnormal */
{ {
u.x *= two114; u.x *= scale;
e = u.u32.exp; e = u.u32.exp;
*eptr -= 114; *eptr -= scale_exp;
} }
*eptr += e - 16382; *eptr += e - (LDBL_MAX_EXP - 2);
u.u32.exp = 0x3ffe; /* 0 */ u.u32.exp = LDBL_MAX_EXP - 2; /* -1 */
return u.x; return u.x;
} }
#endif /* End of 80-bit or 128-bit long double */ #endif /* !_LDBL_EQ_DBL */