From owner-freebsd-numerics@FreeBSD.ORG Fri Aug 23 17:15:27 2013 Return-Path: Delivered-To: freebsd-numerics@FreeBSD.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) (using TLSv1 with cipher ADH-AES256-SHA (256/256 bits)) (No client certificate requested) by hub.freebsd.org (Postfix) with ESMTP id 13FDA758 for ; Fri, 23 Aug 2013 17:15:27 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (No client certificate requested) by mx1.freebsd.org (Postfix) with ESMTPS id CE955252C for ; Fri, 23 Aug 2013 17:15:26 +0000 (UTC) Received: from troutmask.apl.washington.edu (localhost.apl.washington.edu [127.0.0.1]) by troutmask.apl.washington.edu (8.14.7/8.14.7) with ESMTP id r7NHFQqn047750; Fri, 23 Aug 2013 10:15:26 -0700 (PDT) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.14.7/8.14.7/Submit) id r7NHFQrY047749; Fri, 23 Aug 2013 10:15:26 -0700 (PDT) (envelope-from sgk) Date: Fri, 23 Aug 2013 10:15:26 -0700 From: Steve Kargl To: Bruce Evans Subject: Re: (2nd time) tweaks to erff() threshold values Message-ID: <20130823171526.GA47736@troutmask.apl.washington.edu> References: <20130822213315.GA6708@troutmask.apl.washington.edu> <20130823202257.Q1593@besplex.bde.org> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: <20130823202257.Q1593@besplex.bde.org> User-Agent: Mutt/1.5.21 (2010-09-15) Cc: freebsd-numerics@FreeBSD.org X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.14 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Fri, 23 Aug 2013 17:15:27 -0000 On Fri, Aug 23, 2013 at 09:12:33PM +1000, Bruce Evans wrote: > > @ if(ix < 0x3f580000) { /* |x|<0.84375 */ > @ if(ix < 0x38800000) { /* |x|<2**-14 */ > @ if (ix < 0x04000000) /* |x|<0x1p-119 */ > @ return (8*x+efx8*x)/8; /*avoid spurious underflow */ To keep the diff between s_erf.c and s_erff.c small, do you want my to update s_erf.c with the last line above? -- Steve