From owner-cvs-src@FreeBSD.ORG  Sat Feb  9 12:53:16 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 0DF6016A468;
	Sat,  9 Feb 2008 12:53:16 +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 00DF113C45E;
	Sat,  9 Feb 2008 12:53:16 +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 m19CrFGO022358;
	Sat, 9 Feb 2008 12:53:15 GMT (envelope-from bde@repoman.freebsd.org)
Received: (from bde@localhost)
	by repoman.freebsd.org (8.14.1/8.14.1/Submit) id m19CrFkb022357;
	Sat, 9 Feb 2008 12:53:15 GMT (envelope-from bde)
Message-Id: <200802091253.m19CrFkb022357@repoman.freebsd.org>
From: Bruce Evans <bde@FreeBSD.org>
Date: Sat, 9 Feb 2008 12:53:15 +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_expm1f.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: Sat, 09 Feb 2008 12:53:16 -0000

bde         2008-02-09 12:53:15 UTC

  FreeBSD src repository

  Modified files:
    lib/msun/src         s_expm1f.c 
  Log:
  As usual, use a minimax polynomial that is specialized for float
  precision.  The new polynomial has degree 4 instead of 10, and a maximum
  error of 2**-30.04 ulps instead of 2**-33.15.  This doesn't affect the
  final error significantly; the maximum error was and is about 0.5015
  ulps on i386 -O1, and the number of cases with an error of > 0.5 ulps
  is increased from 13851 to 14407.
  
  Note that the error is only this close to 0.5 ulps due to excessive
  extra precision caused by compiler bugs on i386.  The extra precision
  could be obtained intentionally, and is useful for keeping the error
  of the hyperbolic float functions below 1 ulp, since these functions
  are implemented using expm1f.  My recent change for scaling by 2**k
  had the unintentional side effect of retaining extra precision for
  longer, so callers of expm1f see errors of more like 0.0015 ulps than
  0.5015 ulps, and for the hyperbolic functions this reduces the maximum
  error from nearly about 2 ulps to about 0.75 ulps.
  
  This is about 10% faster on i386 (A64).  expm1* is still very slow,
  but now the float version is actually significantly faster.  The
  algorithm is very sophisticated but not very good except on machines
  with fast division.
  
  Revision  Changes    Path
  1.9       +8 -7      src/lib/msun/src/s_expm1f.c