From owner-freebsd-numerics@FreeBSD.ORG  Sat Sep 22 21:47:38 2012
Return-Path: <owner-freebsd-numerics@FreeBSD.ORG>
Delivered-To: freebsd-numerics@freebsd.org
Received: from mx1.freebsd.org (mx1.freebsd.org [69.147.83.52])
	by hub.freebsd.org (Postfix) with ESMTP id CC78A106566B
	for <freebsd-numerics@freebsd.org>;
	Sat, 22 Sep 2012 21:47:38 +0000 (UTC)
	(envelope-from brde@optusnet.com.au)
Received: from mail01.syd.optusnet.com.au (mail01.syd.optusnet.com.au
	[211.29.132.182])
	by mx1.freebsd.org (Postfix) with ESMTP id 5BF5F8FC14
	for <freebsd-numerics@freebsd.org>;
	Sat, 22 Sep 2012 21:47:37 +0000 (UTC)
Received: from c122-106-157-84.carlnfd1.nsw.optusnet.com.au
	(c122-106-157-84.carlnfd1.nsw.optusnet.com.au [122.106.157.84])
	by mail01.syd.optusnet.com.au (8.13.1/8.13.1) with ESMTP id
	q8MLlTQc012562
	(version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=NO);
	Sun, 23 Sep 2012 07:47:30 +1000
Date: Sun, 23 Sep 2012 07:47:29 +1000 (EST)
From: Bruce Evans <brde@optusnet.com.au>
X-X-Sender: bde@besplex.bde.org
To: Stephen Montgomery-Smith <stephen@missouri.edu>
In-Reply-To: <505E27CE.3060107@missouri.edu>
Message-ID: <20120923073807.K2059@besplex.bde.org>
References: <5017111E.6060003@missouri.edu> <5055ECA8.2080008@missouri.edu>
	<20120917022614.R2943@besplex.bde.org> <50562213.9020400@missouri.edu>
	<20120917060116.G3825@besplex.bde.org> <50563C57.60806@missouri.edu>
	<20120918012459.V5094@besplex.bde.org> <5057A932.3000603@missouri.edu>
	<5057F24B.7020605@missouri.edu> <20120918162105.U991@besplex.bde.org>
	<20120918232850.N2144@besplex.bde.org>
	<20120919010613.T2493@besplex.bde.org>
	<505BD9B4.8020801@missouri.edu> <20120921172402.W945@besplex.bde.org>
	<20120921212525.W1732@besplex.bde.org> <505C7490.90600@missouri.edu>
	<20120922042112.E3044@besplex.bde.org> <505CBF14.70908@missouri.edu>
	<505CC11A.5030502@missouri.edu> <20120922081607.F3613@besplex.bde.org>
	<20120922091625.Y3828@besplex.b! de.org>
	<505D1037.8010202@missouri.edu>
	<20120922142349.X4599@besplex.bde.org>
	<20120923044814.S1465@besplex.bde.org>
	<505E2575.6030302@missouri.edu> <505E27CE.3060107@missouri.edu>
MIME-Version: 1.0
Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed
Cc: freebsd-numerics@freebsd.org
Subject: Re: Complex arg-trig functions
X-BeenThere: freebsd-numerics@freebsd.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: "Discussions of high quality implementation of libm functions."
	<freebsd-numerics.freebsd.org>
List-Unsubscribe: <http://lists.freebsd.org/mailman/listinfo/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: Sat, 22 Sep 2012 21:47:38 -0000

On Sat, 22 Sep 2012, Stephen Montgomery-Smith wrote:

> 1. Your recent optimizations seem to have given an overall 3% time saving in 
> my timing tests.  That's pretty good in my opinion.

Hopefully more for large and small args :-).

> 2.  In my accuracy tests for casin(h), I have never seen the double or long 
> double have an error greater than 4 ULP.  But for the float case I have seen 
> 4.15 ULP.

I haven't seen any larger than 3.4.  What is the worst case you found?
Errors found for float precision tend to be because the density of bad
cases is higher so it is easier to test more of them accidentally.  I
did do some non-random testing for all float cases in narrow strips
about x or y = 0 or 1, but not for all combinations of this with all
functions.

Bruce