From owner-cvs-src@FreeBSD.ORG  Mon Feb 11 05:20:02 2008
Return-Path: <owner-cvs-src@FreeBSD.ORG>
Delivered-To: cvs-src@FreeBSD.org
Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:4f8:fff6::34])
	by hub.freebsd.org (Postfix) with ESMTP id 8B78C16A417;
	Mon, 11 Feb 2008 05:20:02 +0000 (UTC) (envelope-from bde@FreeBSD.org)
Received: from repoman.freebsd.org (repoman.freebsd.org
	[IPv6:2001:4f8:fff6::29])
	by mx1.freebsd.org (Postfix) with ESMTP id 753D513C448;
	Mon, 11 Feb 2008 05:20:02 +0000 (UTC) (envelope-from bde@FreeBSD.org)
Received: from repoman.freebsd.org (localhost [127.0.0.1])
	by repoman.freebsd.org (8.14.1/8.14.1) with ESMTP id m1B5K2GS041237;
	Mon, 11 Feb 2008 05:20:02 GMT (envelope-from bde@repoman.freebsd.org)
Received: (from bde@localhost)
	by repoman.freebsd.org (8.14.1/8.14.1/Submit) id m1B5K2Xm041236;
	Mon, 11 Feb 2008 05:20:02 GMT (envelope-from bde)
Message-Id: <200802110520.m1B5K2Xm041236@repoman.freebsd.org>
From: Bruce Evans <bde@FreeBSD.org>
Date: Mon, 11 Feb 2008 05:20:02 +0000 (UTC)
To: src-committers@FreeBSD.org, cvs-src@FreeBSD.org, cvs-all@FreeBSD.org
X-FreeBSD-CVS-Branch: HEAD
Cc: 
Subject: cvs commit: src/lib/msun/src s_exp2f.c
X-BeenThere: cvs-src@freebsd.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: CVS commit messages for the src tree <cvs-src.freebsd.org>
List-Unsubscribe: <http://lists.freebsd.org/mailman/listinfo/cvs-src>,
	<mailto:cvs-src-request@freebsd.org?subject=unsubscribe>
List-Archive: <http://lists.freebsd.org/pipermail/cvs-src>
List-Post: <mailto:cvs-src@freebsd.org>
List-Help: <mailto:cvs-src-request@freebsd.org?subject=help>
List-Subscribe: <http://lists.freebsd.org/mailman/listinfo/cvs-src>,
	<mailto:cvs-src-request@freebsd.org?subject=subscribe>
X-List-Received-Date: Mon, 11 Feb 2008 05:20:02 -0000

bde         2008-02-11 05:20:02 UTC

  FreeBSD src repository

  Modified files:
    lib/msun/src         s_exp2f.c 
  Log:
  Use double precision for z and thus for the entire calculation of
  exp2(i/TBLSIZE) * p(z) instead of only for the final multiplication
  and addition.  This fixes the code to match the comment that the maximum
  error is 0.5010 ulps (except on machines that evaluate float expressions
  in extra precision, e.g., i386's, where the evaluation was already
  in extra precision).
  
  Fix and expand the comment about use of double precision.
  
  The relative roundoff error from evaluating p(z) in non-extra precision
  was about 16 times larger than in exp2() because the interval length
  is 16 times smaller.  Its maximum was at least P1 * (1.0 ulps) *
  max(|z|) ~= log(2) * 1.0 * 1/32 ~= 0.0217 ulps (1.0 ulps from the
  addition in (1 + P1*z) with a cancelation error when z ~= -1/32).  The
  actual final maximum was 0.5313 ulps, of which 0.0303 ulps must have
  come from the additional roundoff error in p(z).  I can't explain why
  the additional roundoff error was almost 3/2 times larger than the rough
  estimate.
  
  Revision  Changes    Path
  1.6       +4 -3      src/lib/msun/src/s_exp2f.c