From owner-freebsd-current@FreeBSD.ORG Mon May 28 23:44:45 2012 Return-Path: Delivered-To: freebsd-current@freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:4f8:fff6::34]) by hub.freebsd.org (Postfix) with ESMTP id 12ED9106566B for ; Mon, 28 May 2012 23:44:45 +0000 (UTC) (envelope-from stephen@missouri.edu) Received: from wilberforce.math.missouri.edu (wilberforce.math.missouri.edu [128.206.184.213]) by mx1.freebsd.org (Postfix) with ESMTP id C5A718FC17 for ; Mon, 28 May 2012 23:44:43 +0000 (UTC) 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 q4SNigTh048304; Mon, 28 May 2012 18:44:42 -0500 (CDT) (envelope-from stephen@missouri.edu) Message-ID: <4FC40DEA.8030703@missouri.edu> Date: Mon, 28 May 2012 18:44:42 -0500 From: Stephen Montgomery-Smith User-Agent: Mozilla/5.0 (X11; Linux i686; rv:12.0) Gecko/20120430 Thunderbird/12.0.1 MIME-Version: 1.0 To: Steve Kargl References: <4FC30090.4070003@gwdg.de> <4D8CF7D2-CBEE-438E-A9E7-9C47A8892622@FreeBSD.org> <4FC36FE1.9080908@gwdg.de> <4FC38B81.6000302@gwdg.de> <4FC3A154.8030702@missouri.edu> <20120528203159.GA76340@troutmask.apl.washington.edu> <4FC3EBDA.2080502@missouri.edu> <20120528221731.GA76723@troutmask.apl.washington.edu> <4FC40449.3040602@missouri.edu> <20120528233035.GA77157@troutmask.apl.washington.edu> In-Reply-To: <20120528233035.GA77157@troutmask.apl.washington.edu> Content-Type: text/plain; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: 7bit Cc: freebsd-current@freebsd.org Subject: Re: Use of C99 extra long double math functions after r236148 X-BeenThere: freebsd-current@freebsd.org X-Mailman-Version: 2.1.5 Precedence: list List-Id: Discussions about the use of FreeBSD-current List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Mon, 28 May 2012 23:44:45 -0000 On 05/28/2012 06:30 PM, Steve Kargl wrote: > >> From clog.c in http://www.netlib.org/cephes/c9x-complex > > double complex > ccosh (z) > double complex z; > { > double complex w; > double x, y; > > x = creal(z); > y = cimag(z); > w = cosh (x) * cos (y) + (sinh (x) * sin (y)) * I; > return (w); > } > > See math_private.h about the above. > I looked in math_private.h - I presume you meant lib/msun/src/math_private.h. I wasn't able to find anything about ccosh there. I think that for a rough and ready ccosh, this is high enough quality for a math/cephes port. I do agree that it might not be high enough quality to make FreeBSD base. (Although I don't think I have ever been in a situation where I would have been tripped up by a transcendental function that responded incorrectly to exceptional input.) > And, finally, Yes, it is very nice. > > Who's writing the code to test the implementations? That is > better much the problem. Without testing, one might get an > implementation that appears to work until it doesn't! It took > me 3+ years to get sqrtl() into libm, but bde and das (and > myself) wanted to make sure the code worked. Fair enough if we are talking about the base system. > I haven't looked at glibc code in years, because I hack on libm > when I can. I do not want to run into questions about whether > my code is tainted by the gpl. > They had similar lists of exceptions.