From owner-freebsd-bugs@FreeBSD.ORG  Sun Jul 29 02:31:27 2012
Return-Path: <owner-freebsd-bugs@FreeBSD.ORG>
Delivered-To: freebsd-bugs@FreeBSD.org
Received: from mx1.freebsd.org (mx1.freebsd.org [69.147.83.52])
	by hub.freebsd.org (Postfix) with ESMTP id 9A21C106566B
	for <freebsd-bugs@FreeBSD.org>; Sun, 29 Jul 2012 02:31:27 +0000 (UTC)
	(envelope-from brde@optusnet.com.au)
Received: from mail17.syd.optusnet.com.au (mail17.syd.optusnet.com.au
	[211.29.132.198])
	by mx1.freebsd.org (Postfix) with ESMTP id 19B098FC1A
	for <freebsd-bugs@FreeBSD.org>; Sun, 29 Jul 2012 02:31:26 +0000 (UTC)
Received: from c122-106-171-246.carlnfd1.nsw.optusnet.com.au
	(c122-106-171-246.carlnfd1.nsw.optusnet.com.au [122.106.171.246])
	by mail17.syd.optusnet.com.au (8.13.1/8.13.1) with ESMTP id
	q6T2VIhb016913
	(version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=NO);
	Sun, 29 Jul 2012 12:31:19 +1000
Date: Sun, 29 Jul 2012 12:31:18 +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: <201207281550.q6SFoBdC073014@freefall.freebsd.org>
Message-ID: <20120729121443.M1008@besplex.bde.org>
References: <201207281550.q6SFoBdC073014@freefall.freebsd.org>
MIME-Version: 1.0
Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed
Cc: freebsd-bugs@FreeBSD.org
Subject: Re: bin/170206: complex arcsinh, log, etc.
X-BeenThere: freebsd-bugs@freebsd.org
X-Mailman-Version: 2.1.5
Precedence: list
List-Id: Bug reports <freebsd-bugs.freebsd.org>
List-Unsubscribe: <http://lists.freebsd.org/mailman/listinfo/freebsd-bugs>,
	<mailto:freebsd-bugs-request@freebsd.org?subject=unsubscribe>
List-Archive: <http://lists.freebsd.org/pipermail/freebsd-bugs>
List-Post: <mailto:freebsd-bugs@freebsd.org>
List-Help: <mailto:freebsd-bugs-request@freebsd.org?subject=help>
List-Subscribe: <http://lists.freebsd.org/mailman/listinfo/freebsd-bugs>,
	<mailto:freebsd-bugs-request@freebsd.org?subject=subscribe>
X-List-Received-Date: Sun, 29 Jul 2012 02:31:27 -0000

On Sat, 28 Jul 2012, Stephen Montgomery-Smith wrote:

> OK.  This clog really seems to work.
>
> x*x + y*y - 1 is computed with a ULP less than 0.8.  The rest of the
> errors seem to be due to the implementation of log1p.  The ULP of the
> final answer seems to be never bigger than a little over 2.

I really don't like this version.  If we want bignums and slowness, we
can use lib[g]mp or even exec a mostly-interactive calculator (bc for
portability, closed-source bigware for unportability).

This version would have been useful in development to verify that the
doubled or tripled floating point was as exact as intended.  Now writing
it in your favourite interactive calculator is easiest, except for the
problem of comparing the results.

Bruce