(2nd time) tweaks to erff() threshold values
    Bruce Evans 
    brde at optusnet.com.au
       
    Mon Aug 26 14:30:25 UTC 2013
    
    
  
On Fri, 23 Aug 2013, Steve Kargl wrote:
> 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?
Yes, including the style change from *0.125 to /8 and any similar style
changes you can find (or don't change this anywhere), and the comment
on the threshold for underflow.
Bruce
    
    
More information about the freebsd-numerics
mailing list