From owner-freebsd-numerics@FreeBSD.ORG Sun Aug 12 23:02:26 2012 Return-Path: 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 84CDB1065676 for ; Sun, 12 Aug 2012 23:02:26 +0000 (UTC) (envelope-from peter@rulingia.com) Received: from vps.rulingia.com (host-122-100-2-194.octopus.com.au [122.100.2.194]) by mx1.freebsd.org (Postfix) with ESMTP id 1D8BE8FC0A for ; Sun, 12 Aug 2012 23:02:25 +0000 (UTC) Received: from server.rulingia.com (c220-239-249-137.belrs5.nsw.optusnet.com.au [220.239.249.137]) by vps.rulingia.com (8.14.5/8.14.5) with ESMTP id q7CN2PDT075606 (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=OK) for ; Mon, 13 Aug 2012 09:02:26 +1000 (EST) (envelope-from peter@rulingia.com) X-Bogosity: Ham, spamicity=0.000000 Received: from server.rulingia.com (localhost.rulingia.com [127.0.0.1]) by server.rulingia.com (8.14.5/8.14.5) with ESMTP id q7CN2JI8021121 (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=NO) for ; Mon, 13 Aug 2012 09:02:19 +1000 (EST) (envelope-from peter@server.rulingia.com) Received: (from peter@localhost) by server.rulingia.com (8.14.5/8.14.5/Submit) id q7CN2JhM021120 for freebsd-numerics@freebsd.org; Mon, 13 Aug 2012 09:02:19 +1000 (EST) (envelope-from peter) Resent-From: Peter Jeremy Resent-Date: Mon, 13 Aug 2012 09:02:19 +1000 Resent-Message-ID: <20120812230219.GB20453@server.rulingia.com> Resent-To: freebsd-numerics@freebsd.org Received: from vps.rulingia.com (host-122-100-2-194.octopus.com.au [122.100.2.194]) by server.rulingia.com (8.14.5/8.14.5) with ESMTP id q6O2f28V064779 (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=OK) for ; Tue, 24 Jul 2012 12:41:03 +1000 (EST) (envelope-from stephen@missouri.edu) Received: from wilberforce.math.missouri.edu (wilberforce.math.missouri.edu [128.206.184.213]) by vps.rulingia.com (8.14.5/8.14.5) with ESMTP id q6O2exHq017785 (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=NO) for ; Tue, 24 Jul 2012 12:41:02 +1000 (EST) (envelope-from stephen@missouri.edu) Received: from [127.0.0.1] (wilberforce.math.missouri.edu [128.206.184.213]) by wilberforce.math.missouri.edu (8.14.5/8.14.5) with ESMTP id q6O2ZgQU092684; Mon, 23 Jul 2012 21:35:42 -0500 (CDT) (envelope-from stephen@missouri.edu) Message-ID: <500E09FE.8020607@missouri.edu> From: Stephen Montgomery-Smith Mail-Followup-To: freebsd-numerics@freebsd.org User-Agent: Mozilla/5.0 (X11; Linux i686; rv:14.0) Gecko/20120714 Thunderbird/14.0 MIME-Version: 1.0 To: Bruce Evans References: <5004A5C7.1040405@missouri.edu> <5004DEA9.1050001@missouri.edu> <20120717040118.GA86840@troutmask.apl.washington.edu> <20120717042125.GF66913@server.rulingia.com> <20120717043848.GB87001@troutmask.apl.washington.edu> <20120717225328.GA86902@server.rulingia.com> <20120717232740.GA95026@troutmask.apl.washington.edu> <20120718001337.GA87817@server.rulingia.com> <20120718123627.D1575@besplex.bde.org> <20120722121219.GC73662@server.rulingia.com> <20120722220031.GA7791@server.rulingia.com> <20120723141319.P1189@besplex.bde.org> <500CD98E.9080103@missouri.edu> <500D6182.8010003@missouri.edu> <20120724100014.I934@besplex.bde.org> In-Reply-To: <20120724100014.I934@besplex.bde.org> Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit Cc: Diane Bruce , John Baldwin , David Chisnall , Bruce Evans , Steve Kargl , David Schultz , Peter Jeremy , Warner Losh Subject: Re: Use of C99 extra long double math functions after r236148 X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , Date: Sun, 12 Aug 2012 23:02:26 -0000 X-Original-Date: Mon, 23 Jul 2012 21:35:42 -0500 X-List-Received-Date: Sun, 12 Aug 2012 23:02:26 -0000 I am committing the software - as is - to you guys. I feel that I worked extremely hard to get an error less than about 3.8 ULP (from 10,000,000 test cases) for the complex arc-trig and arc-hyperbolic functions. My skill is writing numerical software that works, with an awareness of where various numerical issues might arise. I am very proud of the work I have done in the last few days. Now this software requires a different skill set, getting it to conform to various styles, looking for smart efficiencies to overcome defects in the compiler, and checking that NaN's work properly. This is a skill set which you guys have in far greater amounts than I have. And it is something that doesn't engage my interest. If you find problems, like unwanted overflows/underflows or numerical errors occurring in edge cases, then by all means get back to me, and I will attempt to fix it. Thus far, the only fault I have found so far is in catanh, where I forgot to set the "inexact flag" when computing "hp" and "hm" when |y| is very small. I will still continue working hard on getting the numerics in clog to work correctly. Stephen