(2nd time) tweaks to erff() threshold values

Steve Kargl sgk at troutmask.apl.washington.edu
Fri Aug 23 17:15:27 UTC 2013


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


More information about the freebsd-numerics mailing list