From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: (qmail 7700 invoked by alias); 2 Dec 2019 09:05:12 -0000 Mailing-List: contact newlib-cvs-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: newlib-cvs-owner@sourceware.org Received: (qmail 7660 invoked by uid 9078); 2 Dec 2019 09:05:12 -0000 Date: Mon, 02 Dec 2019 09:05:00 -0000 Message-ID: <20191202090512.7658.qmail@sourceware.org> Content-Type: text/plain; charset="us-ascii" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit From: Corinna Vinschen To: newlib-cvs@sourceware.org Subject: [newlib-cygwin] libm: switch sf_log1p from double error routines to float X-Act-Checkin: newlib-cygwin X-Git-Author: Keith Packard X-Git-Refname: refs/heads/master X-Git-Oldrev: 8574f8a1e4b02a4b47e2434111c4e2c382a9f046 X-Git-Newrev: 7a526cdc28a3c4acce98e8a99b06562452c90d07 X-SW-Source: 2019-q4/txt/msg00018.txt.bz2 https://sourceware.org/git/gitweb.cgi?p=newlib-cygwin.git;h=7a526cdc28a3c4acce98e8a99b06562452c90d07 commit 7a526cdc28a3c4acce98e8a99b06562452c90d07 Author: Keith Packard Date: Fri Nov 29 11:23:20 2019 -0800 libm: switch sf_log1p from double error routines to float sf_log1p was using __math_divzero and __math_invalid, which drag in a pile of double-precision code. Switch to using the single-precision variants. This also required making those available in __OBSOLETE_MATH mode. Signed-off-by: Keith Packard Diff: --- newlib/libm/common/math_errf.c | 2 -- newlib/libm/common/sf_log1p.c | 4 ++-- 2 files changed, 2 insertions(+), 4 deletions(-) diff --git a/newlib/libm/common/math_errf.c b/newlib/libm/common/math_errf.c index 762fc27..53c68b1 100644 --- a/newlib/libm/common/math_errf.c +++ b/newlib/libm/common/math_errf.c @@ -75,7 +75,6 @@ __math_oflowf (uint32_t sign) return xflowf (sign, 0x1p97f); } -#if !__OBSOLETE_MATH HIDDEN float __math_divzerof (uint32_t sign) { @@ -89,4 +88,3 @@ __math_invalidf (float x) float y = (x - x) / (x - x); return isnan (x) ? y : with_errnof (y, EDOM); } -#endif /* !__OBSOLETE_MATH */ diff --git a/newlib/libm/common/sf_log1p.c b/newlib/libm/common/sf_log1p.c index d867680..d326e00 100644 --- a/newlib/libm/common/sf_log1p.c +++ b/newlib/libm/common/sf_log1p.c @@ -56,9 +56,9 @@ static float zero = 0.0; if (hx < 0x3ed413d7) { /* x < 0.41422 */ if(ax>=0x3f800000) { /* x <= -1.0 */ if(x==(float)-1.0) - return __math_divzero (1); /* log1p(-1)=-inf */ + return __math_divzerof (1); /* log1p(-1)=-inf */ else - return __math_invalid (x); /* log1p(x<-1)=NaN */ + return __math_invalidf (x); /* log1p(x<-1)=NaN */ } if(ax<0x31000000) { /* |x| < 2**-29 */ if(two25+x>zero /* raise inexact */