From owner-freebsd-numerics@FreeBSD.ORG Sun Aug 12 23:04:36 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 E4AC01065673 for ; Sun, 12 Aug 2012 23:04:36 +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 69F628FC17 for ; Sun, 12 Aug 2012 23:04:36 +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 q7CN4atI075670 (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=OK) for ; Mon, 13 Aug 2012 09:04:36 +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 q7CN4TVP021260 (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=NO) for ; Mon, 13 Aug 2012 09:04:29 +1000 (EST) (envelope-from peter@server.rulingia.com) Received: (from peter@localhost) by server.rulingia.com (8.14.5/8.14.5/Submit) id q7CN4T6h021259 for freebsd-numerics@freebsd.org; Mon, 13 Aug 2012 09:04:29 +1000 (EST) (envelope-from peter) Resent-From: Peter Jeremy Resent-Date: Mon, 13 Aug 2012 09:04:29 +1000 Resent-Message-ID: <20120812230429.GP20453@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 q6SLxkGG072404 (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=OK) for ; Sun, 29 Jul 2012 07:59:47 +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 q6SLxiSU072407 (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA bits=256 verify=NO) for ; Sun, 29 Jul 2012 07:59:46 +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 q6SLxN3m068529; Sat, 28 Jul 2012 16:59:24 -0500 (CDT) (envelope-from stephen@missouri.edu) Message-ID: <501460BB.30806@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: Peter Jeremy References: <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> <500DAD41.5030104@missouri.edu> <20120724113214.G934@besplex.bde.org> <501204AD.30605@missouri.edu> <20120727032611.GB25690@server.rulingia.com> <20120728125824.GA26553@server.rulingia.com> In-Reply-To: <20120728125824.GA26553@server.rulingia.com> Content-Type: text/plain; charset=UTF-8; format=flowed Content-Transfer-Encoding: 8bit Cc: Diane Bruce , Bruce Evans , John Baldwin , David Chisnall , Bruce Evans , Steve Kargl , David Schultz , 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:04:37 -0000 X-Original-Date: Sat, 28 Jul 2012 16:59:23 -0500 X-List-Received-Date: Sun, 12 Aug 2012 23:04:37 -0000 On 07/28/2012 07:58 AM, Peter Jeremy wrote: > On 2012-Jul-27 13:26:11 +1000, Peter Jeremy wrote: >> I've been writing a test harness to vet the special case handling of >> all the complex functions (excluding cpow so far). Basically, it's >> just Appendix G.6 of WG14/N1256 turned into a C array, plus code to >> actually run the tests & interpret the results. So far, it's about >> 1100 lines of which about 1/3 is the test cases and is intended to run >> on x86/armle/sparc and FreeBSD/Linux/Solaris (I'm using Solaris and, >> to a lesser extent, Linux as a cross-check on my interpretation of the >> text). Once I'm happy with it, I'll circulate it. I was initially >> hoping to make it commitable but 8-char tabs and 80-char lines would >> require lots of line wrapping that would make it harder for me to >> follow. > > My test harness can be found at http://www.rulingia.com/~peter/ctest.c > There are no special compilation options, it just needs to be linked > with '-lm' (and '-ldl' on Linux). For normal use, just run the > executable - it will report any failures. For "finite" arguments, it > currently uses 3π/4 and 32769 other random numbers (the latter is > S_COUNT+1). > > It has two test modes for internal testing and debugging: > '-v' verifies that all the argument & result strings are valid and > that there's no duplication of argument vectors (for this purpose, it > doesn't consider '0' as finite and will incorrectly report '1' as an > invalid argument). > '-r' prints all the double-precision test vectors used. This should > generate 3604951 lines of output. > > The output should be reasonably self-explanatory except: > - double-precision function names are printed with a trailing 'd' > - an expected sign of '?' means "don't care". > > It reports no errors on OpenSolaris but does report a number of what > appear to be valid errors on Linux. > > Whilst I was debugging the code, I found the following elisp useful > for post-processing the output: > > (progn (downcase-region (point-min) (point-max)) > (repl-regexp "^ [ ]c" "..c") > (repl-regexp "^ c" ".c") > (repl-regexp " *0x[0-9a-f]+ *" " ") > (repl-regexp " *0x[0-9a-f]+$" "") > (repl-regexp "infinit[y]" "inf") > (repl-regexp "0\\.0+e\\+0+\\>" "zer") > (repl-regexp "1\\.0+e\\+0+\\>" "one") > (repl-regexp "3\\.14159[0-9]+e\\+00" "pi.") > (repl-regexp "1\\.57079[0-9]+e\\+00" "p_2") > (repl-regexp "7\\.85398[0-9]+e\\-01" "p_4") > (repl-regexp "2\\.35619[0-9]+e\\+00" "3p4") > (repl-regexp "[0-9]\\.[0-9]+e[-\\+][0-9]+" "fin") > (repl-regexp "^ *\012" "") > (repl-regexp "\012 *=" " =") > (repl-regexp "\012 *expected: *" " # ") > (repl-regexp "\012 *want *" " # ") > (repl-regexp " +" " ") > (repl-regexp "-\\+" " ") > (repl-regexp " +$" "") > (repl-regexp "\\([^)]\\)$" "\\1 %%") > (repl-regexp "^\\(.*= \\)\\(.\\)\\(...\\)\\( .*# \\)\\(.\\)\\3\\(.*\\) %%" "\\1\\2\\3\\4\\5\\3\\6 \\2\\5") > (repl-regexp "\\([^)]\\)$" "\\1 %%") > (repl-regexp "^\\(.*= .... \\)\\(.\\)\\(...\\)\\( # .... \\)\\(.\\)\\3\\(.*\\) .." "\\1\\2\\3\\4\\5\\3\\6 \\2\\5") > (repl-regexp "^\\(......\\)f\\(:.*\012\\)\\1d\\2\\1l\\2" "\\1x\\2") > (repl-regexp "^\\(......\\)d\\(:.*\012\\)\\1f\\2\\1l\\2" "\\1x\\2") > ) > > This turns the output into a series of lines like: > ..ctanx: +3p4 +inf = +zer +one # -zer +one +- ++ > fn ^ Argument Result Expected XX YY > +- precision (f/d/l) or 'x' if all 3 affected > > XX and YY are the real and imaginary actual and expected result signs > or '%' if the category differs between expected and actual. The above > line (from the Linux output) means that > ctan(3π/4 + I*Inf) returns (+0 + I*1) instead of (-0 + I*1) > '+-' means that the signs of the real parts differ > '++' means that the signs of the imaginary parts are both '+' > > Please let me know if you find any errors or have any comments. > It is a really nice program. I tried it on the clog and catrig functions. I was able to get the catrig functions to completely comply with your program. See the diff at the end of http://www.freebsd.org/cgi/query-pr.cgi?pr=bin/170206 The clog program was already working, because Bruce had fixed it up. I forgot - does it check the fenv settings as well? It would be great if it does.