From owner-freebsd-numerics@FreeBSD.ORG  Mon Aug 26 14:30:25 2013
Return-Path: <owner-freebsd-numerics@FreeBSD.ORG>
Delivered-To: freebsd-numerics@freebsd.org
Received: from mx1.freebsd.org (mx1.freebsd.org [8.8.178.115])
 (using TLSv1 with cipher ADH-AES256-SHA (256/256 bits))
 (No client certificate requested)
 by hub.freebsd.org (Postfix) with ESMTP id 6A9AAC1B
 for <freebsd-numerics@freebsd.org>; Mon, 26 Aug 2013 14:30:25 +0000 (UTC)
 (envelope-from brde@optusnet.com.au)
Received: from mail110.syd.optusnet.com.au (mail110.syd.optusnet.com.au
 [211.29.132.97]) by mx1.freebsd.org (Postfix) with ESMTP id 2F95A282E
 for <freebsd-numerics@freebsd.org>; Mon, 26 Aug 2013 14:30:24 +0000 (UTC)
Received: from c122-106-156-23.carlnfd1.nsw.optusnet.com.au
 (c122-106-156-23.carlnfd1.nsw.optusnet.com.au [122.106.156.23])
 by mail110.syd.optusnet.com.au (Postfix) with ESMTPS id E78DE783074;
 Tue, 27 Aug 2013 00:30:13 +1000 (EST)
Date: Tue, 27 Aug 2013 00:30:12 +1000 (EST)
From: Bruce Evans <brde@optusnet.com.au>
X-X-Sender: bde@besplex.bde.org
To: Steve Kargl <sgk@troutmask.apl.washington.edu>
Subject: Re: (2nd time) tweaks to erff() threshold values
In-Reply-To: <20130823171526.GA47736@troutmask.apl.washington.edu>
Message-ID: <20130827002821.K2328@besplex.bde.org>
References: <20130822213315.GA6708@troutmask.apl.washington.edu>
 <20130823202257.Q1593@besplex.bde.org>
 <20130823171526.GA47736@troutmask.apl.washington.edu>
MIME-Version: 1.0
Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed
X-Optus-CM-Score: 0
X-Optus-CM-Analysis: v=2.1 cv=bpB1Wiqi c=1 sm=1 tr=0
 a=ebeQFi2P/qHVC0Yw9JDJ4g==:117 a=PO7r1zJSAAAA:8 a=F6uwjYSw4mAA:10
 a=kj9zAlcOel0A:10 a=JzwRw_2MAAAA:8 a=bRjWgPqadKYA:10
 a=0N4hs-jxxM66dZRz4AkA:9 a=CjuIK1q_8ugA:10
Cc: freebsd-numerics@freebsd.org, Bruce Evans <brde@optusnet.com.au>
X-BeenThere: freebsd-numerics@freebsd.org
X-Mailman-Version: 2.1.14
Precedence: list
List-Id: "Discussions of high quality implementation of libm functions."
 <freebsd-numerics.freebsd.org>
List-Unsubscribe: <http://lists.freebsd.org/mailman/options/freebsd-numerics>, 
 <mailto:freebsd-numerics-request@freebsd.org?subject=unsubscribe>
List-Archive: <http://lists.freebsd.org/pipermail/freebsd-numerics>
List-Post: <mailto:freebsd-numerics@freebsd.org>
List-Help: <mailto:freebsd-numerics-request@freebsd.org?subject=help>
List-Subscribe: <http://lists.freebsd.org/mailman/listinfo/freebsd-numerics>, 
 <mailto:freebsd-numerics-request@freebsd.org?subject=subscribe>
X-List-Received-Date: Mon, 26 Aug 2013 14:30:25 -0000

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