From owner-cvs-src@FreeBSD.ORG  Tue Jun  1 12:03:31 2004
Return-Path: <owner-cvs-src@FreeBSD.ORG>
Delivered-To: cvs-src@freebsd.org
Received: from mx1.FreeBSD.org (mx1.freebsd.org [216.136.204.125])
	by hub.freebsd.org (Postfix) with ESMTP
	id E23BF16A4CE; Tue,  1 Jun 2004 12:03:31 -0700 (PDT)
Received: from repoman.freebsd.org (repoman.freebsd.org [216.136.204.115])
	by mx1.FreeBSD.org (Postfix) with ESMTP
	id DBB5243D5E; Tue,  1 Jun 2004 12:03:31 -0700 (PDT)
	(envelope-from bde@FreeBSD.org)
Received: from repoman.freebsd.org (localhost [127.0.0.1])
	by repoman.freebsd.org (8.12.11/8.12.11) with ESMTP id i51J3VZV066312;
	Tue, 1 Jun 2004 12:03:31 -0700 (PDT)
	(envelope-from bde@repoman.freebsd.org)
Received: (from bde@localhost)
	by repoman.freebsd.org (8.12.11/8.12.11/Submit) id i51J3VHL066311;
	Tue, 1 Jun 2004 12:03:31 -0700 (PDT)
	(envelope-from bde)
Message-Id: <200406011903.i51J3VHL066311@repoman.freebsd.org>
From: Bruce Evans <bde@FreeBSD.org>
Date: Tue, 1 Jun 2004 12:03:31 -0700 (PDT)
To: src-committers@FreeBSD.org, cvs-src@FreeBSD.org,
	cvs-all@FreeBSD.org
X-FreeBSD-CVS-Branch: HEAD
Subject: cvs commit: src/lib/msun/src e_powf.c
X-BeenThere: cvs-src@freebsd.org
X-Mailman-Version: 2.1.1
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: Tue, 01 Jun 2004 19:03:32 -0000

bde         2004/06/01 12:03:31 PDT

  FreeBSD src repository

  Modified files:
    lib/msun/src         e_powf.c 
  Log:
  Fixed another precision bug in powf().  This one is in the computation
  [t=p_l+p_h High].  We multiply t by lg2_h, and want the result to be
  exact.  For the bogus float case of the high-low decomposition trick,
  we normally discard the lowest 12 bits of the fraction for the high
  part, keeping 12 bits of precision.  That was used for t here, but it
  doesnt't work because for some reason we only discard the lowest 9
  bits in the fraction for lg2_h.  Discard another 3 bits of the fraction
  for t to compensate.
  
  This bug gave wrong results like:
  
        powf(0.9999999, -2.9999995) = 1.0000002 (should be 1.0000001)
          hex values: 3F7FFFFF C03FFFFE 3F800002 3F800001
  
  As explained in the log for the previous commit, the bug is normally
  masked by doing float calculations in extra precision on i386's, but
  is easily detected by ucbtest on systems that don't have accidental
  extra precision.
  
  This completes fixing all the bugs in powf() that were routinely found
  by ucbtest.
  
  Revision  Changes    Path
  1.11      +1 -1      src/lib/msun/src/e_powf.c