From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: from conuserg-08.nifty.com (conuserg-08.nifty.com [210.131.2.75]) by sourceware.org (Postfix) with ESMTPS id CE9EA3858D3C for ; Thu, 2 Dec 2021 10:22:16 +0000 (GMT) DMARC-Filter: OpenDMARC Filter v1.4.1 sourceware.org CE9EA3858D3C Authentication-Results: sourceware.org; dmarc=fail (p=none dis=none) header.from=nifty.ne.jp Authentication-Results: sourceware.org; spf=fail smtp.mailfrom=nifty.ne.jp Received: from localhost.localdomain (z221123.dynamic.ppp.asahi-net.or.jp [110.4.221.123]) (authenticated) by conuserg-08.nifty.com with ESMTP id 1B2ALttX007541; Thu, 2 Dec 2021 19:22:01 +0900 DKIM-Filter: OpenDKIM Filter v2.10.3 conuserg-08.nifty.com 1B2ALttX007541 DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=nifty.ne.jp; s=dec2015msa; t=1638440522; bh=wbjclgG0i1CtGe7lX4MlURaPoQqNnJb1vr6HBZCHBHs=; h=From:To:Cc:Subject:Date:From; b=rb2+Eg5zON2WmsxArQOfF9O0Mk3NQw9g5ebacHlT16UavlCgmwW06wwQHNJOrn4oT PVWl5LIQiYy30WT1NqkQ7L7K0ZfsoiCDf4BLc06zJO+sVkPaxKOR4Ub+9H0dDhGMv8 v3yRj3maqXnMYuVicAl1mDrKIFDmfHq/U1IojbMYn4npscl1BtfkxBOmFSAw4/GKKY usFjRMs8VtjAcJ78sc/cN7wIFe2aRtfpPjWKP4RlcWKhGUZKm3WsGChdh6hn0AxQSO JRSpiub75Uq7lG5rVv5yaqMs68z5uoVcgAxwEw143BvS9N3mEwPSUj3f5vWbMeejbZ 9FTZPOJaYInBA== X-Nifty-SrcIP: [110.4.221.123] From: Takashi Yano To: newlib@sourceware.org Subject: [PATCH v2] frexpl: Support smaller long double of LDBL_MANT_DIG == 53. Date: Thu, 2 Dec 2021 19:21:46 +0900 Message-Id: <20211202102146.13630-1-takashi.yano@nifty.ne.jp> X-Mailer: git-send-email 2.33.0 MIME-Version: 1.0 Content-Transfer-Encoding: 8bit X-Spam-Status: No, score=-10.4 required=5.0 tests=BAYES_00, DKIM_SIGNED, DKIM_VALID, DKIM_VALID_AU, DKIM_VALID_EF, GIT_PATCH_0, RCVD_IN_DNSWL_NONE, SPF_HELO_NONE, SPF_PASS, TXREP autolearn=ham autolearn_force=no version=3.4.4 X-Spam-Checker-Version: SpamAssassin 3.4.4 (2020-01-24) on server2.sourceware.org X-BeenThere: newlib@sourceware.org X-Mailman-Version: 2.1.29 Precedence: list List-Id: Newlib mailing list List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Thu, 02 Dec 2021 10:22:20 -0000 - 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. --- newlib/libm/common/frexpl.c | 54 ++++++++++++++++++++++++++++--------- 1 file changed, 42 insertions(+), 12 deletions(-) diff --git a/newlib/libm/common/frexpl.c b/newlib/libm/common/frexpl.c index 22c811689..8da2823ca 100644 --- a/newlib/libm/common/frexpl.c +++ b/newlib/libm/common/frexpl.c @@ -33,16 +33,42 @@ POSSIBILITY OF SUCH DAMAGE. #include "local.h" /* 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 frexpl (long double x, int *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 { long double x; struct { @@ -68,7 +94,9 @@ union ldbl { # endif } 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 { long double x; struct { @@ -96,9 +124,11 @@ union ldbl { # endif } u32; }; +# else +# error Unsupported long double format. # endif -static const double two114 = 0x1p114; +static const int scale_exp = LDBL_MANT_DIG + 1; long double frexpl (long double x, int *eptr) @@ -107,16 +137,16 @@ frexpl (long double x, int *eptr) u.x = x; int e = u.u32.exp; *eptr = 0; - if (e == 0x7fff || x == 0) + if (e == (LDBL_MAX_EXP*2 - 1) || x == 0) return x; /* inf,nan,0 */ if (e == 0) /* subnormal */ { - u.x *= two114; + u.x *= scale; e = u.u32.exp; - *eptr -= 114; + *eptr -= scale_exp; } - *eptr += e - 16382; - u.u32.exp = 0x3ffe; /* 0 */ + *eptr += e - (LDBL_MAX_EXP - 2); + u.u32.exp = LDBL_MAX_EXP - 2; /* -1 */ return u.x; } -#endif /* End of 80-bit or 128-bit long double */ +#endif /* !_LDBL_EQ_DBL */ -- 2.33.0