Skip site navigation (1)Skip section navigation (2)
Date:      Fri, 23 Aug 2013 10:15:26 -0700
From:      Steve Kargl <sgk@troutmask.apl.washington.edu>
To:        Bruce Evans <brde@optusnet.com.au>
Cc:        freebsd-numerics@FreeBSD.org
Subject:   Re: (2nd time) tweaks to erff() threshold values
Message-ID:  <20130823171526.GA47736@troutmask.apl.washington.edu>
In-Reply-To: <20130823202257.Q1593@besplex.bde.org>
References:  <20130822213315.GA6708@troutmask.apl.washington.edu> <20130823202257.Q1593@besplex.bde.org>

Next in thread | Previous in thread | Raw E-Mail | Index | Archive | Help
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



Want to link to this message? Use this URL: <http://docs.FreeBSD.org/cgi/mid.cgi?20130823171526.GA47736>