From owner-cvs-all@FreeBSD.ORG  Sun Nov 13 00:08:23 2005
Return-Path: <owner-cvs-all@FreeBSD.ORG>
X-Original-To: cvs-all@FreeBSD.org
Delivered-To: cvs-all@FreeBSD.org
Received: from mx1.FreeBSD.org (mx1.freebsd.org [216.136.204.125])
	by hub.freebsd.org (Postfix) with ESMTP id C2C6716A41F;
	Sun, 13 Nov 2005 00:08:23 +0000 (GMT) (envelope-from bde@FreeBSD.org)
Received: from repoman.freebsd.org (repoman.freebsd.org [216.136.204.115])
	by mx1.FreeBSD.org (Postfix) with ESMTP id 81BDB43D4C;
	Sun, 13 Nov 2005 00:08:23 +0000 (GMT) (envelope-from bde@FreeBSD.org)
Received: from repoman.freebsd.org (localhost [127.0.0.1])
	by repoman.freebsd.org (8.13.1/8.13.1) with ESMTP id jAD08Nlj059875;
	Sun, 13 Nov 2005 00:08:23 GMT (envelope-from bde@repoman.freebsd.org)
Received: (from bde@localhost)
	by repoman.freebsd.org (8.13.1/8.13.1/Submit) id jAD08NTn059874;
	Sun, 13 Nov 2005 00:08:23 GMT (envelope-from bde)
Message-Id: <200511130008.jAD08NTn059874@repoman.freebsd.org>
From: Bruce Evans <bde@FreeBSD.org>
Date: Sun, 13 Nov 2005 00:08:23 +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 e_coshf.c
X-BeenThere: cvs-all@freebsd.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: CVS commit messages for the entire tree <cvs-all.freebsd.org>
List-Unsubscribe: <http://lists.freebsd.org/mailman/listinfo/cvs-all>,
	<mailto:cvs-all-request@freebsd.org?subject=unsubscribe>
List-Archive: <http://lists.freebsd.org/pipermail/cvs-all>
List-Post: <mailto:cvs-all@freebsd.org>
List-Help: <mailto:cvs-all-request@freebsd.org?subject=help>
List-Subscribe: <http://lists.freebsd.org/mailman/listinfo/cvs-all>,
	<mailto:cvs-all-request@freebsd.org?subject=subscribe>
X-List-Received-Date: Sun, 13 Nov 2005 00:08:23 -0000

bde         2005-11-13 00:08:23 UTC

  FreeBSD src repository

  Modified files:
    lib/msun/src         e_coshf.c 
  Log:
  Fixed some magic numbers.
  
  The threshold for not being tiny was confusing and too small.  Use the
  usual 2**-12 threshold and simplify the algorithm slightly so that
  this threshold works (now use the threshold for sinhf() instead of one
  for 1+expm1()).  This is just a small optimization.
  
  The magic number 22 is log(DBL_EPSILON)/2 plus slop.  This is bogus
  for float precision.  Use 9 (~log(FLT_EPSILON)/2 plus less slop than
  for double precision).
  
  The threshold for switching from returning exp(x)/2 to returning
  exp(x/2)^2/2 was a little smaller than necessary.  This was not quite
  harmless since the exp(x/2)^2/2 case is inaccurate.  Fixing it happens
  to avoid accuracy problems for 2*6 of the 2*151 args that were handled
  by the exp(x)/2 case.  This leaves accuracy problems for about 2*19997
  args near the overflow threshold (~89); the maximum error there is
  2.5029 ulps.
  
  There are also accuracy probles for args in +-[0.5*ln2, 9] -- 2*188885
  args with errors of more than 1 ulp, with a maximum error of 1.384 ulps.
  
  Fixed a syntax error and naming errors in pseudo-code in comments.
  
  Revision  Changes    Path
  1.7       +6 -6      src/lib/msun/src/e_coshf.c