From owner-freebsd-numerics@freebsd.org Mon Nov 6 19:49:58 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 23DE0E666DF for ; Mon, 6 Nov 2017 19:49:58 +0000 (UTC) (envelope-from mike.d.ft402@gmail.com) Received: from mail-wm0-x22d.google.com (mail-wm0-x22d.google.com [IPv6:2a00:1450:400c:c09::22d]) (using TLSv1.2 with cipher ECDHE-RSA-AES128-GCM-SHA256 (128/128 bits)) (Client CN "smtp.gmail.com", Issuer "Google Internet Authority G2" (verified OK)) by mx1.freebsd.org (Postfix) with ESMTPS id B5AF680F13 for ; Mon, 6 Nov 2017 19:49:57 +0000 (UTC) (envelope-from mike.d.ft402@gmail.com) Received: by mail-wm0-x22d.google.com with SMTP id t139so16625369wmt.1 for ; Mon, 06 Nov 2017 11:49:57 -0800 (PST) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20161025; h=from:date:to:subject:message-id:mime-version:content-disposition; bh=9QoVz8Y0aVWfJVeVRgz7MYUXBppqQuaPe9UUke29TVQ=; b=WYvY1YEazEJI7SC48ZjhaBY3p4o7HZU/nEnxd8AtN+AojiU2rspreCdEN72/hwQJdX liKiTWuvvCID6JQUjBnnecgD/GB2uS4s3MnKNg5FA4PWaQyJXcTDrdLv7l1KOXANAKkW qKagu18KCS2OCKLcVfZ7BHZ8u9oTA0tH5O+q8ImrEZUuW739M67tDfSEO04YEGeoRlJ+ pxVTOSDm2tdYNpUVCVVrxXztjoa7TN29rw0eCui09gyCSHrTJbYZNYJNSrw99AzzhIIN UJNjV4pGMpu2g3E2CnMZe7GMshPe7V0Dik1YNVuKGr9kQ/n2TnG0AeHtc9R5IO3900ck M0Lg== X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20161025; h=x-gm-message-state:from:date:to:subject:message-id:mime-version :content-disposition; bh=9QoVz8Y0aVWfJVeVRgz7MYUXBppqQuaPe9UUke29TVQ=; b=AcgRjvJkC0BP31tnvgoLsAUf3ZDFvCW4UQoR+QritO05OFSTVQ+71LvvUJxJ8UmJ7D 8bx4HQQh2jFTa53DaW30ig8JpOl/7Qllh8/qAnfn3wOKJt/KJlPBx7cp5nJAELYA/w3B KPOUKMcLow15grEYe5ipYndv6zyTlasPkUWygSH45/CYBaR7n/TcTds/S2hCIjvw4Blk uSb5jKb+D6a8KvrV48HQC4gGV7pFBP1HtnIKH4qpqNcLf/eAvFDoxsV6I9it5I9qUN8/ fZoIx76SDAZus38pZpyPD72bc+Y0EuJp3O0yYjAaPiOqEbAfZLwpNSg2uIpcLN1eEVXY zLmw== X-Gm-Message-State: AMCzsaUXruUCSaFfGEXazpTab5IVO2QP1UwNBBsgEXZIvFPFVXdrt1Yu SHXtWLrLcZi7XWATwxvdRqlTRA== X-Google-Smtp-Source: ABhQp+SK7LEhkH/y+lnPyuE971IMVOFevGv0Um+E2J58l5cRiMA2WGGro1MR6NlbkQfajk5++uTloA== X-Received: by 10.80.146.47 with SMTP id i44mr22154384eda.198.1509997795619; Mon, 06 Nov 2017 11:49:55 -0800 (PST) Received: from freebird.localdomain (uhh-wlan-fo-134-100-17-28.rrz.uni-hamburg.de. [134.100.17.28]) by smtp.gmail.com with ESMTPSA id e24sm10802722edc.9.2017.11.06.11.49.54 for (version=TLS1_2 cipher=ECDHE-RSA-AES128-GCM-SHA256 bits=128/128); Mon, 06 Nov 2017 11:49:55 -0800 (PST) From: Michael Danilov X-Google-Original-From: Michael Danilov Received: by freebird.localdomain (Postfix, from userid 1001) id CF6F22CC718F; Mon, 6 Nov 2017 20:49:53 +0100 (CET) Date: Mon, 6 Nov 2017 20:49:43 +0100 To: freebsd-numerics@freebsd.org Subject: cpow and clog Message-ID: <20171106194937.GA87725@freebird> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Mon, 06 Nov 2017 19:49:58 -0000 Hello, I would like to have some feedback on my attempt to import OpenBSD code for cpow and clog: https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=221341 https://bugs.freebsd.org/bugzilla/attachment.cgi?id=187693 What happened to the alternative implementation mentioned in the thread below? And what had stopped the developers from just reusing the Net- or OpenBSD code? Far as I see, FreeBSD's msun source resembles those other libm's a lot anyway? https://lists.freebsd.org/pipermail/freebsd-numerics/2012-August/000028.html Kind regards, Mike From owner-freebsd-numerics@freebsd.org Mon Nov 6 20:41:28 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id DB8F8E674C3 for ; Mon, 6 Nov 2017 20:41:28 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "troutmask", Issuer "troutmask" (not verified)) by mx1.freebsd.org (Postfix) with ESMTPS id C1E721585 for ; Mon, 6 Nov 2017 20:41:28 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (localhost [127.0.0.1]) by troutmask.apl.washington.edu (8.15.2/8.15.2) with ESMTPS id vA6KfMce037628 (version=TLSv1.2 cipher=DHE-RSA-AES256-GCM-SHA384 bits=256 verify=NO); Mon, 6 Nov 2017 12:41:22 -0800 (PST) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.15.2/8.15.2/Submit) id vA6KfLi0037627; Mon, 6 Nov 2017 12:41:21 -0800 (PST) (envelope-from sgk) Date: Mon, 6 Nov 2017 12:41:21 -0800 From: Steve Kargl To: Michael Danilov Cc: freebsd-numerics@freebsd.org Subject: Re: cpow and clog Message-ID: <20171106204121.GB37361@troutmask.apl.washington.edu> Reply-To: sgk@troutmask.apl.washington.edu References: <20171106194937.GA87725@freebird> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: <20171106194937.GA87725@freebird> User-Agent: Mutt/1.7.2 (2016-11-26) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Mon, 06 Nov 2017 20:41:29 -0000 On Mon, Nov 06, 2017 at 08:49:43PM +0100, Michael Danilov wrote: > Hello, > > I would like to have some feedback on my attempt to import OpenBSD > code for cpow and clog: > > https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=221341 > https://bugs.freebsd.org/bugzilla/attachment.cgi?id=187693 > > What happened to the alternative implementation mentioned in the thread below? bde has an implementation of clog[fl]. He may someday commit it. I don't know if anyone ever worked on cpow[fl]. I stopped working on powl and tgammal when I returned my commit bit due to differences with "higher-ranking" committers. > And what had stopped the developers from just reusing the Net-i > or OpenBSD code? How have you tested the NetBSD and/or OpenBSD code? What is the quality? Have the long double clogl and cpowl been tested on both ld80 and ld128 hardware? See FreeBSD's lib/msun/src/math_private.h for a discussion of possible issues of using I from complex.h in this code. > Far as I see, FreeBSD's msun source resembles those other libm's a > lot anyway? Well, of course, the various libm's look alike. They all started life from Sun Microsystems fdlibm code. You can get the last release from Netlib. http://www.netlib.org/fdlibm/index.html -- Steve From owner-freebsd-numerics@freebsd.org Mon Nov 6 21:33:09 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 6F5C1E67E21 for ; Mon, 6 Nov 2017 21:33:09 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "troutmask", Issuer "troutmask" (not verified)) by mx1.freebsd.org (Postfix) with ESMTPS id 48B1E2A4A for ; Mon, 6 Nov 2017 21:33:09 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (localhost [127.0.0.1]) by troutmask.apl.washington.edu (8.15.2/8.15.2) with ESMTPS id vA6LX8Wg073810 (version=TLSv1.2 cipher=DHE-RSA-AES256-GCM-SHA384 bits=256 verify=NO); Mon, 6 Nov 2017 13:33:08 -0800 (PST) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.15.2/8.15.2/Submit) id vA6LX8Hh073809; Mon, 6 Nov 2017 13:33:08 -0800 (PST) (envelope-from sgk) Date: Mon, 6 Nov 2017 13:33:08 -0800 From: Steve Kargl To: Michael Danilov Cc: freebsd-numerics@freebsd.org Subject: Re: cpow and clog Message-ID: <20171106213308.GA73787@troutmask.apl.washington.edu> Reply-To: sgk@troutmask.apl.washington.edu References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: <20171106204121.GB37361@troutmask.apl.washington.edu> User-Agent: Mutt/1.7.2 (2016-11-26) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Mon, 06 Nov 2017 21:33:09 -0000 On Mon, Nov 06, 2017 at 12:41:21PM -0800, Steve Kargl wrote: > On Mon, Nov 06, 2017 at 08:49:43PM +0100, Michael Danilov wrote: > > Hello, > > > > I would like to have some feedback on my attempt to import OpenBSD > > code for cpow and clog: > > > > https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=221341 > > https://bugs.freebsd.org/bugzilla/attachment.cgi?id=187693 > > > > What happened to the alternative implementation mentioned in the thread below? > > bde has an implementation of clog[fl]. He may someday > commit it. Note, I packaged Bruce's code with manpages here https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=216863 -- Steve 20170425 https://www.youtube.com/watch?v=VWUpyCsUKR4 20161221 https://www.youtube.com/watch?v=IbCHE-hONow From owner-freebsd-numerics@freebsd.org Tue Nov 7 10:39:04 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 778B7E54E05 for ; Tue, 7 Nov 2017 10:39:04 +0000 (UTC) (envelope-from dimpase@gmail.com) Received: from mail-wm0-x232.google.com (mail-wm0-x232.google.com [IPv6:2a00:1450:400c:c09::232]) (using TLSv1.2 with cipher ECDHE-RSA-AES128-GCM-SHA256 (128/128 bits)) (Client CN "smtp.gmail.com", Issuer "Google Internet Authority G2" (verified OK)) by mx1.freebsd.org (Postfix) with ESMTPS id 06567782DF for ; Tue, 7 Nov 2017 10:39:04 +0000 (UTC) (envelope-from dimpase@gmail.com) Received: by mail-wm0-x232.google.com with SMTP id r68so2768699wmr.3 for ; Tue, 07 Nov 2017 02:39:03 -0800 (PST) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20161025; h=from:date:to:cc:subject:message-id:references:mime-version :content-disposition:in-reply-to:user-agent; bh=7wTnPdG3S/gop3YBebaVo1rc2nrLcpSukWecaAxqj7A=; b=n/AJKfsMoQsTLyjO5K9ZEZS6kJQK1mrLaImoBdQjoXHJnZ5/3W3LBSZ4/nMcqalFyG Tnpib1/gq8wwQGKsnCGuqbmiHR0u+5uA/XL9evD+uIRmr6BdE+jIf9XPWuuR/hS9lEBH AXBz0l7Jdnp1LDudRjUgkR8PV93+ZLzxU7pyK8/lViYTmaHgl300wDFcax2udXuV9P6s gPLkN0+8TwNgYsDb/m9t0Awn4skhT4Wteg369UC7O4Md1ml2L/WiTnQN8MxfHqdMzkMF ZE9oCblRyqUAWa5zZPt74BSBkEBYoPM70n4wPBhYtY/F4NoDyyLLV86hctmrXPws3EWe OA0A== X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20161025; h=x-gm-message-state:from:date:to:cc:subject:message-id:references :mime-version:content-disposition:in-reply-to:user-agent; bh=7wTnPdG3S/gop3YBebaVo1rc2nrLcpSukWecaAxqj7A=; b=jWmOqNTPI3cj910Ji70MMnDTUK7D6lWBse7/oQLd5ao6FmVfhyB4iejVyHxslor0GP CHEpdNOVkg3GgGvEWbsSAzqeGdVfKdPVNh/W+9/wq/1P8mnSxJM//t4aqIvSVNOwHB7M 2wwzknJbwSrNVllKhfQvz3ACiwYPxBiskyEe0b7iR3+hxQy2L8LQxSnR0ntuux8i2N/l Fd4t5aNZ+d06tT7ELSBRrBYtInG4sqVk7VqTh1xvJM/0OOWAoH50koozgVOnBLewkAKh J/u64xrkyemf/HdtlRV5HhFPsDiIFH5mBijjYM8U/q0s34ucw4/fqsG/810hiltLaejX MOdg== X-Gm-Message-State: AJaThX5q0EfUdAY5SNgS5gyta/iXVml/Ll7OxQ4vcBaNFg3lWURwYpHu 2L5+isu0p0X3iUwhpND+6NbIiRuz X-Google-Smtp-Source: ABhQp+SFXNedc58c+JCTvMdn7PkGNo4XnEa6fUPM9n1vPCKvg8i0fRfO2g1MRt39TXjfPfbJeCh8tw== X-Received: by 10.28.5.148 with SMTP id 142mr1030147wmf.142.1510051142256; Tue, 07 Nov 2017 02:39:02 -0800 (PST) Received: from hilbert.cs.ox.ac.uk (cpc69402-oxfd27-2-0-cust392.4-3.cable.virginm.net. [82.14.197.137]) by smtp.gmail.com with ESMTPSA id 69sm2274422wmp.37.2017.11.07.02.39.00 (version=TLS1_2 cipher=ECDHE-RSA-AES128-GCM-SHA256 bits=128/128); Tue, 07 Nov 2017 02:39:00 -0800 (PST) From: dimpase@gmail.com X-Google-Original-From: dimpase+freebsd@gmail.com Received: by hilbert.cs.ox.ac.uk (Postfix, from userid 1000) id D73E174E96E; Tue, 7 Nov 2017 10:38:59 +0000 (GMT) Date: Tue, 7 Nov 2017 10:38:59 +0000 To: Steve Kargl Cc: Michael Danilov , freebsd-numerics@freebsd.org Subject: Re: cpow and clog Message-ID: <20171107103858.GA8468@hilbert> References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: <20171106204121.GB37361@troutmask.apl.washington.edu> X-PGP-Key: http://keys.gnupg.net/pks/lookup?op=vindex&search=dimpase@cs.ox.ac.uk User-Agent: Mutt/1.8.3 (2017-05-23) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 07 Nov 2017 10:39:04 -0000 On Mon, Nov 06, 2017 at 12:41:21PM -0800, Steve Kargl wrote: > On Mon, Nov 06, 2017 at 08:49:43PM +0100, Michael Danilov wrote: > > I would like to have some feedback on my attempt to import OpenBSD > > code for cpow and clog: > > > > https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=221341 > > https://bugs.freebsd.org/bugzilla/attachment.cgi?id=187693 > > > > What happened to the alternative implementation mentioned in the thread below? > > bde has an implementation of clog[fl]. He may someday > commit it. I don't know if anyone ever worked on cpow[fl]. > I stopped working on powl and tgammal when I returned my > commit bit due to differences with "higher-ranking" committers. > > > And what had stopped the developers from just reusing the Net-i > > or OpenBSD code? > > How have you tested the NetBSD and/or OpenBSD code? What is the > quality? Have the long double clogl and cpowl been tested on both > ld80 and ld128 hardware? See FreeBSD's lib/msun/src/math_private.h > for a discussion of possible issues of using I from complex.h in this > code. I would like to point out that various FreeBSD ports already contain implementations of the functions in question. For instance, the current numpy includes these implementations - they say it's taken from msun's FreeBSD back in 2013, mostly (sic!). https://github.com/numpy/numpy/blob/master/numpy/core/src/npymath/npy_math_complex.c.src Another such port, semi-obsolete, is sage (a.k.a. sagemath), something that I am keen on reviving. I am sure there are more such ports (e.g. in the C++-land), surely you ought to know better. Sorry for being blunt, but IMHO the attitude on this list appears to be to let the numerics stack on FreeBSD die a slow death. Indeed, most people hate to reinvent the wheel. It's really no fun at all to scramble to get these missing implementations somehow, there are certainly much better ways to use one's time and brainpower. On this list people prefer to point at some private code in uncertain shape, and hope that somehow by some magic FreeBSD will have the best humanely possible implementation of the complex transcendental functions... Why don't you first of all try to provide *some* reasonably working implementation (thus allowing porters not to have to reinvent this wheel, badly, for $n$-th time over, and then having *fun* making sure the tools know where to get these functions), and only then try to improve it? Cheers, Dima http://users.ox.ac.uk/~coml0531/ > > > Far as I see, FreeBSD's msun source resembles those other libm's a > > lot anyway? > > Well, of course, the various libm's look alike. They all started > life from Sun Microsystems fdlibm code. You can get the last > release from Netlib. > > http://www.netlib.org/fdlibm/index.html > > -- > Steve > _______________________________________________ > freebsd-numerics@freebsd.org mailing list > https://lists.freebsd.org/mailman/listinfo/freebsd-numerics > To unsubscribe, send any mail to "freebsd-numerics-unsubscribe@freebsd.org" From owner-freebsd-numerics@freebsd.org Tue Nov 7 15:31:50 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 72B65E5BCD5 for ; Tue, 7 Nov 2017 15:31:50 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "troutmask", Issuer "troutmask" (not verified)) by mx1.freebsd.org (Postfix) with ESMTPS id 52BED6EF71 for ; Tue, 7 Nov 2017 15:31:50 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (localhost [127.0.0.1]) by troutmask.apl.washington.edu (8.15.2/8.15.2) with ESMTPS id vA7FVnXs044844 (version=TLSv1.2 cipher=DHE-RSA-AES256-GCM-SHA384 bits=256 verify=NO); Tue, 7 Nov 2017 07:31:49 -0800 (PST) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.15.2/8.15.2/Submit) id vA7FVnXW044843; Tue, 7 Nov 2017 07:31:49 -0800 (PST) (envelope-from sgk) Date: Tue, 7 Nov 2017 07:31:49 -0800 From: Steve Kargl To: dimpase@gmail.com Cc: Michael Danilov , freebsd-numerics@freebsd.org Subject: Re: cpow and clog Message-ID: <20171107153149.GA44102@troutmask.apl.washington.edu> Reply-To: sgk@troutmask.apl.washington.edu References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171107103858.GA8468@hilbert> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: <20171107103858.GA8468@hilbert> User-Agent: Mutt/1.7.2 (2016-11-26) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 07 Nov 2017 15:31:50 -0000 On Tue, Nov 07, 2017 at 10:38:59AM +0000, dimpase@gmail.com wrote: > On Mon, Nov 06, 2017 at 12:41:21PM -0800, Steve Kargl wrote: > > On Mon, Nov 06, 2017 at 08:49:43PM +0100, Michael Danilov wrote: > > > I would like to have some feedback on my attempt to import OpenBSD > > > code for cpow and clog: > > > > > > https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=221341 > > > https://bugs.freebsd.org/bugzilla/attachment.cgi?id=187693 > > > > > > What happened to the alternative implementation mentioned in the thread below? > > > > bde has an implementation of clog[fl]. He may someday > > commit it. I don't know if anyone ever worked on cpow[fl]. > > I stopped working on powl and tgammal when I returned my > > commit bit due to differences with "higher-ranking" committers. > > > > > And what had stopped the developers from just reusing the Net-i > > > or OpenBSD code? > > > > How have you tested the NetBSD and/or OpenBSD code? What is the > > quality? Have the long double clogl and cpowl been tested on both > > ld80 and ld128 hardware? See FreeBSD's lib/msun/src/math_private.h > > for a discussion of possible issues of using I from complex.h in this > > code. > > I would like to point out that various FreeBSD ports already contain > implementations of the functions in question. What is the quality of those implementations? Have you tested the code? Have you read the comment in math_private.h concerning the use of I in computations? > Sorry for being blunt, but IMHO the attitude on this list appears to be > to let the numerics stack on FreeBSD die a slow death. You're talking to the one person who has spent decades trying to improve libm through bug fixes and implementing missing functionality. > Indeed, most people hate to reinvent the wheel. It's > really no fun at all to scramble to get these missing implementations > somehow, there are certainly much better ways to use one's time and > brainpower. On this list people prefer to point at some private code in > uncertain shape, and hope that somehow by some magic FreeBSD will have > the best humanely possible implementation of the complex transcendental > functions... On this list, which is mostly Bruce and me it seems, people care about the quality of the code. People, which is mostly Bruce and me, spend quite a bit of time benchmarking proposed patches go ensure that users aren't fed wrong results. > Why don't you first of all try to provide *some* reasonably > working implementation (thus allowing porters not to have to reinvent > this wheel, badly, for $n$-th time over, and then having *fun* making > sure the tools know where to get these functions), and only then try to > improve it? Who is 'you' here? I recently provided implementations of sinpi[fl], cospi[fl], and tanpi[fl] for ISO/IEC TS 18661-4. Patches were submitted to this list (check the archive). Bruce is the only person to comment on the patch. I took the time to fix the lang/python* ports, which break with this patch, and sent the patches to freebsd-ports list where the patches are lost in time. Check out the history of msun/src/imprecise.c. This is a 'workaround' committed by David Chisnall so that FreeBSD could claim C++ conformance. David has never lifted a finger to fix that kludge. I am the person who has clean up that mess (until I stopped working on powl and tgammal). This is exactly the the situation that people on this list are trying to avoid when someone wants to simply "provide *some* reasonably working implementation". Instead of complaining about the people on this list, which is mostly Bruce and me, why not help. So, I'll repeat myself here. What is the quality of those implementations? Have you tested the code? Have you read the comment in math_private.h concerning the use of I in computations? -- Steve From owner-freebsd-numerics@freebsd.org Tue Nov 7 17:54:01 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 34E93E5FFCE for ; Tue, 7 Nov 2017 17:54:01 +0000 (UTC) (envelope-from dimpase@gmail.com) Received: from mail-wm0-x236.google.com (mail-wm0-x236.google.com [IPv6:2a00:1450:400c:c09::236]) (using TLSv1.2 with cipher ECDHE-RSA-AES128-GCM-SHA256 (128/128 bits)) (Client CN "smtp.gmail.com", Issuer "Google Internet Authority G2" (verified OK)) by mx1.freebsd.org (Postfix) with ESMTPS id B1A027B5F8 for ; Tue, 7 Nov 2017 17:54:00 +0000 (UTC) (envelope-from dimpase@gmail.com) Received: by mail-wm0-x236.google.com with SMTP id n74so15263667wmi.1 for ; Tue, 07 Nov 2017 09:54:00 -0800 (PST) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20161025; h=from:date:to:cc:subject:message-id:references:mime-version :content-disposition:in-reply-to:user-agent; bh=rUMXPg/MEbuzBCGeMC+ZnqCHyYxS7vaASPYeyuBvBQE=; b=DHHeaK1Uquh2FGCsd7EXfuuKelrG5eaYvAAfy4tygkS5ReRh2KnZkitKK5cGGsQ22L BkGf3MegKoXcnaQeul4KcIZmX6Vg8E08JjMWB9qxLIDE00ZPWNH+ytyxvF6M+NjclOrw pIthWUuwgmDVZaR23TJWqGLeaWyulS0QsuGhGBWn7j357B6sAohdjjCfhN8/fV+QNbpQ KgcQ9YJqEHFuAX/HWTNBo3eVQxuPIaILWK8kQm3u8Y8lTmTQ9ZHtH3biXeFazZAd5cwl 8kSXA+6tUtXkA75x2bYqamav0qs/HrdZR9V1EGljyoh2SgpJyjGGbWq30lLk4Vbejsl2 oKlg== X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20161025; h=x-gm-message-state:from:date:to:cc:subject:message-id:references :mime-version:content-disposition:in-reply-to:user-agent; bh=rUMXPg/MEbuzBCGeMC+ZnqCHyYxS7vaASPYeyuBvBQE=; b=Tg5CHbQwUFBthlkcreDi+9K+aPL1hjiAIyYGS4yUKKVsCmplfym0cAmAEqneGsXA5R SHAag4B6EFqKGV9Kg13bmrxOt3Hu63JI1of/hIqWxRg7nFmkcUrJo+lhhb7HvomZemRq aOYih0CTtdi3HDzWdeRbBbvZ+Usgl0L3op7DFw19AwK1JOFGXttF3A7oPq6cIqgP6y4t 36o2e/JwXZi1v13YF1OjFBPzQihqnIv+bK0fR+Nd+rycg1zBziSsJSkDy5ILUEfGRHk5 /jca4oRu+pCPIWRM6gP7KtlBmp5JdybRbG0z3thB8bgnqKFE1cmbs4TU7jFGNZ0g5YoA Hj6g== X-Gm-Message-State: AJaThX4ezZ9xeZv12NUwHWo+srUudraXe519pMkW+yDx7DrNQnMl7LXB etKRDq62LF1CXGXCBNx8OMw+hHpz X-Google-Smtp-Source: ABhQp+QYtLHnid7Gz02hMB8CM48/jNSuZ7LtzAYLk14bT2QQjGwzwn8vZc2n5+u4Rd6exGu53w2C6Q== X-Received: by 10.80.164.183 with SMTP id w52mr28590edb.44.1510077238622; Tue, 07 Nov 2017 09:53:58 -0800 (PST) Received: from hilbert.cs.ox.ac.uk (cpc69402-oxfd27-2-0-cust392.4-3.cable.virginm.net. [82.14.197.137]) by smtp.gmail.com with ESMTPSA id h51sm1947596eda.56.2017.11.07.09.53.57 (version=TLS1_2 cipher=ECDHE-RSA-AES128-GCM-SHA256 bits=128/128); Tue, 07 Nov 2017 09:53:57 -0800 (PST) From: dimpase@gmail.com X-Google-Original-From: dimpase+freebsd@gmail.com Received: by hilbert.cs.ox.ac.uk (Postfix, from userid 1000) id 809DC74E96E; Tue, 7 Nov 2017 17:53:56 +0000 (GMT) Date: Tue, 7 Nov 2017 17:53:56 +0000 To: Steve Kargl Cc: Michael Danilov , freebsd-numerics@freebsd.org Subject: Re: cpow and clog Message-ID: <20171107175354.GA26733@hilbert> References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171107103858.GA8468@hilbert> <20171107153149.GA44102@troutmask.apl.washington.edu> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: <20171107153149.GA44102@troutmask.apl.washington.edu> X-PGP-Key: http://keys.gnupg.net/pks/lookup?op=vindex&search=dimpase@cs.ox.ac.uk User-Agent: Mutt/1.8.3 (2017-05-23) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 07 Nov 2017 17:54:01 -0000 On Tue, Nov 07, 2017 at 07:31:49AM -0800, Steve Kargl wrote: > On Tue, Nov 07, 2017 at 10:38:59AM +0000, dimpase@gmail.com wrote: > > On Mon, Nov 06, 2017 at 12:41:21PM -0800, Steve Kargl wrote: > > > On Mon, Nov 06, 2017 at 08:49:43PM +0100, Michael Danilov wrote: > > > > I would like to have some feedback on my attempt to import OpenBSD > > > > code for cpow and clog: > > > > > > > > https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=221341 > > > > https://bugs.freebsd.org/bugzilla/attachment.cgi?id=187693 > > > > > > > > What happened to the alternative implementation mentioned in the thread below? > > > > > > bde has an implementation of clog[fl]. He may someday > > > commit it. I don't know if anyone ever worked on cpow[fl]. > > > I stopped working on powl and tgammal when I returned my > > > commit bit due to differences with "higher-ranking" committers. > > > > > > > And what had stopped the developers from just reusing the Net-i > > > > or OpenBSD code? > > > > > > How have you tested the NetBSD and/or OpenBSD code? What is the > > > quality? Have the long double clogl and cpowl been tested on both > > > ld80 and ld128 hardware? See FreeBSD's lib/msun/src/math_private.h > > > for a discussion of possible issues of using I from complex.h in this > > > code. > > > > I would like to point out that various FreeBSD ports already contain > > implementations of the functions in question. > > What is the quality of those implementations? Have you tested the > code? Do you have a testsuite for these functions? Perhaps it can be found somewhere, do you know? We are able to make bits of code needed for conformance to pass our tests in Sagemath on FreeBSD 11 (x86_64 only, although I probably can get FreeBSD running on aarch64 somehow) with various clang versions, and with gcc6. Basically I am talking about these functions from the Cephes library: http://files.sagemath.org/spkg/upstream/cephes/index.html > Have you read the comment in math_private.h concerning the > use of I in computations? I presume you refer to "...x+I*y is * currently unusable in general since gcc introduces many overflow, * underflow, sign and efficiency bugs..." We are building FreeBSD with clang, and not gcc, I don't know how relevant it is, for clang is a fast-moving target. By the way, https://wiki.freebsd.org/Numerics still mentiones FreeBSD 10 as "current"... > > > Sorry for being blunt, but IMHO the attitude on this list appears to be > > to let the numerics stack on FreeBSD die a slow death. > > You're talking to the one person who has spent decades trying to > improve libm through bug fixes and implementing missing functionality. > > > Indeed, most people hate to reinvent the wheel. It's > > really no fun at all to scramble to get these missing implementations > > somehow, there are certainly much better ways to use one's time and > > brainpower. On this list people prefer to point at some private code in > > uncertain shape, and hope that somehow by some magic FreeBSD will have > > the best humanely possible implementation of the complex transcendental > > functions... > > On this list, which is mostly Bruce and me it seems, people care about > the quality of the code. People, which is mostly Bruce and me, spend > quite a bit of time benchmarking proposed patches go ensure that users > aren't fed wrong results. > > > Why don't you first of all try to provide *some* reasonably > > working implementation (thus allowing porters not to have to reinvent > > this wheel, badly, for $n$-th time over, and then having *fun* making > > sure the tools know where to get these functions), and only then try to > > improve it? > > Who is 'you' here? I recently provided implementations of sinpi[fl], cospi[fl], > and tanpi[fl] for ISO/IEC TS 18661-4. Patches were submitted to this list > (check the archive). Well, this involves browsing entries of https://lists.freebsd.org/pipermail/freebsd-numerics/ Would you mind at least mention year and month? > Bruce is the only person to comment on the patch. I took > the time to fix the lang/python* ports, which break with this patch, and sent > the patches to freebsd-ports list where the patches are lost in time. Nowadays one would rather expect a git branch with the changes, not patches, as you might know. > > Check out the history of msun/src/imprecise.c. This is a 'workaround' > committed by David Chisnall so that FreeBSD could claim C++ conformance. > David has never lifted a finger to fix that kludge. I am the person > who has clean up that mess (until I stopped working on powl and tgammal). > This is exactly the the situation that people on this list are trying to > avoid when someone wants to simply "provide *some* reasonably working > implementation". See, it was great work, it has allowed him to publish stuff, e.g.: http://www.cl.cam.ac.uk/research/security/ctsrd/pdfs/201403-asiabsdcon2014-llvmbsd.pdf :-) His homepage however says that he's "Former Core Team member" of FreeBSD. If https://wiki.freebsd.org/Numerics is the right place to look at, this issue ought to be mentioned there, don't you think so? > > Instead of complaining about the people on this list, which is mostly > Bruce and me, why not help. So, I'll repeat myself here. What is the > quality of those implementations? Have you tested the code? Have you > read the comment in math_private.h concerning the use of I in computations? While I do have some time to spend on this in the coming year or so, I would like to understand how the developent process works in this particular case, who have rights to commit changes, how tests are run, whether there is any CI system to be used, etc. It seems that one has to start by setting up a set of tests for these functions. Dima > > -- > Steve From owner-freebsd-numerics@freebsd.org Tue Nov 7 19:24:57 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 83BE2E62CA1 for ; Tue, 7 Nov 2017 19:24:57 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "troutmask", Issuer "troutmask" (not verified)) by mx1.freebsd.org (Postfix) with ESMTPS id 626CF35AA for ; Tue, 7 Nov 2017 19:24:57 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (localhost [127.0.0.1]) by troutmask.apl.washington.edu (8.15.2/8.15.2) with ESMTPS id vA7JOuKL046728 (version=TLSv1.2 cipher=DHE-RSA-AES256-GCM-SHA384 bits=256 verify=NO); Tue, 7 Nov 2017 11:24:56 -0800 (PST) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.15.2/8.15.2/Submit) id vA7JOtsm046727; Tue, 7 Nov 2017 11:24:55 -0800 (PST) (envelope-from sgk) Date: Tue, 7 Nov 2017 11:24:55 -0800 From: Steve Kargl To: dimpase@gmail.com Cc: Michael Danilov , freebsd-numerics@freebsd.org Subject: Re: cpow and clog Message-ID: <20171107192455.GA46182@troutmask.apl.washington.edu> Reply-To: sgk@troutmask.apl.washington.edu References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171107103858.GA8468@hilbert> <20171107153149.GA44102@troutmask.apl.washington.edu> <20171107175354.GA26733@hilbert> MIME-Version: 1.0 Content-Type: text/plain; charset=utf-8 Content-Disposition: inline Content-Transfer-Encoding: 8bit In-Reply-To: <20171107175354.GA26733@hilbert> User-Agent: Mutt/1.7.2 (2016-11-26) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Tue, 07 Nov 2017 19:24:57 -0000 On Tue, Nov 07, 2017 at 05:53:56PM +0000, dimpase@gmail.com wrote: > On Tue, Nov 07, 2017 at 07:31:49AM -0800, Steve Kargl wrote: > > On Tue, Nov 07, 2017 at 10:38:59AM +0000, dimpase@gmail.com wrote: > > > On Mon, Nov 06, 2017 at 12:41:21PM -0800, Steve Kargl wrote: > > > > > > > > How have you tested the NetBSD and/or OpenBSD code? What is the > > > > quality? Have the long double clogl and cpowl been tested on both > > > > ld80 and ld128 hardware? See FreeBSD's lib/msun/src/math_private.h > > > > for a discussion of possible issues of using I from complex.h in this > > > > code. > > > > > > I would like to point out that various FreeBSD ports already contain > > > implementations of the functions in question. > > > > What is the quality of those implementations? Have you tested the > > code? > Do you have a testsuite for these functions? Not for these functions. That's why I'm asking if someone has tested them. I'm still working on non-complex functions, I can actually use in my *real* job. > Perhaps it can be found somewhere, do you know? You could ask the Netbsd and Openbsd projects if there are tests available or you could attempt to write a test. Google 'n1256.pdf'. Look at 7.3.7.2 and Annex G. clog has a branch cut along the negative real axis. If one is close to the branch cut, does clog return the correct answer or does it (occasionally) jump across the branch cut. clog needs to satisfy a number conditions for exceptional values (See Annex G). Does the NetBSd/OpenBSD code meet those conditions. > We are able to make bits of code needed for conformance to pass > our tests in Sagemath > on FreeBSD 11 (x86_64 only, although I probably can get FreeBSD running > on aarch64 somehow) with various clang versions, and with gcc6. > Basically I am talking about these functions from the Cephes library: > http://files.sagemath.org/spkg/upstream/cephes/index.html > > > Have you read the comment in math_private.h concerning the > > use of I in computations? > I presume you refer to > > "...x+I*y is > * currently unusable in general since gcc introduces many overflow, > * underflow, sign and efficiency bugs..." > > We are building FreeBSD with clang, and not gcc, I don't know how > relevant it is, for clang is a fast-moving target. Last time I checked, clang had the same issue as gcc. The issue was reported to clang developers 7 years ago. https://bugs.llvm.org/show_bug.cgi?id=8532 > By the way, https://wiki.freebsd.org/Numerics still mentiones FreeBSD 10 > as "current"... David Das use to keep the wiki up-to-date. I haven't seen das@ in a long time. Feel free to fix the wiki. > > > Sorry for being blunt, but IMHO the attitude on this list appears to be > > > to let the numerics stack on FreeBSD die a slow death. > > > > You're talking to the one person who has spent decades trying to > > improve libm through bug fixes and implementing missing functionality. > > > > > Indeed, most people hate to reinvent the wheel. It's > > > really no fun at all to scramble to get these missing implementations > > > somehow, there are certainly much better ways to use one's time and > > > brainpower. On this list people prefer to point at some private code in > > > uncertain shape, and hope that somehow by some magic FreeBSD will have > > > the best humanely possible implementation of the complex transcendental > > > functions... > > > > On this list, which is mostly Bruce and me it seems, people care about > > the quality of the code. People, which is mostly Bruce and me, spend > > quite a bit of time benchmarking proposed patches go ensure that users > > aren't fed wrong results. > > > > > Why don't you first of all try to provide *some* reasonably > > > working implementation (thus allowing porters not to have to reinvent > > > this wheel, badly, for $n$-th time over, and then having *fun* making > > > sure the tools know where to get these functions), and only then try to > > > improve it? > > > > Who is 'you' here? I recently provided implementations of sinpi[fl], > > cospi[fl], and tanpi[fl] for ISO/IEC TS 18661-4. Patches were submitted > > to this list (check the archive). > > Well, this involves browsing entries of > https://lists.freebsd.org/pipermail/freebsd-numerics/ > Would you mind at least mention year and month? Google "FreeBSD sinpi" the top two results are 218514 – [LIBM] implementations of sinpi[fl], cospi ... - FreeBSD Bugzilla https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=218514 [PATCH] lang/python27 -- Fix namespace collision - FreeBSD lists https://lists.freebsd.org/pipermail/freebsd-ports/2017-July/109675.html It is also not too difficult to click on the various *recent* months in the mail archive, but since you ask. https://lists.freebsd.org/pipermail/freebsd-numerics/2017-April/thread.html https://lists.freebsd.org/pipermail/freebsd-numerics/2017-May/thread.html > > Bruce is the only person to comment on the patch. I took > > the time to fix the lang/python* ports, which break with this > > patch, and sent the patches to freebsd-ports list where the > > patches are lost in time. > > Nowadays one would rather expect a git branch with the changes, not > patches, as you might know. I don't use git. cd /usr/ports svn diff lang/python27 > python27.diff svn diff lang/python36 > python36.diff Mail patches to freebsd-ports and freebsd-python. Watch patches wither. > > Check out the history of msun/src/imprecise.c. This is a 'workaround' > > committed by David Chisnall so that FreeBSD could claim C++ conformance. > > David has never lifted a finger to fix that kludge. I am the person > > who has clean up that mess (until I stopped working on powl and tgammal). > > This is exactly the the situation that people on this list are trying to > > avoid when someone wants to simply "provide *some* reasonably working > > implementation". > > See, it was great work, it has allowed him to publish stuff, e.g.: > http://www.cl.cam.ac.uk/research/security/ctsrd/pdfs/201403-asiabsdcon2014-llvmbsd.pdf > :-) > > His homepage however says that he's "Former Core Team member" of FreeBSD. > If https://wiki.freebsd.org/Numerics is the right place to look at, this > issue ought to be mentioned there, don't you think so? I objected when it was first suggested. However, I wasn't a member of Core, so you see were those objections went. Feel free to update the wiki. > > Instead of complaining about the people on this list, which is mostly > > Bruce and me, why not help. So, I'll repeat myself here. What is the > > quality of those implementations? Have you tested the code? Have you > > read the comment in math_private.h concerning the use of I in computations? > > While I do have some time to spend on this in the coming year or so, I > would like to understand how the developent process works in this > particular case, For libm stuff, the development process in the recent past has been Steve (or Stephen Montegomery or Peter Jeremy) writes code for a function. Steve does a significant amount of testing. Steve submits patch for freebsd-numerics. Bruce notices patch and grabs it. Bruce tests patch, fixes patch, tests fixed patch, and then sends Steve a long email with suggested improvement and/or issues with the original code. Steve understands about 10% of Bruce's email and tries to fix the code. Steve then does significant testing of the new code. Steve creates a new patch and submits it to FreeBSD-numerics. Rinse-and-repeat. It should also be noted that I no longer have a commit bit and Bruce rarely commits, so patch to libm lnger in bugzilla. > who have rights to commit changes, how tests are run, > whether there is any CI system to be used, etc. > It seems that one has to start by setting up a set of tests for these > functions. At a minimum grab n1256.pdf. Annex G has G.6.3.2 The clog functions * clog(conj(z)) = conj(clog(z)). * clog(-0 + i*0) returns -inf + i*pi and raises the divide-by-zero floating-point exception. * clog(+0 + i*0) returns -inf + i*0 and raises the divide-by-zero floating-point exception. * clog(x + i*inf) returns +inf + i*pi/2, for finite x. * clog(x + i*NaN) returns NaN + i*NaN and optionally raises the invalid floating-point exception, for finite x. * clog(-inf + i*y) returns +inf + i*pi, for finite positive-signed y. * clog(+inf + i*y) returns +inf + i*0, for finite positive-signed y. * clog(-inf + i*inf) returns +inf + i*3*pi/4. * clog(+inf + i*inf) returns +inf + i*pi/4. * clog(+-inf + i*NaN) returns +inf + i*NaN. * clog(NaN + iy) returns NaN + i*NaN and optionally raises the invalid floating-point exception, for finite y. * clog(NaN + i*inf) returns +inf + i*NaN. * clog(NaN + i*NaN) returns NaN + i*NaN. Do the NetBSD/OpenBSD meet the above conditions? For finite x+I*y, what is the accuracy of the return results? -- Steve From owner-freebsd-numerics@freebsd.org Wed Nov 8 01:50:44 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id BF2FAE6A987 for ; Wed, 8 Nov 2017 01:50:44 +0000 (UTC) (envelope-from stephen@missouri.edu) Received: from um-nip4-missouri-out.um.umsystem.edu (um-nip4-missouri-out.um.umsystem.edu [198.209.49.177]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "um-tip1.um.umsystem.edu", Issuer "InCommon RSA Server CA" (verified OK)) by mx1.freebsd.org (Postfix) with ESMTPS id 772C78036D for ; Wed, 8 Nov 2017 01:50:44 +0000 (UTC) (envelope-from stephen@missouri.edu) X-IronPort-Anti-Spam-Filtered: true X-IronPort-Anti-Spam-Result: =?us-ascii?q?A2H7AgA3YgJa/xY40cZdGwEBAQEDAQEBC?= =?us-ascii?q?QEBAYM0ZG4nB4N2mUWBViaYWgEJI4FegzoCGoRdQxQBAQEBAQEBAQEDaB0Lgmp?= =?us-ascii?q?HWAEBAQEBASMCDV4BBSMRVQIBCBgCAiYCAgIwFRACBA0IAQGKHxCoaYIniBcBB?= =?us-ascii?q?4J6AQEIAQEBAR8FgQ+CIYIHg2YLgnaDQoFTF4J+gmMFohYCh2SgT4xniQwCAgI?= =?us-ascii?q?CCQIagTk2IoFxfBNJgmWEXncBiwIBgRABAQE?= X-IPAS-Result: =?us-ascii?q?A2H7AgA3YgJa/xY40cZdGwEBAQEDAQEBCQEBAYM0ZG4nB4N?= =?us-ascii?q?2mUWBViaYWgEJI4FegzoCGoRdQxQBAQEBAQEBAQEDaB0LgmpHWAEBAQEBASMCD?= =?us-ascii?q?V4BBSMRVQIBCBgCAiYCAgIwFRACBA0IAQGKHxCoaYIniBcBB4J6AQEIAQEBAR8?= =?us-ascii?q?FgQ+CIYIHg2YLgnaDQoFTF4J+gmMFohYCh2SgT4xniQwCAgICCQIagTk2IoFxf?= =?us-ascii?q?BNJgmWEXncBiwIBgRABAQE?= Received: from ex2-n14.um.umsystem.edu ([198.209.56.22]) by um-nip4-exch-relay.um.umsystem.edu with ESMTP; 07 Nov 2017 19:49:33 -0600 Received: from EX2-N14.um.umsystem.edu (198.209.56.22) by EX2-N14.um.umsystem.edu (198.209.56.22) with Microsoft SMTP Server (version=TLS1_2, cipher=TLS_ECDHE_RSA_WITH_AES_128_CBC_SHA256_P256) id 15.1.1261.35; Tue, 7 Nov 2017 19:49:32 -0600 Received: from EX2-N14.um.umsystem.edu ([198.209.56.22]) by EX2-N14.um.umsystem.edu ([198.209.56.22]) with mapi id 15.01.1261.035; Tue, 7 Nov 2017 19:49:32 -0600 From: "Montgomery-Smith, Stephen" To: "freebsd-numerics@freebsd.org" Subject: Re: cpow and clog Thread-Topic: cpow and clog Thread-Index: AQHTVzhwY4EyR5x3Wk2DmT17BVFyU6MINVmAgAHoboA= Date: Wed, 8 Nov 2017 01:49:32 +0000 Message-ID: <4330a817-ae31-6336-7aa6-3e9981b02ab7@missouri.edu> References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> In-Reply-To: <20171106204121.GB37361@troutmask.apl.washington.edu> Accept-Language: en-US Content-Language: en-US X-MS-Has-Attach: X-MS-TNEF-Correlator: user-agent: Mozilla/5.0 (X11; Linux x86_64; rv:52.0) Gecko/20100101 Thunderbird/52.4.0 x-originating-ip: [128.206.49.160] Content-Type: text/plain; charset="utf-8" Content-ID: <77F2E25B9A66B341ACFE9E67457E717C@missouri.edu> Content-Transfer-Encoding: base64 MIME-Version: 1.0 X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 01:50:44 -0000 T24gMTEvMDYvMjAxNyAwMjo0MSBQTSwgU3RldmUgS2FyZ2wgd3JvdGU6DQo+IE9uIE1vbiwgTm92 IDA2LCAyMDE3IGF0IDA4OjQ5OjQzUE0gKzAxMDAsIE1pY2hhZWwgRGFuaWxvdiB3cm90ZToNCj4+ IEhlbGxvLA0KPj4NCj4+IEkgd291bGQgbGlrZSB0byBoYXZlIHNvbWUgZmVlZGJhY2sgb24gbXkg YXR0ZW1wdCB0byBpbXBvcnQgT3BlbkJTRA0KPj4gY29kZSBmb3IgY3BvdyBhbmQgY2xvZzoNCj4+ DQo+PiBodHRwczovL2J1Z3MuZnJlZWJzZC5vcmcvYnVnemlsbGEvc2hvd19idWcuY2dpP2lkPTIy MTM0MQ0KPj4gaHR0cHM6Ly9idWdzLmZyZWVic2Qub3JnL2J1Z3ppbGxhL2F0dGFjaG1lbnQuY2dp P2lkPTE4NzY5Mw0KDQpUaGUgcHJvYmxlbSB3aXRoIHRoaXMgaW1wbGVtZW50YXRpb24gb2YgY2xv ZyBpcyB0aGF0IHRoZXkgd2lsbCBwcm9kdWNlDQp2ZXJ5IGluYWNjdXJhdGUgcmVzdWx0cyB3aGVu IHx6fCBpcyBjbG9zZSB0byAxLiAgQSBwcm9wZXIgaW1wbGVtZW50YXRpb24NCmlzIGdpdmVuIGlu IHRoZSBwYXBlcjoNClQuIEUuIEh1bGwsIFRob21hcyBGLiBGYWlyZ3JpZXZlLCBQaW5nIFRhayBQ ZXRlciBUYW5nLCDigJxJbXBsZW1lbnRpbmcNCkNvbXBsZXggRWxlbWVudGFyeSBGdW5jdGlvbnMg VXNpbmcgRXhjZXB0aW9uIEhhbmRsaW5n4oCdLCBBQ00gVHJhbnNhY3Rpb25zDQpvbiBNYXRoZW1h dGljYWwgU29mdHdhcmUsIFZvbHVtZSAyMCAoMTk5NCksIHBwIDIxNeKAkzI0NCwgQ29ycmlnZW5k YSwgcDU1My4NCg0KVGhleSB1c2UgYSBwcmVjaXNpb24gd2hpY2ggaXMgdHdpY2UgdGhhdCBvZiB0 aGUgcmVxdWlyZWQgcHJlY2lzaW9uLiAgRm9yDQpleGFtcGxlLCBpZiB5b3Ugd2FudCBsb2coc3Fy dCh4XjIreV4yKSkgYWNjdXJhdGUgdG8gYWJvdXQgMTcgZGlnaXRzLCB5b3UNCndpbGwgbmVlZCB0 byBjb21wdXRlIGl0IHRvIDM0IGRpZ2l0cy4NCg0KSWYgYWxsIHlvdSBjYXJlIGFib3V0IGlzIGdl dHRpbmcgYSBhIGdvb2QgcmVsYXRpdmUgZXJyb3IgZm9yIHRoZSB0b3RhbA0KYW5zd2VyIChyYXRo ZXIgdGhhbiBpbnNpc3Rpbmcgb24gYSBnb29kIHJlbGF0aXZlIGVycm9yIGZvciB0aGUgcmVhbCBh bmQNCmltYWdpbmFyeSBwYXJ0cyBzZXBhcmF0ZWx5KSwgdGhlbiBmb3IgfHgrSSp5fCBjbG9zZSB0 byBvbmUsIGNvbXB1dGUgdGhlDQpyZWFsIHBhcnQgdXNpbmcNCg0KMC41KmxvZ3AxKCh4LTEpKih4 KzEpK3kqeSkNCg0KDQoNCj4+DQo+PiBXaGF0IGhhcHBlbmVkIHRvIHRoZSBhbHRlcm5hdGl2ZSBp bXBsZW1lbnRhdGlvbiBtZW50aW9uZWQgaW4gdGhlIHRocmVhZCBiZWxvdz8NCj4gDQo+IGJkZSBo YXMgYW4gaW1wbGVtZW50YXRpb24gb2YgY2xvZ1tmbF0uICBIZSBtYXkgc29tZWRheSANCj4gY29t bWl0IGl0LiAgSSBkb24ndCBrbm93IGlmIGFueW9uZSBldmVyIHdvcmtlZCBvbiBjcG93W2ZsXS4N Cj4gSSBzdG9wcGVkIHdvcmtpbmcgb24gcG93bCBhbmQgdGdhbW1hbCB3aGVuIEkgcmV0dXJuZWQg bXkNCj4gY29tbWl0IGJpdCBkdWUgdG8gZGlmZmVyZW5jZXMgd2l0aCAiaGlnaGVyLXJhbmtpbmci IGNvbW1pdHRlcnMuDQo+IA0KDQpiZGUncyBpbXBsZW1lbnRhdGlvbiB3YXMgZmFudGFzdGljLiAg SGUgbWFuYWdlZCB0byB1c2UgdHdvIGRvdWJsZQ0KdmFyaWFibGVzIHRvIHJlcHJlc2VudCBhIHNp bmdsZSBxdWFkcnVwbGUgcHJlY2lzaW9uIHJlYWwgKHVzaW5nIGENCnRlY2huaXF1ZSB3aGljaCBJ IGJlbGlldmUgaXMgc3RhbmRhcmQsIGJ1dCBpdCB3YXMgbmV3IHRvIG1lKS4NCg0KV2h5IGhlIG5l dmVyIGNvbW1pdHRlZCBpdCwgSSBkb24ndCBrbm93LiAgUGVvcGxlIG9uIHRoaXMgbGlzdCB0YWtl IHdheQ0KdG9vIGxvbmcgdG8gY29tbWl0IGV4Y2VsbGVudCBzb2Z0d2FyZS4NCg0KRm9yIGV4YW1w bGUsIEkgaGFkIGNvZGUgZm9yIGNhY29zLCBjYXRhbiwgZXRjLCByZWFkeSBmb3IgYSBsb25nIHRp bWUNCmJlZm9yZSBzb21lb25lIGVsc2Ugc3RlcHBlZCBpbiBhbmQgY29tbWl0dGVkIGl0LiAgQW5k IEkgdGVzdGVkIGl0IHRvDQpkZWF0aC4gIEkgZXZlbiBmb3VuZCBhIGJ1ZyBpbg0KVC4gRS4gSHVs bCwgVGhvbWFzIEYuIEZhaXJncmlldmUsIFBpbmcgVGFrIFBldGVyIFRhbmcsIOKAnEltcGxlbWVu dGluZyB0aGUNCmNvbXBsZXggYXJjc2luIGFuZCBhcmNjb3NpbmUgZnVuY3Rpb25zIHVzaW5nIGV4 Y2VwdGlvbiBoYW5kbGluZ+KAnSwgQUNNDQpUcmFuc2FjdGlvbnMgb24gTWF0aGVtYXRpY2FsIFNv ZnR3YXJlLCBWb2x1bWUgMjMgKDE5OTcpIHBwIDI5OeKAkzMzNS4NCg== From owner-freebsd-numerics@freebsd.org Wed Nov 8 01:50:54 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 52DA3E6A9A4 for ; Wed, 8 Nov 2017 01:50:54 +0000 (UTC) (envelope-from stephen@missouri.edu) Received: from um-nip4-missouri-out.um.umsystem.edu (um-nip4-missouri-out.um.umsystem.edu [198.209.49.177]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "um-tip1.um.umsystem.edu", Issuer "InCommon RSA Server CA" (verified OK)) by mx1.freebsd.org (Postfix) with ESMTPS id E35E080389 for ; Wed, 8 Nov 2017 01:50:53 +0000 (UTC) (envelope-from stephen@missouri.edu) X-IronPort-Anti-Spam-Filtered: true X-IronPort-Anti-Spam-Result: =?us-ascii?q?A2EBAgA3YgJa/yI40cZdGgEBAQEBAgEBA?= =?us-ascii?q?QEIAQEBAYQYbicHg3aLE5AumFoBCSOFGAIahSAUAQEBAQEBAQEBA2gdC4JqR1g?= =?us-ascii?q?BAQEBAQEjAg1eAQUjEVUCAQgYAgImAgICMBUQAgQNCAEBih8QqGmCJ4gXAQeCe?= =?us-ascii?q?gEBCAEBAQEfBYEPgiGFbYMBg0KBUxeCfoJjBaIWAodkjyyGBIsfjGeJDAICAgI?= =?us-ascii?q?JAhqBOTYigXF8E4MuhF53AYsCAYEQAQEB?= X-IPAS-Result: =?us-ascii?q?A2EBAgA3YgJa/yI40cZdGgEBAQEBAgEBAQEIAQEBAYQYbic?= =?us-ascii?q?Hg3aLE5AumFoBCSOFGAIahSAUAQEBAQEBAQEBA2gdC4JqR1gBAQEBAQEjAg1eA?= =?us-ascii?q?QUjEVUCAQgYAgImAgICMBUQAgQNCAEBih8QqGmCJ4gXAQeCegEBCAEBAQEfBYE?= =?us-ascii?q?PgiGFbYMBg0KBUxeCfoJjBaIWAodkjyyGBIsfjGeJDAICAgIJAhqBOTYigXF8E?= =?us-ascii?q?4MuhF53AYsCAYEQAQEB?= Received: from ex-n22.um.umsystem.edu (HELO EX2-N22.um.umsystem.edu) ([198.209.56.34]) by um-nip4-exch-relay.um.umsystem.edu with ESMTP; 07 Nov 2017 19:50:52 -0600 Received: from EX2-N14.um.umsystem.edu (198.209.56.22) by EX2-N22.um.umsystem.edu (198.209.56.34) with Microsoft SMTP Server (version=TLS1_2, cipher=TLS_ECDHE_RSA_WITH_AES_128_CBC_SHA256_P256) id 15.1.1261.35; Tue, 7 Nov 2017 19:50:52 -0600 Received: from EX2-N14.um.umsystem.edu ([198.209.56.22]) by EX2-N14.um.umsystem.edu ([198.209.56.22]) with mapi id 15.01.1261.035; Tue, 7 Nov 2017 19:50:52 -0600 From: "Montgomery-Smith, Stephen" To: "freebsd-numerics@freebsd.org" Subject: Re: cpow and clog Thread-Topic: cpow and clog Thread-Index: AQHTVzhwY4EyR5x3Wk2DmT17BVFyU6MINVmAgAAOeACAAdpVAA== Date: Wed, 8 Nov 2017 01:50:52 +0000 Message-ID: References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171106213308.GA73787@troutmask.apl.washington.edu> In-Reply-To: <20171106213308.GA73787@troutmask.apl.washington.edu> Accept-Language: en-US Content-Language: en-US X-MS-Has-Attach: X-MS-TNEF-Correlator: user-agent: Mozilla/5.0 (X11; Linux x86_64; rv:52.0) Gecko/20100101 Thunderbird/52.4.0 x-originating-ip: [128.206.49.160] Content-Type: text/plain; charset="utf-8" Content-ID: <75E7808CF8357540B6EA1AF6F339ECA8@missouri.edu> Content-Transfer-Encoding: base64 MIME-Version: 1.0 X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 01:50:54 -0000 T24gMTEvMDYvMjAxNyAwMzozMyBQTSwgU3RldmUgS2FyZ2wgd3JvdGU6DQo+IE9uIE1vbiwgTm92 IDA2LCAyMDE3IGF0IDEyOjQxOjIxUE0gLTA4MDAsIFN0ZXZlIEthcmdsIHdyb3RlOg0KPj4gT24g TW9uLCBOb3YgMDYsIDIwMTcgYXQgMDg6NDk6NDNQTSArMDEwMCwgTWljaGFlbCBEYW5pbG92IHdy b3RlOg0KPj4+IEhlbGxvLA0KPj4+DQo+Pj4gSSB3b3VsZCBsaWtlIHRvIGhhdmUgc29tZSBmZWVk YmFjayBvbiBteSBhdHRlbXB0IHRvIGltcG9ydCBPcGVuQlNEDQo+Pj4gY29kZSBmb3IgY3BvdyBh bmQgY2xvZzoNCj4+Pg0KPj4+IGh0dHBzOi8vYnVncy5mcmVlYnNkLm9yZy9idWd6aWxsYS9zaG93 X2J1Zy5jZ2k/aWQ9MjIxMzQxDQo+Pj4gaHR0cHM6Ly9idWdzLmZyZWVic2Qub3JnL2J1Z3ppbGxh L2F0dGFjaG1lbnQuY2dpP2lkPTE4NzY5Mw0KPj4+DQo+Pj4gV2hhdCBoYXBwZW5lZCB0byB0aGUg YWx0ZXJuYXRpdmUgaW1wbGVtZW50YXRpb24gbWVudGlvbmVkIGluIHRoZSB0aHJlYWQgYmVsb3c/ DQo+Pg0KPj4gYmRlIGhhcyBhbiBpbXBsZW1lbnRhdGlvbiBvZiBjbG9nW2ZsXS4gIEhlIG1heSBz b21lZGF5IA0KPj4gY29tbWl0IGl0Lg0KPiANCj4gTm90ZSwgSSBwYWNrYWdlZCBCcnVjZSdzIGNv ZGUgd2l0aCBtYW5wYWdlcyBoZXJlDQo+IA0KPiBodHRwczovL2J1Z3MuZnJlZWJzZC5vcmcvYnVn emlsbGEvc2hvd19idWcuY2dpP2lkPTIxNjg2Mw0KDQpQbGVhc2UsIHNvbWVvbmUgc2hvdWxkIGNv bW1pdCBpdC4NCg== From owner-freebsd-numerics@freebsd.org Wed Nov 8 05:18:43 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 4F4FEE6E269 for ; Wed, 8 Nov 2017 05:18:43 +0000 (UTC) (envelope-from lists@eitanadler.com) Received: from mail-yw0-x22c.google.com (mail-yw0-x22c.google.com [IPv6:2607:f8b0:4002:c05::22c]) (using TLSv1.2 with cipher ECDHE-RSA-AES128-GCM-SHA256 (128/128 bits)) (Client CN "smtp.gmail.com", Issuer "Google Internet Authority G2" (verified OK)) by mx1.freebsd.org (Postfix) with ESMTPS id 16D186A8F9 for ; Wed, 8 Nov 2017 05:18:42 +0000 (UTC) (envelope-from lists@eitanadler.com) Received: by mail-yw0-x22c.google.com with SMTP id w5so1332372ywg.11 for ; Tue, 07 Nov 2017 21:18:42 -0800 (PST) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=eitanadler.com; s=0xdeadbeef; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :cc; bh=ddbMLIz+q6z8Yt6sVsAg0B+YketT4/BZpMsA3lmzqHA=; b=qcphM2rW7FID+vsbBVOgvHt2Ht3Yo+I4+xL/ehE46qbiu91dAr9CcCmCv5GQ6NOCEu eSVVGECQdNZj7Pe1WGyVSKQPm7d0dNqUGIpXDfnL4M0SEaESsGeitENtPxCaI2Lw24s2 lYBhbUa7uIE9FCi69L3pIPd6vXz4jfQf8nXWs= X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20161025; h=x-gm-message-state:mime-version:in-reply-to:references:from:date :message-id:subject:to:cc; bh=ddbMLIz+q6z8Yt6sVsAg0B+YketT4/BZpMsA3lmzqHA=; b=rbLNnc1Qv4y+ucAiHUlL3r5wGAXhMTdnFqzHrrZX9fPG21rnb0oPMjLRfgiTwaExWb BQFw8s24kDgIHOg63FUGTIrLvwFnNZnyG8yvujyR7uYYDSvWkgH2pDDDnB1v5hSyf6La 9KwnWUIP7mw5EbWm9k41zbcAppNT41qw5tJPh1UFFHau7a5r7HLO662rq0f877Vqgcnh QPOGh3eVV71r9BNCew/EBMLpF/1hpF5d/NRGXTo1WmBhInOTJi5hY3Nx3vDEX91o7Zk8 dz17YYGdenWnq1Yxps037Q2fomlcJTd1wwWt6EqWNsU0Hr/bY4TuNVflpSjs1+8j6bgu 5CIw== X-Gm-Message-State: AJaThX4/F92OU7zyJCHwGJ39SHJ/jS/wseUf4CmNLyYH61wy2Pw1xuLz BREqi4cBrzBQSTrVitbmfb2ZlsZFP5Ikv/u3nvGGmA== X-Google-Smtp-Source: ABhQp+RXSx7LZrr7E+WYWGsrLPJnOBhDL8BArPpvCUs5/KrxlSeP9l9bOvFNM20FEHMVblDpsuPoiayrZKKHw7z92ko= X-Received: by 10.37.199.137 with SMTP id w131mr661376ybe.118.1510118322018; Tue, 07 Nov 2017 21:18:42 -0800 (PST) MIME-Version: 1.0 Received: by 10.37.56.66 with HTTP; Tue, 7 Nov 2017 21:18:11 -0800 (PST) In-Reply-To: References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171106213308.GA73787@troutmask.apl.washington.edu> From: Eitan Adler Date: Tue, 7 Nov 2017 21:18:11 -0800 Message-ID: Subject: Re: cpow and clog To: "Montgomery-Smith, Stephen" Cc: "freebsd-numerics@freebsd.org" Content-Type: text/plain; charset="UTF-8" X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 05:18:43 -0000 On 7 November 2017 at 17:50, Montgomery-Smith, Stephen wrote: >> Note, I packaged Bruce's code with manpages here >> >> https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=216863 > > Please, someone should commit it. :( patching file include/complex.h Hunk #1 succeeded at 99 (offset 12 lines). patching file lib/libc/softfloat/bits64/softfloat-macros patching file lib/msun/Makefile Hunk #2 succeeded at 134 (offset 3 lines). Hunk #3 succeeded at 171 (offset 7 lines). patching file lib/msun/Symbol.map Hunk #1 FAILED at 294. 1 out of 1 hunk FAILED -- saving rejects to file lib/msun/Symbol.map.rej patching file lib/msun/man/clog.3 patching file lib/msun/man/complex.3 Hunk #1 succeeded at 24 with fuzz 2. patching file lib/msun/src/math_private.h Hunk #2 FAILED at 309. 1 out of 2 hunks FAILED -- saving rejects to file lib/msun/src/math_private.h.rej patching file lib/msun/src/s_clog.c patching file lib/msun/src/s_clogf.c patching file lib/msun/src/s_clogl.c -- Eitan Adler From owner-freebsd-numerics@freebsd.org Wed Nov 8 06:15:05 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 05D7BE6F012 for ; Wed, 8 Nov 2017 06:15:05 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "troutmask", Issuer "troutmask" (not verified)) by mx1.freebsd.org (Postfix) with ESMTPS id D73AE6C72C for ; Wed, 8 Nov 2017 06:15:04 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (localhost [127.0.0.1]) by troutmask.apl.washington.edu (8.15.2/8.15.2) with ESMTPS id vA86F2Gv051071 (version=TLSv1.2 cipher=DHE-RSA-AES256-GCM-SHA384 bits=256 verify=NO); Tue, 7 Nov 2017 22:15:02 -0800 (PST) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.15.2/8.15.2/Submit) id vA86F1wp051070; Tue, 7 Nov 2017 22:15:01 -0800 (PST) (envelope-from sgk) Date: Tue, 7 Nov 2017 22:15:01 -0800 From: Steve Kargl To: Eitan Adler Cc: "Montgomery-Smith, Stephen" , "freebsd-numerics@freebsd.org" Subject: Re: cpow and clog Message-ID: <20171108061501.GA51059@troutmask.apl.washington.edu> Reply-To: sgk@troutmask.apl.washington.edu References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171106213308.GA73787@troutmask.apl.washington.edu> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: User-Agent: Mutt/1.7.2 (2016-11-26) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 06:15:05 -0000 On Tue, Nov 07, 2017 at 09:18:11PM -0800, Eitan Adler wrote: > On 7 November 2017 at 17:50, Montgomery-Smith, Stephen > wrote: > > >> Note, I packaged Bruce's code with manpages here > >> > >> https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=216863 > > > > Please, someone should commit it. > > :( Well, yes, a patch that sits in bugzilla tends to rot. > patching file include/complex.h > Hunk #1 succeeded at 99 (offset 12 lines). > patching file lib/libc/softfloat/bits64/softfloat-macros > patching file lib/msun/Makefile > Hunk #2 succeeded at 134 (offset 3 lines). > Hunk #3 succeeded at 171 (offset 7 lines). > patching file lib/msun/Symbol.map > Hunk #1 FAILED at 294. Probably a trivial edit to deal with clog[fl] placed in different order of Symbol.map. > 1 out of 1 hunk FAILED -- saving rejects to file lib/msun/Symbol.map.rej > patching file lib/msun/man/clog.3 > patching file lib/msun/man/complex.3 > Hunk #1 succeeded at 24 with fuzz 2. > patching file lib/msun/src/math_private.h > Hunk #2 FAILED at 309. > 1 out of 2 hunks FAILED -- saving rejects to file > lib/msun/src/math_private.h.rej Don't know if someone has made change sot math_private.h. I guess I'll have to rebase my patch -- Steve From owner-freebsd-numerics@freebsd.org Wed Nov 8 11:00:45 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 9316EE4EF75 for ; Wed, 8 Nov 2017 11:00:45 +0000 (UTC) (envelope-from dimpase@gmail.com) Received: from mail-lf0-x22a.google.com (mail-lf0-x22a.google.com [IPv6:2a00:1450:4010:c07::22a]) (using TLSv1.2 with cipher ECDHE-RSA-AES128-GCM-SHA256 (128/128 bits)) (Client CN "smtp.gmail.com", Issuer "Google Internet Authority G2" (verified OK)) by mx1.freebsd.org (Postfix) with ESMTPS id 1259374CC6 for ; Wed, 8 Nov 2017 11:00:45 +0000 (UTC) (envelope-from dimpase@gmail.com) Received: by mail-lf0-x22a.google.com with SMTP id a16so2640657lfk.0 for ; Wed, 08 Nov 2017 03:00:44 -0800 (PST) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20161025; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :cc; bh=WvHD1gYrGifVQp3Y87Pj5K0qOGPpFI29CPv4W+CgW0c=; b=eFeA8Lgdx1rHwa2rLCuqp6/DUNBnoxeUwceUSu1C7x0EVl6N/58fNYW1xkia/grfTO k2uZzqCsg1SMkGu+0cl6aAVJuCgg9EqzwIq+E3Zc7r6kxpiFGgeslwHl1iRdON8eQTIY 9DdGMFMFlNFFW4gbsEb2w9iaCZfq8q9hkdD8aRCztpuAZeI2CNAld1ZfD3IT7ZcBqnXj o5BLk7c2t+urx+3Ic7ySNEDnmUm+jA5RhZJwXA0gviKRIOVyCUsxIamFo48Iq5lJmpXj uocqNqjdECOgbjj1fGbWdlE6eVfN13L5F1Nwd6C52hLltpzeeyPTA3FY5kzlFD9xN0Da aEKA== X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20161025; h=x-gm-message-state:mime-version:in-reply-to:references:from:date :message-id:subject:to:cc; bh=WvHD1gYrGifVQp3Y87Pj5K0qOGPpFI29CPv4W+CgW0c=; b=Nua86Yh9m1kw2Ozp8bdX4nOLpCS9i7ntqNOqn37K3vqj68zHJLwtg9CpiVZpbL8se0 dZ4dXnxoDCjPw3uBVHCUp5k7tMA76npJKQ+goIKdRcGUpWfuSTsChLkPoYUOXf2So509 OYgIp+RluronEsV4hGUWIYHGSUAQi+KYJyPNi6GkXHW4AtOzOSCCGdE3t+Xa9CFqHo70 G6D4W9egnxvKrLoiIuy5YEAWBCiGA3S8WedSN6tOdUsIj10yEXTPyQle3MDhEo3sduth +cNnRyAvyU1yxSYtY+W7ygDdY0mbHraO/mhTIT4EV7fVhBbV562rTQ5MYe9wvoMp7qbo raBQ== X-Gm-Message-State: AJaThX5lnp9RnMzUVKzkLhBsBI9VydKQCBVA1Sz7sBXG08q21duTj3l+ SZ7WxXYmYwLQYtfpLDq9WCzZE8liDtcXEojx8Ws/LA== X-Google-Smtp-Source: ABhQp+RCotYF2MjM3+U82JIRUK0e4xQG8UNFh5GYKeH3SqLG7UD/A2WK8K7Rj3BdXYoiAntGBcHZ0XtD9CVxk6kRDo0= X-Received: by 10.46.83.67 with SMTP id t3mr38089ljd.135.1510138840952; Wed, 08 Nov 2017 03:00:40 -0800 (PST) MIME-Version: 1.0 Received: by 10.46.29.68 with HTTP; Wed, 8 Nov 2017 03:00:40 -0800 (PST) Received: by 10.46.29.68 with HTTP; Wed, 8 Nov 2017 03:00:40 -0800 (PST) In-Reply-To: <20171108061501.GA51059@troutmask.apl.washington.edu> References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171106213308.GA73787@troutmask.apl.washington.edu> <20171108061501.GA51059@troutmask.apl.washington.edu> From: Dima Pasechnik Date: Wed, 8 Nov 2017 11:00:40 +0000 Message-ID: Subject: Re: cpow and clog To: Steve Kargl Cc: Eitan Adler , "Montgomery-Smith, Stephen" , "freebsd-numerics@freebsd.org" Content-Type: text/plain; charset="UTF-8" X-Content-Filtered-By: Mailman/MimeDel 2.1.23 X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 11:00:45 -0000 On 8 Nov 2017 07:15, "Steve Kargl" wrote: On Tue, Nov 07, 2017 at 09:18:11PM -0800, Eitan Adler wrote: > On 7 November 2017 at 17:50, Montgomery-Smith, Stephen > wrote: > > >> Note, I packaged Bruce's code with manpages here > >> > >> https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=216863 > > > > Please, someone should commit it. > > :( Well, yes, a patch that sits in bugzilla tends to rot. A patch not tagged by a version control system surely rots very quickly, as noone except the submitter knows the exact state of the source tree it is made against, especially if the tree is also not under a revision control system. And surely the submitter's memory is not perfect, either. If your patches were made on a version control system, rebasing them would have been mist probably trivial, too (assuming it is just some permutation of commuting chunks of code causing patch to fail, as it most often the case). That is, one reason for these patches to move so slowly is simply very outdated workflow. Dima PS. The current freebsd workflow, using svn instead of git (as basically everyone else does nowadays) is also oudated, but by a much smaller margin. Still, at least if used properly, patches are perfectly traceable to the state of the source tree they are created at. PPS. yes, in case you may wonder, I have learnt to program using punch cards, almost 40 years ago :-) > patching file include/complex.h > Hunk #1 succeeded at 99 (offset 12 lines). > patching file lib/libc/softfloat/bits64/softfloat-macros > patching file lib/msun/Makefile > Hunk #2 succeeded at 134 (offset 3 lines). > Hunk #3 succeeded at 171 (offset 7 lines). > patching file lib/msun/Symbol.map > Hunk #1 FAILED at 294. Probably a trivial edit to deal with clog[fl] placed in different order of Symbol.map. > 1 out of 1 hunk FAILED -- saving rejects to file lib/msun/Symbol.map.rej > patching file lib/msun/man/clog.3 > patching file lib/msun/man/complex.3 > Hunk #1 succeeded at 24 with fuzz 2. > patching file lib/msun/src/math_private.h > Hunk #2 FAILED at 309. > 1 out of 2 hunks FAILED -- saving rejects to file > lib/msun/src/math_private.h.rej Don't know if someone has made change sot math_private.h. I guess I'll have to rebase my patch -- Steve _______________________________________________ freebsd-numerics@freebsd.org mailing list https://lists.freebsd.org/mailman/listinfo/freebsd-numerics To unsubscribe, send any mail to "freebsd-numerics-unsubscribe@freebsd.org" From owner-freebsd-numerics@freebsd.org Wed Nov 8 11:12:54 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 18FF0E4F4D0 for ; Wed, 8 Nov 2017 11:12:54 +0000 (UTC) (envelope-from bugzilla-noreply@freebsd.org) Received: from kenobi.freebsd.org (kenobi.freebsd.org [IPv6:2001:1900:2254:206a::16:76]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client did not present a certificate) by mx1.freebsd.org (Postfix) with ESMTPS id 0762E75480 for ; Wed, 8 Nov 2017 11:12:54 +0000 (UTC) (envelope-from bugzilla-noreply@freebsd.org) Received: from bugs.freebsd.org ([127.0.1.118]) by kenobi.freebsd.org (8.15.2/8.15.2) with ESMTP id vA8BCqOe004519 for ; Wed, 8 Nov 2017 11:12:53 GMT (envelope-from bugzilla-noreply@freebsd.org) From: bugzilla-noreply@freebsd.org To: freebsd-numerics@FreeBSD.org Subject: [Bug 222247] libmsun tests does not compile during a buildworld Date: Wed, 08 Nov 2017 11:12:53 +0000 X-Bugzilla-Reason: CC X-Bugzilla-Type: changed X-Bugzilla-Watch-Reason: None X-Bugzilla-Product: Base System X-Bugzilla-Component: bin X-Bugzilla-Version: 11.1-STABLE X-Bugzilla-Keywords: X-Bugzilla-Severity: Affects Only Me X-Bugzilla-Who: bsdpr@phoe.frmug.org X-Bugzilla-Status: Closed X-Bugzilla-Resolution: FIXED X-Bugzilla-Priority: --- X-Bugzilla-Assigned-To: rlibby@freebsd.org X-Bugzilla-Flags: X-Bugzilla-Changed-Fields: resolution bug_status Message-ID: In-Reply-To: References: Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: quoted-printable X-Bugzilla-URL: https://bugs.freebsd.org/bugzilla/ Auto-Submitted: auto-generated MIME-Version: 1.0 X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 11:12:54 -0000 https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=3D222247 Bertrand Petit changed: What |Removed |Added ---------------------------------------------------------------------------- Resolution|--- |FIXED Status|In Progress |Closed --=20 You are receiving this mail because: You are on the CC list for the bug.= From owner-freebsd-numerics@freebsd.org Wed Nov 8 14:33:08 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id C06F9E5518D for ; Wed, 8 Nov 2017 14:33:08 +0000 (UTC) (envelope-from brde@optusnet.com.au) Received: from mail107.syd.optusnet.com.au (mail107.syd.optusnet.com.au [211.29.132.53]) by mx1.freebsd.org (Postfix) with ESMTP id 6F2F27D0C8 for ; Wed, 8 Nov 2017 14:33:08 +0000 (UTC) (envelope-from brde@optusnet.com.au) Received: from [192.168.0.104] (c110-21-101-228.carlnfd1.nsw.optusnet.com.au [110.21.101.228]) by mail107.syd.optusnet.com.au (Postfix) with ESMTPS id DDF00D463B5; Thu, 9 Nov 2017 01:33:00 +1100 (AEDT) Date: Thu, 9 Nov 2017 01:33:00 +1100 (EST) From: Bruce Evans X-X-Sender: bde@besplex.bde.org To: "Montgomery-Smith, Stephen" cc: "freebsd-numerics@freebsd.org" Subject: Re: cpow and clog In-Reply-To: <4330a817-ae31-6336-7aa6-3e9981b02ab7@missouri.edu> Message-ID: <20171108224608.B1089@besplex.bde.org> References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <4330a817-ae31-6336-7aa6-3e9981b02ab7@missouri.edu> MIME-Version: 1.0 X-Optus-CM-Score: 0 X-Optus-CM-Analysis: v=2.2 cv=bc8baKHB c=1 sm=1 tr=0 a=PalzARQSbocsUSjMRkwAPg==:117 a=PalzARQSbocsUSjMRkwAPg==:17 a=nlC_4_pT8q9DhB4Ho9EA:9 a=6I5d2MoRAAAA:8 a=zEAycw12O3tiVqPjVFMA:9 a=45ClL6m2LaAA:10 a=IjZwj45LgO3ly-622nXo:22 Content-Type: TEXT/PLAIN; charset=X-UNKNOWN; format=flowed Content-Transfer-Encoding: QUOTED-PRINTABLE X-Content-Filtered-By: Mailman/MimeDel 2.1.23 X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 14:33:08 -0000 On Wed, 8 Nov 2017, Montgomery-Smith, Stephen wrote: > On 11/06/2017 02:41 PM, Steve Kargl wrote: >> On Mon, Nov 06, 2017 at 08:49:43PM +0100, Michael Danilov wrote: >>> Hello, >>> >>> I would like to have some feedback on my attempt to import OpenBSD >>> code for cpow and clog: >>> >>> https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=3D221341 >>> https://bugs.freebsd.org/bugzilla/attachment.cgi?id=3D187693 > > The problem with this implementation of clog is that they will produce > very inaccurate results when |z| is close to 1. A proper implementation > is given in the paper: > T. E. Hull, Thomas F. Fairgrieve, Ping Tak Peter Tang, =E2=80=9CImplement= ing > Complex Elementary Functions Using Exception Handling=E2=80=9D, ACM Trans= actions > on Mathematical Software, Volume 20 (1994), pp 215=E2=80=93244, Corrigend= a, p553. >=20 > They use a precision which is twice that of the required precision. For > example, if you want log(sqrt(x^2+y^2)) accurate to about 17 digits, you > will need to compute it to 34 digits. I don't see how that can be enough. I found that about 3 times as much precision is needed for the worst case in double precision. Closer to 160 bits than 150, and 2 times as much much is only 106 bits. This defeated my intuition, but it was easy to find the worst case by searching a few billion cases near z =3D 1 to improve my intuition, then a few more in a smaller search area to find the worst case. The authors of the paper were/are real numerical analysts, but CPUs were too slow in 1004 for them to routinely search billions of cases to check their algorithms. > If all you care about is getting a a good relative error for the total > answer (rather than insisting on a good relative error for the real and > imaginary parts separately), then for |x+I*y| close to one, compute the > real part using > > 0.5*logp1((x-1)*(x+1)+y*y) z =3D 1 is a very non-exceptional value, so clog() should be very accurate near it. I just noticed that draft C99 (n869) Annex F (IEEE754 binding) is bad enough to say that cpow(z, c) is implicitly specified by the simplistic formula: cpow(z, c) =3D cexp(c * clog(z)) (this is in a section with similar specifications like: ccos(z) =3D ccosh(I * z)) which are not incorrect since they don't lose precision or give spurious exceptions (provided I * z is not evaluated simplistically)). This bug is missing without Annex F, by allowing any quality of implementation. This bug is fixed in C99 through draft C11 (n1570) Annex by not requiring anything extra except raising all relevant exception and saying in a footnote that the above low quality implementation conforms since it raises all relevant exceptions (it might raise spurious exceptions but this conforms). Anyway, if you use the simplistic formula for cpow(), then clog() should be very accurate to avoid exponentiating a large error. Even the best possible minimal error of 0.5 ulps is not good. >>> What happened to the alternative implementation mentioned in the thread= below? >> >> bde has an implementation of clog[fl]. He may someday >> commit it. I don't know if anyone ever worked on cpow[fl]. >> I stopped working on powl and tgammal when I returned my >> commit bit due to differences with "higher-ranking" committers. > > bde's implementation was fantastic. He managed to use two double > variables to represent a single quadruple precision real (using a quad-double > technique which I believe is standard, but it was new to me). Thanks. I forgot the details, but just checked some. It uses more like 4 double variables to represent a single tri-double precision variable. Full quad-double precision would require something like 3 ripple carries per addition and I managed to reduce this something like only 1 ripple carry per addition by doing bi-double additions in a delicate order. This loses some precision but retains enough. This is least delicate near z =3D 1, but it is more accurate when it works, and not much slower, so it is used for approx. 0.5 < |z|^2 < 3. I don't know if this is standard. Other parts of the algorithm are standar= d. > Why he never committed it, I don't know. People on this list take way > too long to commit excellent software. > > For example, I had code for cacos, catan, etc, ready for a long time > before someone else stepped in and committed it. And I tested it to > death. I even found a bug in > T. E. Hull, Thomas F. Fairgrieve, Ping Tak Peter Tang, =E2=80=9CImplement= ing the > complex arcsin and arccosine functions using exception handling=E2=80=9D,= ACM > Transactions on Mathematical Software, Volume 23 (1997) pp 299=E2=80=9333= 5. We were waiting for you to commit your code :-). It is still not complete. It is missing many engineering details: - file names not in a normal namespace (ca* instead of normal e_c* or s_c*)= =2E East to fix, especially before committing. - multiple functions per file. Not so easy to fix. I didn't manage to fix it for my log*l(). I tried many not so good splittings using #include's of a .c file with ifdefs, and different kernel functions, and inline functions. Inline functions are now optimized well enough by clang (but static ones waste space and extern ones are hard to manage). My clog*() file has the same problem. Steve split it, but I haven't got around to even testing that. - switching the rounding precision for long double functions on i386. I couldn't see a clean way to do this. It can be done by sprinkling ENTERIT() and RETURNI() macros. ENTERIT() (ENTERI() with the return type specified) is also needed for clogl(), but hasn't been committed yet. - comment stripping and/or generating source files from a template. You used this for development. The scripts to do this were not very suitabl= e for committing, but perhaps better than nothing. IIRC, the template was the double version with a few extras. The double version has full comments, and all comments are stripped in the float and double versions= , including in the committed versions. But the normal style is only strip large comments and adjust small comments in the float and double versions. My clog*() handles this problem in the opposite way not strip= ping any comments and minimizing differences in comments. It also minimizes difference in code by manual inspection/rewrite instead of a template. Bruce From owner-freebsd-numerics@freebsd.org Wed Nov 8 15:03:05 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 596F7E55DA8 for ; Wed, 8 Nov 2017 15:03:05 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "troutmask", Issuer "troutmask" (not verified)) by mx1.freebsd.org (Postfix) with ESMTPS id 396C37E227 for ; Wed, 8 Nov 2017 15:03:05 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (localhost [127.0.0.1]) by troutmask.apl.washington.edu (8.15.2/8.15.2) with ESMTPS id vA8F337b053970 (version=TLSv1.2 cipher=DHE-RSA-AES256-GCM-SHA384 bits=256 verify=NO); Wed, 8 Nov 2017 07:03:03 -0800 (PST) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.15.2/8.15.2/Submit) id vA8F325h053969; Wed, 8 Nov 2017 07:03:02 -0800 (PST) (envelope-from sgk) Date: Wed, 8 Nov 2017 07:03:02 -0800 From: Steve Kargl To: Dima Pasechnik Cc: Eitan Adler , "Montgomery-Smith, Stephen" , "freebsd-numerics@freebsd.org" Subject: Re: cpow and clog Message-ID: <20171108150302.GA53567@troutmask.apl.washington.edu> Reply-To: sgk@troutmask.apl.washington.edu References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171106213308.GA73787@troutmask.apl.washington.edu> <20171108061501.GA51059@troutmask.apl.washington.edu> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: User-Agent: Mutt/1.7.2 (2016-11-26) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 15:03:05 -0000 On Wed, Nov 08, 2017 at 11:00:40AM +0000, Dima Pasechnik wrote: > > A patch not tagged by a version control system surely rots very quickly, > as noone except the submitter knows the exact state of the source tree it > is made against, especially if the tree is also not under a revision > control system. > And surely the submitter's memory is not perfect, either. The patch is against trunk on the day it was submitted. The patch was created under the officical revision control system used by FreeBSD. -- Steve From owner-freebsd-numerics@freebsd.org Wed Nov 8 15:58:28 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id E115CE56F52 for ; Wed, 8 Nov 2017 15:58:28 +0000 (UTC) (envelope-from lwhsu.freebsd@gmail.com) Received: from mail-wm0-f44.google.com (mail-wm0-f44.google.com [74.125.82.44]) (using TLSv1.2 with cipher ECDHE-RSA-AES128-GCM-SHA256 (128/128 bits)) (Client CN "smtp.gmail.com", Issuer "Google Internet Authority G2" (verified OK)) by mx1.freebsd.org (Postfix) with ESMTPS id A1AA37FB13 for ; Wed, 8 Nov 2017 15:58:28 +0000 (UTC) (envelope-from lwhsu.freebsd@gmail.com) Received: by mail-wm0-f44.google.com with SMTP id t139so11856897wmt.1 for ; Wed, 08 Nov 2017 07:58:28 -0800 (PST) X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20161025; h=x-gm-message-state:mime-version:references:in-reply-to:from:date :message-id:subject:to:cc; bh=DqaYIpAFrMyxhe4mhO3Ya+2wy+D19IifIG3s8Tk3YfY=; b=mvLqdTl2dg1BS+EJkNRGxl8wEvQyHuuui5Abm7lkJS6sDHITooML4dCEGVsY8kZwPJ 4H96zm0CbpDiAIHAK90DzNxlZ7MKkiHN6xKgzimhEbm94vvWJZswvPq8Q9E5Aknlif18 NeE8Cn6w1EMoSj0U3y6U3o82mCHm0RbXTS88iUtv+da49aKWp1T08mp1EXAmu7oCTVLb Gs+rFL/3lKL9h5+y+Y0B2im823u61O9nXTRY4ElFPyp2y8wRq4tKWqi3H75j4z9jIYTu JGBL+u0nGkPu0QRKLw6ek4oznMBfrx3hd8oHhykjMrGgPFhnowYhDFiGFlv8rKwwzoox zyPg== X-Gm-Message-State: AJaThX7bHQFD4vC6ewds10EKuB2+6Reu3+OSvcTlL7GpDfwMKopwNY68 7v5n1hhJMKPQ311jvXpmT3qJ4IZRFAvr2bjPWqs= X-Google-Smtp-Source: ABhQp+S5D4oK1lJoheIw5+Eebb4P+Q/6pUmt3aOovjDpw4K3qTZqm8le2+Wp+wARMfD729B0qpgrDJ07VNUxlNn2nVc= X-Received: by 10.28.35.9 with SMTP id j9mr707559wmj.109.1510156700679; Wed, 08 Nov 2017 07:58:20 -0800 (PST) MIME-Version: 1.0 References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171106213308.GA73787@troutmask.apl.washington.edu> <20171108061501.GA51059@troutmask.apl.washington.edu> <20171108150302.GA53567@troutmask.apl.washington.edu> In-Reply-To: <20171108150302.GA53567@troutmask.apl.washington.edu> From: Li-Wen Hsu Date: Wed, 08 Nov 2017 15:58:10 +0000 Message-ID: Subject: Re: cpow and clog To: Dima Pasechnik , sgk@troutmask.apl.washington.edu Cc: Eitan Adler , "Montgomery-Smith, Stephen" , "freebsd-numerics@freebsd.org" Content-Type: text/plain; charset="UTF-8" Content-Transfer-Encoding: quoted-printable X-Content-Filtered-By: Mailman/MimeDel 2.1.23 X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 15:58:29 -0000 On 2017=E5=B9=B411=E6=9C=888=E6=97=A5 =E9=80=B1=E4=B8=89 at 23:03 Steve Kar= gl wrote: > On Wed, Nov 08, 2017 at 11:00:40AM +0000, Dima Pasechnik wrote: > > > > A patch not tagged by a version control system surely rots very quickl= y, > > as noone except the submitter knows the exact state of the source tree = it > > is made against, especially if the tree is also not under a revision > > control system. > > And surely the submitter's memory is not perfect, either. > > The patch is against trunk on the day it was submitted. > The patch was created under the officical revision control > system used by FreeBSD. We also have an official patch review system: https://reviews.freebsd.org , which is not perfect, but still useable. This might help Dima's concern. Li-Wen --=20 Li-Wen Hsu https://lwhsu.org From owner-freebsd-numerics@freebsd.org Wed Nov 8 16:30:55 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id 0CBABE578D2 for ; Wed, 8 Nov 2017 16:30:55 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "troutmask", Issuer "troutmask" (not verified)) by mx1.freebsd.org (Postfix) with ESMTPS id DEAD280B48; Wed, 8 Nov 2017 16:30:54 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (localhost [127.0.0.1]) by troutmask.apl.washington.edu (8.15.2/8.15.2) with ESMTPS id vA8GUqLM055317 (version=TLSv1.2 cipher=DHE-RSA-AES256-GCM-SHA384 bits=256 verify=NO); Wed, 8 Nov 2017 08:30:52 -0800 (PST) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.15.2/8.15.2/Submit) id vA8GUqag055316; Wed, 8 Nov 2017 08:30:52 -0800 (PST) (envelope-from sgk) Date: Wed, 8 Nov 2017 08:30:52 -0800 From: Steve Kargl To: Li-Wen Hsu Cc: Dima Pasechnik , Eitan Adler , "Montgomery-Smith, Stephen" , "freebsd-numerics@freebsd.org" Subject: Re: cpow and clog Message-ID: <20171108163052.GA55224@troutmask.apl.washington.edu> Reply-To: sgk@troutmask.apl.washington.edu References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171106213308.GA73787@troutmask.apl.washington.edu> <20171108061501.GA51059@troutmask.apl.washington.edu> <20171108150302.GA53567@troutmask.apl.washington.edu> MIME-Version: 1.0 Content-Type: text/plain; charset=utf-8 Content-Disposition: inline Content-Transfer-Encoding: 8bit In-Reply-To: User-Agent: Mutt/1.7.2 (2016-11-26) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 16:30:55 -0000 On Wed, Nov 08, 2017 at 03:58:10PM +0000, Li-Wen Hsu wrote: > On 2017年11月8日 週三 at 23:03 Steve Kargl > wrote: > > > On Wed, Nov 08, 2017 at 11:00:40AM +0000, Dima Pasechnik wrote: > > > > > > A patch not tagged by a version control system surely rots very quickly, > > > as noone except the submitter knows the exact state of the source tree it > > > is made against, especially if the tree is also not under a revision > > > control system. > > > And surely the submitter's memory is not perfect, either. > > > > The patch is against trunk on the day it was submitted. > > The patch was created under the officical revision control > > system used by FreeBSD. > > > We also have an official patch review system: https://reviews.freebsd.org , > which is not perfect, but still useable. This might help Dima's concern. > Looks like it requires "yet another account". -- Steve From owner-freebsd-numerics@freebsd.org Wed Nov 8 17:11:23 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id ED53BE5829A for ; Wed, 8 Nov 2017 17:11:23 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "troutmask", Issuer "troutmask" (not verified)) by mx1.freebsd.org (Postfix) with ESMTPS id C1F261D73 for ; Wed, 8 Nov 2017 17:11:23 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (localhost [127.0.0.1]) by troutmask.apl.washington.edu (8.15.2/8.15.2) with ESMTPS id vA8HBMLP056154 (version=TLSv1.2 cipher=DHE-RSA-AES256-GCM-SHA384 bits=256 verify=NO); Wed, 8 Nov 2017 09:11:22 -0800 (PST) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.15.2/8.15.2/Submit) id vA8HBMKD056153; Wed, 8 Nov 2017 09:11:22 -0800 (PST) (envelope-from sgk) Date: Wed, 8 Nov 2017 09:11:22 -0800 From: Steve Kargl To: Dima Pasechnik Cc: "Montgomery-Smith, Stephen" , Eitan Adler , "freebsd-numerics@freebsd.org" Subject: Re: cpow and clog Message-ID: <20171108171122.GA56093@troutmask.apl.washington.edu> Reply-To: sgk@troutmask.apl.washington.edu References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171106213308.GA73787@troutmask.apl.washington.edu> <20171108061501.GA51059@troutmask.apl.washington.edu> <20171108150302.GA53567@troutmask.apl.washington.edu> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: <20171108150302.GA53567@troutmask.apl.washington.edu> User-Agent: Mutt/1.7.2 (2016-11-26) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Wed, 08 Nov 2017 17:11:24 -0000 On Wed, Nov 08, 2017 at 07:03:02AM -0800, Steve Kargl wrote: > On Wed, Nov 08, 2017 at 11:00:40AM +0000, Dima Pasechnik wrote: > > > > A patch not tagged by a version control system surely rots very quickly, > > as noone except the submitter knows the exact state of the source tree it > > is made against, especially if the tree is also not under a revision > > control system. > > And surely the submitter's memory is not perfect, either. > > The patch is against trunk on the day it was submitted. > The patch was created under the officical revision control > system used by FreeBSD. > Hmmm. Just looked at the patch in bugzilla. The top four lines are Index: include/complex.h =================================================================== --- include/complex.h (revision 313777) +++ include/complex.h (working copy) So, the patchs tagged with the revision number. HTH -- Steve From owner-freebsd-numerics@freebsd.org Thu Nov 9 01:23:01 2017 Return-Path: Delivered-To: freebsd-numerics@mailman.ysv.freebsd.org Received: from mx1.freebsd.org (mx1.freebsd.org [IPv6:2001:1900:2254:206a::19:1]) by mailman.ysv.freebsd.org (Postfix) with ESMTP id A1411E6118A for ; Thu, 9 Nov 2017 01:23:01 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (troutmask.apl.washington.edu [128.95.76.21]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "troutmask", Issuer "troutmask" (not verified)) by mx1.freebsd.org (Postfix) with ESMTPS id 7CF0A6FCB0 for ; Thu, 9 Nov 2017 01:23:01 +0000 (UTC) (envelope-from sgk@troutmask.apl.washington.edu) Received: from troutmask.apl.washington.edu (localhost [127.0.0.1]) by troutmask.apl.washington.edu (8.15.2/8.15.2) with ESMTPS id vA91N0E4006428 (version=TLSv1.2 cipher=DHE-RSA-AES256-GCM-SHA384 bits=256 verify=NO); Wed, 8 Nov 2017 17:23:00 -0800 (PST) (envelope-from sgk@troutmask.apl.washington.edu) Received: (from sgk@localhost) by troutmask.apl.washington.edu (8.15.2/8.15.2/Submit) id vA91Mxvg006427; Wed, 8 Nov 2017 17:22:59 -0800 (PST) (envelope-from sgk) Date: Wed, 8 Nov 2017 17:22:59 -0800 From: Steve Kargl To: Dima Pasechnik Cc: "Montgomery-Smith, Stephen" , Eitan Adler , "freebsd-numerics@freebsd.org" Subject: Re: cpow and clog Message-ID: <20171109012259.GA6402@troutmask.apl.washington.edu> Reply-To: sgk@troutmask.apl.washington.edu References: <20171106194937.GA87725@freebird> <20171106204121.GB37361@troutmask.apl.washington.edu> <20171106213308.GA73787@troutmask.apl.washington.edu> <20171108061501.GA51059@troutmask.apl.washington.edu> <20171108150302.GA53567@troutmask.apl.washington.edu> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Disposition: inline In-Reply-To: <20171108150302.GA53567@troutmask.apl.washington.edu> User-Agent: Mutt/1.7.2 (2016-11-26) X-BeenThere: freebsd-numerics@freebsd.org X-Mailman-Version: 2.1.23 Precedence: list List-Id: "Discussions of high quality implementation of libm functions." List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , X-List-Received-Date: Thu, 09 Nov 2017 01:23:01 -0000 On Wed, Nov 08, 2017 at 07:03:02AM -0800, Steve Kargl wrote: > On Wed, Nov 08, 2017 at 11:00:40AM +0000, Dima Pasechnik wrote: > > > > A patch not tagged by a version control system surely rots very quickly, > > as noone except the submitter knows the exact state of the source tree it > > is made against, especially if the tree is also not under a revision > > control system. > > And surely the submitter's memory is not perfect, either. > > The patch is against trunk on the day it was submitted. > The patch was created under the officical revision control > system used by FreeBSD. > I just updated my patch in bugzilla. https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=216863 I have no idea what Dima's problems were with the original patch that was attached to bugzilla as I applied that patch to up-to-date trunk without issue. svn info Path: . Working Copy Root Path: /usr/home/kargl/freebsd/src URL: svn+ssh://kargl@svn.freebsd.org/base/head Relative URL: ^/head Repository Root: svn+ssh://kargl@svn.freebsd.org/base Repository UUID: ccf9f872-aa2e-dd11-9fc8-001c23d0bc1f Revision: 325568 HTH -- Steve 20170425 https://www.youtube.com/watch?v=VWUpyCsUKR4 20161221 https://www.youtube.com/watch?v=IbCHE-hONow