From owner-freebsd-numerics@freebsd.org Wed Feb 8 11:03:40 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 93A90CD63C7 for ; Wed, 8 Feb 2017 11:03:40 +0000 (UTC) (envelope-from mokhi64@gmail.com) Received: from mail-io0-x230.google.com (mail-io0-x230.google.com [IPv6:2607:f8b0:4001:c06::230]) (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 621B91BD6 for ; Wed, 8 Feb 2017 11:03:40 +0000 (UTC) (envelope-from mokhi64@gmail.com) Received: by mail-io0-x230.google.com with SMTP id j13so112817361iod.3 for ; Wed, 08 Feb 2017 03:03:40 -0800 (PST) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20161025; h=mime-version:from:date:message-id:subject:to; bh=5UEQ0H+g5lUNJOlIDx5s+yZDhW/9VmEijNLffdEEsTA=; b=VbkkHFe8GqqwHIhd7k/urxHhp4rxcxo4wBVRqy9u58xxthO/sRXVPVQDq+nKma86KU iMKOq7e+7BVbWlNFAlpzvtMRmheDs6DvPJ2S79mpB8j9QcpgFv8sKM3Zxk1sPo0eYQ/C feV61cpo35s4/Mi0jnD94CxWcWLxAXty4x8CJGc2aMlCbQK4JVN6VSuxmVRIhFeUBmeO chNFOvPdSFtaTdrkWEJddyUkhWwILUh0Us+qN2QB6fIhMN3QP2AuQJ+tKpdWeYgis3xM 3+IfIJAXmeXf5GAQExnsx67MjHOY289rDTgE6nxN8ns9LE2bhP+MsL2c75r7Z1p691nh i3gA== X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20161025; h=x-gm-message-state:mime-version:from:date:message-id:subject:to; bh=5UEQ0H+g5lUNJOlIDx5s+yZDhW/9VmEijNLffdEEsTA=; b=NsUNe3FFzcCNELOvsHobjV8Tx+6igsvAqt7PIBdJ7RMV7kYykrWKqglTkQKsp1Mwp8 X0HMZzoBZgMc+/ZFgkQRR/AQ78jsn8kxtu87b5TBK5IqW7ASblppGWwGQ93tqrgDPds+ y9j4+29vq57og4pms7jCEkPihI8X0WSqbRTj+pEpt4+rMXo3dsDuLiuve6ZdWq54Yhch b9nEWG5WPp0EfNBHWyUFWUjtG4i6r3xb5KpDgL4MVt1OyLRJYrTBfZlK1xrcdMB4U7pY ZMdy/IJawmSyAWbeIP9hFo7JpOjtOSjKrO6csQYdVDxcrP0XxJHhyZ7hWF1QFAeeilb8 PzeQ== X-Gm-Message-State: AMke39nzWuaWdBrhmu7ZiDfV9U2GOkP+rgE4I2kHP1uSXYxX5yNOTgF4wVf7fQsVS/QiJQXPgRV+5wlLdn3uQw== X-Received: by 10.107.174.219 with SMTP id n88mr7994281ioo.132.1486551819632; Wed, 08 Feb 2017 03:03:39 -0800 (PST) MIME-Version: 1.0 Received: by 10.107.146.134 with HTTP; Wed, 8 Feb 2017 03:03:39 -0800 (PST) From: mokhi Date: Wed, 8 Feb 2017 14:33:39 +0330 Message-ID: Subject: C11 conformance of casinl-like functions. To: 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 Feb 2017 11:03:40 -0000 Hi. I recently saw bunch of PRs opened about C11 lack of conformance in FreeBSD on Bugzilla, complaining cosinl, acosinl, ... not implemented. I've searched about these and I found documents about them too[2][3]. Do you think we should implement them? Or standards doesn't have suggestions on these? If yes (you think we should implement them), would you suggest simply `strong aliasing symbols of FUNC_l to FUNC`? or implementing FUNC_l from scratch? I've made some patch based on my idea of aliasing symbols, If you agree I like to work on this and probably can start a review on phabricator for that. Thanks and best wishes, Mokhi. ========================================== [2] https://www.gnu.org/software/gnulib/manual/gnulib.html [3] http://pubs.opengroup.org/onlinepubs/9699919799/functions/casinhl.html and many more like this on opengroup. From owner-freebsd-numerics@freebsd.org Wed Feb 8 11:38:48 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 551C5CD6C8F for ; Wed, 8 Feb 2017 11:38:48 +0000 (UTC) (envelope-from brde@optusnet.com.au) Received: from mail105.syd.optusnet.com.au (mail105.syd.optusnet.com.au [211.29.132.249]) by mx1.freebsd.org (Postfix) with ESMTP id 2059C1097 for ; Wed, 8 Feb 2017 11:38:48 +0000 (UTC) (envelope-from brde@optusnet.com.au) Received: from besplex.bde.org (c122-106-153-191.carlnfd1.nsw.optusnet.com.au [122.106.153.191]) by mail105.syd.optusnet.com.au (Postfix) with ESMTPS id 779261048403; Wed, 8 Feb 2017 22:38:40 +1100 (AEDT) Date: Wed, 8 Feb 2017 22:38:39 +1100 (EST) From: Bruce Evans X-X-Sender: bde@besplex.bde.org To: mokhi cc: freebsd-numerics@freebsd.org Subject: Re: C11 conformance of casinl-like functions. In-Reply-To: Message-ID: <20170208221449.K14261@besplex.bde.org> References: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed X-Optus-CM-Score: 0 X-Optus-CM-Analysis: v=2.2 cv=BKLDlBYG c=1 sm=1 tr=0 a=Tj3pCpwHnMupdyZSltBt7Q==:117 a=Tj3pCpwHnMupdyZSltBt7Q==:17 a=kj9zAlcOel0A:10 a=6I5d2MoRAAAA:8 a=ofqOCum91nHaNWPVwLQA:9 a=CjuIK1q_8ugA:10 a=IjZwj45LgO3ly-622nXo:22 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 Feb 2017 11:38:48 -0000 On Wed, 8 Feb 2017, mokhi wrote: > I recently saw bunch of PRs opened about C11 lack of conformance in > FreeBSD on Bugzilla, complaining cosinl, acosinl, ... not implemented. I think you mean acosl, asinl, ... These were implemented quite well in 2012-2013, but not quite finished, and not committed. Only the float and double version were committed. The raw versions are still available in https://people.freebsd.org/~stephen/catrig*.c. These have rotted and require some editing. Compare with the committed parts to see most of the necessary editing. > I've searched about these and I found documents about them too[2][3]. > Do you think we should implement them? Or standards doesn't have > suggestions on these? > > If yes (you think we should implement them), would you suggest simply > `strong aliasing symbols of FUNC_l to FUNC`? They cannot be implemented by aliasing (except on arches which don't have real long doubles). Bad versions can be implemented using wrappers. Bad versions are worse than none IMO. > or implementing FUNC_l > from scratch? It would be hard to do better than the 2012-2013 implementation. Bruce From owner-freebsd-numerics@freebsd.org Wed Feb 8 12:47:03 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 10385CD6B3D for ; Wed, 8 Feb 2017 12:47:03 +0000 (UTC) (envelope-from mokhi64@gmail.com) Received: from mail-io0-x22c.google.com (mail-io0-x22c.google.com [IPv6:2607:f8b0:4001:c06::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 C26DD18DC for ; Wed, 8 Feb 2017 12:47:02 +0000 (UTC) (envelope-from mokhi64@gmail.com) Received: by mail-io0-x22c.google.com with SMTP id v96so114836706ioi.0 for ; Wed, 08 Feb 2017 04:47:02 -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=puaUB8qBI1pLtgcspESG2lzDx9BCa6kspAF2JdRcsKc=; b=UiK48ioSQSzNTv1/LtDscql9ThApdnLqhuTke/GTn5vSS65QYLeC/w8k24gQYOCvzB cxPbubpNJvqkTH38ZIulCe+XaNUgf0TZC4oCVDz/HGhLUdc5Dj+AoUTGTaE/gw2Itxi/ /9UPvoQXTeuOGFqVLsYOiqe9CFWKk8ZOp3TRw2gEC7ekigWShirh5+jmierLgunz8Wst q1JhBVQRSDAsLvHlsk0HRUtSuJEVBBvljl6+BaJP4BJis6oY1D0kcHi9IhbYrYqdrERj gHxwOoezn5erQqaOOHtKWQqRBXqMxB0mnrpjc/Tgw0aADvggMw2bCivholthvDA72vCO BBhA== 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=puaUB8qBI1pLtgcspESG2lzDx9BCa6kspAF2JdRcsKc=; b=CHqjSpl8SDMSoMQLh1uXvPVAvtEp3uVZ8UXv0Vh6bkqdpaAmem54WIdSDfzZJK6KRZ acAj+de+6fvV8JIJPywOkL15S1vSVGwC+1yfwhbHzPuWO9+ODygfSkRLs6DLdylYlvNR ceQ40EyQaHYHk0o0/O/e5CkswXMJfG5PTUpBZ/uuEwfM1IqeA/eD3m4k9iAKPbzj7YaY mw8whmLtKbwaOcQZbMh8QvfvBC/omj6Mv8nl+OM+KXIwAXtlLZlEdjW8SAPQGBYU4MgX 0B60fGwfSteWgF3CauBHZPY/7/AZa9QRFF0WFaO9cHbOa8UEocvpvvv40cdfaK6Vs89/ AY/A== X-Gm-Message-State: AMke39m4xL38i71PNlgawCSMiUTdLD9RF7VowcTJuumscadi4PAaL0BG/85EVpEKYiMBBQTKm7tTIlZCjIIixw== X-Received: by 10.107.147.134 with SMTP id v128mr15000218iod.75.1486558022186; Wed, 08 Feb 2017 04:47:02 -0800 (PST) MIME-Version: 1.0 Received: by 10.107.146.134 with HTTP; Wed, 8 Feb 2017 04:47:01 -0800 (PST) In-Reply-To: <20170208221449.K14261@besplex.bde.org> References: <20170208221449.K14261@besplex.bde.org> From: mokhi Date: Wed, 8 Feb 2017 16:17:01 +0330 Message-ID: Subject: Re: C11 conformance of casinl-like functions. To: Bruce Evans 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 Feb 2017 12:47:03 -0000 > I think you mean acosl, asinl, ... yeah :D > These were implemented quite well in 2012-2013, but not quite finished, > and not committed. Only the float and double version were committed. > The raw versions are still available in > https://people.freebsd.org/~stephen/catrig*.c. These have rotted and > require some editing. Compare with the committed parts to see most > of the necessary editing. Okay, I'll do. Would you like that I add you on Phabricator for reviewing when my editing is done? Thanks and best wishes, Mokhi. From owner-freebsd-numerics@freebsd.org Wed Feb 8 14:14:00 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 C38C1CD6760 for ; Wed, 8 Feb 2017 14:14:00 +0000 (UTC) (envelope-from brde@optusnet.com.au) Received: from mail105.syd.optusnet.com.au (mail105.syd.optusnet.com.au [211.29.132.249]) by mx1.freebsd.org (Postfix) with ESMTP id 895401799 for ; Wed, 8 Feb 2017 14:13:59 +0000 (UTC) (envelope-from brde@optusnet.com.au) Received: from besplex.bde.org (c122-106-153-191.carlnfd1.nsw.optusnet.com.au [122.106.153.191]) by mail105.syd.optusnet.com.au (Postfix) with ESMTPS id 0DCF6104776E; Thu, 9 Feb 2017 01:13:57 +1100 (AEDT) Date: Thu, 9 Feb 2017 01:13:56 +1100 (EST) From: Bruce Evans X-X-Sender: bde@besplex.bde.org To: mokhi cc: freebsd-numerics@freebsd.org Subject: Re: C11 conformance of casinl-like functions. In-Reply-To: Message-ID: <20170209010941.J16975@besplex.bde.org> References: <20170208221449.K14261@besplex.bde.org> MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed X-Optus-CM-Score: 0 X-Optus-CM-Analysis: v=2.2 cv=c+HbeV1l c=1 sm=1 tr=0 a=Tj3pCpwHnMupdyZSltBt7Q==:117 a=Tj3pCpwHnMupdyZSltBt7Q==:17 a=kj9zAlcOel0A:10 a=6I5d2MoRAAAA:8 a=-w7s6v8T-rG7N3JJM4kA:9 a=CjuIK1q_8ugA:10 a=IjZwj45LgO3ly-622nXo:22 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 Feb 2017 14:14:00 -0000 On Wed, 8 Feb 2017, mokhi wrote: >> I think you mean acosl, asinl, ... > yeah :D > >> These were implemented quite well in 2012-2013, but not quite finished, >> and not committed. Only the float and double version were committed. >> The raw versions are still available in >> https://people.freebsd.org/~stephen/catrig*.c. These have rotted and >> require some editing. Compare with the committed parts to see most >> of the necessary editing. > Okay, I'll do. > Would you like that I add you on Phabricator for reviewing when my > editing is done? Only put it on Phabricator if you want me to not see it. This list was supposed to be for specialized review[er]s but doesn't have enough writers to be very useful. Bruce From owner-freebsd-numerics@freebsd.org Wed Feb 8 14:28:36 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 CB29FCD6AC7 for ; Wed, 8 Feb 2017 14:28:36 +0000 (UTC) (envelope-from mokhi64@gmail.com) Received: from mail-it0-x22b.google.com (mail-it0-x22b.google.com [IPv6:2607:f8b0:4001:c0b::22b]) (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 8A32E1E4A for ; Wed, 8 Feb 2017 14:28:36 +0000 (UTC) (envelope-from mokhi64@gmail.com) Received: by mail-it0-x22b.google.com with SMTP id 203so101209358ith.0 for ; Wed, 08 Feb 2017 06:28:36 -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=H9ZaOUa8xozgiRa+GGt3LdkTaf9nU3RdNbsezEVs1yQ=; b=b9Og33/GTro4GuJnB+TBeFyh+n0rxYtrE0/mghHlA58zsBPgZ9xU2TEhpuSSLUyAan HFC+UJ5ZAWxj3KjwD61a316DOp7UdmH+il3nFt23/nxqgK0oOGetJuKu4KFFLzpTG05a NJIaaX8hgJzZtoGmxwbJ9ITeAFV649gz2jPY0nSXZM/HznHTeETwacCE90LNIVLl2yC8 /6OKndLxWIzoDV3D3eLoca6AYeEVUinam6KijnvhLKF1tOyM5XGRR+1iObm9afOJO9JE li3qtfizq3zPzX22G5xR8QY9Ee7TpiaKv+NnnWJEAaUXbpWvt0lPPcdCOlKHCW6NoEN2 Y+AQ== 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=H9ZaOUa8xozgiRa+GGt3LdkTaf9nU3RdNbsezEVs1yQ=; b=YMaCqDQW5cy8BUx8xGHlRm14cKfv6pTJJT6lf1E+0THkBKxSGNb1i6EVnzXa0otB6U 76jGpJfWv9vqTrjhK/2Hc/G3WVw2MhkOYIi7ydLjKNNoa2VQu2NaM0QNk4H7X+ibsiKV bsQDPKnl2+d4TrOZx8wcQ2p94/rhWbq+Ot1lWeWPjfR/wh3FIE9lD4dxexy/SYgsS8YA 9Vm7xDvClDWmtWp7Br88BgCHRXtSBd7Zk2cJk4c/l+51ixqFDMyvOK2cUHInlt/B66M9 4el3AfQgYzg6bW/oTSQzXQ+YSTBA8hBTd2DBmfhgPt+z4mhFbGtt+PZkXfTAtvZMpnuv 1MMQ== X-Gm-Message-State: AIkVDXIWrWPMk4PeDaDKo14tjGWkkCVdI0WQVdcSg0OjjMJ4c8pCLWm5AD3xnjdTrmbHCR+g+AXPyzF94zzOFQ== X-Received: by 10.36.33.135 with SMTP id e129mr16416715ita.9.1486564116048; Wed, 08 Feb 2017 06:28:36 -0800 (PST) MIME-Version: 1.0 Received: by 10.107.146.134 with HTTP; Wed, 8 Feb 2017 06:28:35 -0800 (PST) In-Reply-To: <20170209010941.J16975@besplex.bde.org> References: <20170208221449.K14261@besplex.bde.org> <20170209010941.J16975@besplex.bde.org> From: mokhi Date: Wed, 8 Feb 2017 17:58:35 +0330 Message-ID: Subject: Re: C11 conformance of casinl-like functions. To: Bruce Evans 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 Feb 2017 14:28:36 -0000 > Only put it on Phabricator if you want me to not see it. :))) So, would you prefer that I email it? There are some people suggest Phabricator too, So if you have no opinion against it, I'll use both ways and sync them ;-) Thanks and best wishes, Mokhi. From owner-freebsd-numerics@freebsd.org Wed Feb 8 14:32: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 2FBB4CD6B0A for ; Wed, 8 Feb 2017 14:32:08 +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 BFD9C130 for ; Wed, 8 Feb 2017 14:32:07 +0000 (UTC) (envelope-from stephen@missouri.edu) X-IronPort-Anti-Spam-Filtered: true X-IronPort-Anti-Spam-Result: =?us-ascii?q?A2HoLgBYKptY/xM40cZdFgYBAQQBAQoBA?= =?us-ascii?q?YNRYSpfAQafRB+XQgEpgkKDNgKCaEQTAQIBAQEBAQEBA18dC4JbPTwBAQEBAQE?= =?us-ascii?q?jAg0GWAEFgQkCAQgYLjIlAgQNCAEBiXAOsh2IEQEHg22IUQiCYoRUhWUFm3ABh?= =?us-ascii?q?myNc443kxM3IX5RE4R+ggB1h3IBgQsBAQE?= X-IPAS-Result: =?us-ascii?q?A2HoLgBYKptY/xM40cZdFgYBAQQBAQoBAYNRYSpfAQafRB+?= =?us-ascii?q?XQgEpgkKDNgKCaEQTAQIBAQEBAQEBA18dC4JbPTwBAQEBAQEjAg0GWAEFgQkCA?= =?us-ascii?q?QgYLjIlAgQNCAEBiXAOsh2IEQEHg22IUQiCYoRUhWUFm3ABhmyNc443kxM3IX5?= =?us-ascii?q?RE4R+ggB1h3IBgQsBAQE?= Received: from ex2-t13.um.umsystem.edu ([198.209.56.19]) by um-nip4-exch-relay.um.umsystem.edu with ESMTP; 08 Feb 2017 08:30:55 -0600 Received: from EX2-T13.um.umsystem.edu (198.209.56.19) by EX2-T13.um.umsystem.edu (198.209.56.19) with Microsoft SMTP Server (version=TLS1_2, cipher=TLS_ECDHE_RSA_WITH_AES_256_CBC_SHA384_P384) id 15.1.669.32; Wed, 8 Feb 2017 08:30:55 -0600 Received: from EX2-T13.um.umsystem.edu ([198.209.56.19]) by EX2-T13.um.umsystem.edu ([198.209.56.19]) with mapi id 15.01.0669.032; Wed, 8 Feb 2017 08:30:55 -0600 From: "Montgomery-Smith, Stephen" To: "freebsd-numerics@freebsd.org" Subject: Re: C11 conformance of casinl-like functions. Thread-Topic: C11 conformance of casinl-like functions. Thread-Index: AQHSgfsFxO+P4JscBEe/hcN6eq2YKqFfYGeAgAATGoCAAB0GAA== Date: Wed, 8 Feb 2017 14:30:55 +0000 Message-ID: <0256e890-05ca-6a1b-9635-88034c544724@missouri.edu> References: <20170208221449.K14261@besplex.bde.org> In-Reply-To: Accept-Language: en-US Content-Language: en-US X-MS-Has-Attach: X-MS-TNEF-Correlator: user-agent: Mozilla/5.0 (X11; FreeBSD amd64; rv:45.0) Gecko/20100101 Thunderbird/45.5.1 x-ms-exchange-messagesentrepresentingtype: 1 x-originating-ip: [128.206.49.160] Content-Type: text/plain; charset="Windows-1252" Content-ID: <9D138D9D14B4024EBA0806C687605728@missouri.edu> Content-Transfer-Encoding: quoted-printable 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 Feb 2017 14:32:08 -0000 On 02/08/17 06:47, mokhi wrote: >> I think you mean acosl, asinl, ... > yeah :D > >> These were implemented quite well in 2012-2013, but not quite finished, >> and not committed. Only the float and double version were committed. >> The raw versions are still available in >> https://people.freebsd.org/~stephen/catrig*.c. These have rotted and >> require some editing. Compare with the committed parts to see most >> of the necessary editing. > Okay, I'll do. > Would you like that I add you on Phabricator for reviewing when my > editing is done? > > Thanks and best wishes, Mokhi. I would love to see these functions implemented. I spent a lot of time=20 on them, and they went through a lot of testing to make sure they were=20 accurate, even in very extreme edge cases. Here are more details about the implementation. http://faculty.missouri.edu/~stephen/software/#catrig= From owner-freebsd-numerics@freebsd.org Wed Feb 8 14:35:29 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 BBB88CD6C48 for ; Wed, 8 Feb 2017 14:35:29 +0000 (UTC) (envelope-from mokhi64@gmail.com) Received: from mail-io0-x232.google.com (mail-io0-x232.google.com [IPv6:2607:f8b0:4001:c06::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 7A1EE24A for ; Wed, 8 Feb 2017 14:35:29 +0000 (UTC) (envelope-from mokhi64@gmail.com) Received: by mail-io0-x232.google.com with SMTP id v96so117219358ioi.0 for ; Wed, 08 Feb 2017 06:35:29 -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=xdItsGxGQ0k6lH9y5hDIuWwiM5cTT/zAu7iVaLEo7kk=; b=FcQgpCsA/HT3etGMaRgHVFD4sGv3X4Wd7XFvLZ0ydOD0I1PYC8r2nRcZbNq7E1WxAJ uj5SvY0aGyzO8CONrY/aYwSmaki31Jhb4ValTouiLYaNriLIL9gVK+8ArQG+G3jd/iwN +5TNSfhBVZNI8dWoFJWfFh0ajszm+a6DkKjXkKfCuV8mK+B/ndfA4mZZJ/7y3bDu7VxC ThnwSTqFZKfYrCXhFM1VEME6Qu01IRJQApxb/iKPaxF1AZp7F5vyerHam9dR+OjcVVQ6 r0TbfX/RoAMl6cgNRQzb8/LVRvURQcFPzMFFQliQZSy6ExRVlKqjJpLC/KWbSfR6MxkA W+0Q== 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=xdItsGxGQ0k6lH9y5hDIuWwiM5cTT/zAu7iVaLEo7kk=; b=P4oorLgfwtx6xb6lkitkWhiOTDQj+UaACHqLtGPT0EepRdlbRfyDQHtPqzNPTzQZTJ vzjlAYYDa/7/B+srZ6bhOj56lp1sW77PqLLFsU2ENfH4Zf7QtbGjqUFOG1z/FJN6ET7n +K3SJFc8ZMlcXr5EQLtTWkvZu3W38oPbXK4ChXibgaiCH5uJvlAaGoGc0y8KRf0BljZp 9iWjP9rATYUlLrwX8RuC/3A2BD7jtCpxOT58Gb/SEcTV+2zRVa7Iy2TpXRbAMrXbvu+f IRK6Cf1SzS7YvOBdrM8/xEO3jhJMFGU6Vnd+swnMd8QlxVSELCcnbva0QNiUuhpNFrIn 3doA== X-Gm-Message-State: AMke39ldYmy810pmFdZCm8U8iKjMuzEWeUN6g/zREJbcrEs+GwjvXEMh9M8X2IXzOWzfRcwwjVUJSSV7zGLozw== X-Received: by 10.107.3.18 with SMTP id 18mr8745140iod.234.1486564528923; Wed, 08 Feb 2017 06:35:28 -0800 (PST) MIME-Version: 1.0 Received: by 10.107.146.134 with HTTP; Wed, 8 Feb 2017 06:35:28 -0800 (PST) In-Reply-To: <0256e890-05ca-6a1b-9635-88034c544724@missouri.edu> References: <20170208221449.K14261@besplex.bde.org> <0256e890-05ca-6a1b-9635-88034c544724@missouri.edu> From: mokhi Date: Wed, 8 Feb 2017 18:05:28 +0330 Message-ID: Subject: Re: C11 conformance of casinl-like functions. To: "Montgomery-Smith, Stephen" Cc: "freebsd-numerics@freebsd.org" , Bruce Evans 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 Feb 2017 14:35:29 -0000 > I would love to see these functions implemented. I spent a lot of time > on them, and they went through a lot of testing to make sure they were > accurate, even in very extreme edge cases. Thanks for your works. I am doing based on them, and editing the things needed to be renewed __/|\__ > Here are more details about the implementation. > > http://faculty.missouri.edu/~stephen/software/#catrig Thanks I'll take a look on them. Thanks and best wishes, Mokhi. From owner-freebsd-numerics@freebsd.org Wed Feb 8 15:10:36 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 08B92CD647E for ; Wed, 8 Feb 2017 15:10:36 +0000 (UTC) (envelope-from brde@optusnet.com.au) Received: from mail108.syd.optusnet.com.au (mail108.syd.optusnet.com.au [211.29.132.59]) by mx1.freebsd.org (Postfix) with ESMTP id C4A011342 for ; Wed, 8 Feb 2017 15:10:35 +0000 (UTC) (envelope-from brde@optusnet.com.au) Received: from besplex.bde.org (c122-106-153-191.carlnfd1.nsw.optusnet.com.au [122.106.153.191]) by mail108.syd.optusnet.com.au (Postfix) with ESMTPS id 0C9741A334D; Thu, 9 Feb 2017 02:10:24 +1100 (AEDT) Date: Thu, 9 Feb 2017 02:10:23 +1100 (EST) From: Bruce Evans X-X-Sender: bde@besplex.bde.org To: "Montgomery-Smith, Stephen" cc: "freebsd-numerics@freebsd.org" Subject: Re: C11 conformance of casinl-like functions. In-Reply-To: <0256e890-05ca-6a1b-9635-88034c544724@missouri.edu> Message-ID: <20170209014517.Y17120@besplex.bde.org> References: <20170208221449.K14261@besplex.bde.org> <0256e890-05ca-6a1b-9635-88034c544724@missouri.edu> MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed X-Optus-CM-Score: 0 X-Optus-CM-Analysis: v=2.2 cv=H7qr+6Qi c=1 sm=1 tr=0 a=Tj3pCpwHnMupdyZSltBt7Q==:117 a=Tj3pCpwHnMupdyZSltBt7Q==:17 a=kj9zAlcOel0A:10 a=6I5d2MoRAAAA:8 a=KKXDUkkEAAAA:8 a=nZGPi9ZBRngCGGTXSC0A:9 a=CjuIK1q_8ugA:10 a=IjZwj45LgO3ly-622nXo:22 a=y4xSRo1Ep5AFAAv36VUM:22 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 Feb 2017 15:10:36 -0000 On Wed, 8 Feb 2017, Montgomery-Smith, Stephen wrote: > On 02/08/17 06:47, mokhi wrote: >>> I think you mean acosl, asinl, ... >> yeah :D >> >>> These were implemented quite well in 2012-2013, but not quite finished, >>> and not committed. Only the float and double version were committed. >>> The raw versions are still available in >>> https://people.freebsd.org/~stephen/catrig*.c. These have rotted and >>> require some editing. Compare with the committed parts to see most >>> of the necessary editing. >> Okay, I'll do. >> Would you like that I add you on Phabricator for reviewing when my >> editing is done? >> >> Thanks and best wishes, Mokhi. > > I would love to see these functions implemented. I spent a lot of time > on them, and they went through a lot of testing to make sure they were > accurate, even in very extreme edge cases. It was you that implemented them :-). Except for many onerous details like testing on all supported arches and committing them. > Here are more details about the implementation. > > http://faculty.missouri.edu/~stephen/software/#catrig It says that the the results are good to 4 ulps except in float precision. I dout that that is correct. Higher precisions are just less likely to get near the corner cases where the errors are larger. Corner cases are probably only near zeros and poles, but I test the real and imaginary parts separately and this gives a dubious size for an ulp (for the error relative to each part instead of relative to the absolute value of the result). This gives 1-dimensional sets of corner cases. The functions try to get a small error by this measure and mostly succeed. Bruce From owner-freebsd-numerics@freebsd.org Wed Feb 8 15:13:13 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 591ADCD65B0 for ; Wed, 8 Feb 2017 15:13:13 +0000 (UTC) (envelope-from brde@optusnet.com.au) Received: from mail105.syd.optusnet.com.au (mail105.syd.optusnet.com.au [211.29.132.249]) by mx1.freebsd.org (Postfix) with ESMTP id 19C9D1640 for ; Wed, 8 Feb 2017 15:13:13 +0000 (UTC) (envelope-from brde@optusnet.com.au) Received: from besplex.bde.org (c122-106-153-191.carlnfd1.nsw.optusnet.com.au [122.106.153.191]) by mail105.syd.optusnet.com.au (Postfix) with ESMTPS id 0CE061041D12; Thu, 9 Feb 2017 02:13:12 +1100 (AEDT) Date: Thu, 9 Feb 2017 02:13:04 +1100 (EST) From: Bruce Evans X-X-Sender: bde@besplex.bde.org To: mokhi cc: freebsd-numerics@freebsd.org Subject: Re: C11 conformance of casinl-like functions. In-Reply-To: Message-ID: <20170209021059.L17120@besplex.bde.org> References: <20170208221449.K14261@besplex.bde.org> <20170209010941.J16975@besplex.bde.org> MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII; format=flowed X-Optus-CM-Score: 0 X-Optus-CM-Analysis: v=2.2 cv=H7qr+6Qi c=1 sm=1 tr=0 a=Tj3pCpwHnMupdyZSltBt7Q==:117 a=Tj3pCpwHnMupdyZSltBt7Q==:17 a=kj9zAlcOel0A:10 a=lqDJHkaq4FsQOm_mrR8A:9 a=CjuIK1q_8ugA:10 a=WdIPYdACYwkA:10 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 Feb 2017 15:13:13 -0000 On Wed, 8 Feb 2017, mokhi wrote: >> Only put it on Phabricator if you want me to not see it. > :))) > So, would you prefer that I email it? Yes. Bruce From owner-freebsd-numerics@freebsd.org Wed Feb 8 17:28:59 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 34298CD699D for ; Wed, 8 Feb 2017 17:28:59 +0000 (UTC) (envelope-from stephen@missouri.edu) Received: from um-kip8-missouri-out.um.umsystem.edu (um-kip8-missouri-out.um.umsystem.edu [209.106.229.245]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (Client CN "um-kip5.um.umsystem.edu", Issuer "InCommon RSA Server CA" (verified OK)) by mx1.freebsd.org (Postfix) with ESMTPS id DAE5F18D0 for ; Wed, 8 Feb 2017 17:28:58 +0000 (UTC) (envelope-from stephen@missouri.edu) X-IronPort-Anti-Spam-Filtered: true X-IronPort-Anti-Spam-Result: =?us-ascii?q?A2FaBwC7VJtY/yI40cZdFgUBAQEDAQEBC?= =?us-ascii?q?QEBAYQyKl8BBo5MkReXQgEpgkKDNgKDLxQBAgEBAQEBAQEDXyiCWz08AQEBAQE?= =?us-ascii?q?BIwINBlgBBXkQAgEIDgouMiUCBA0BBwEBiXAOsj+IDgEHgz4BAQEHAQEBAQEji?= =?us-ascii?q?FGCaoREEIVlBZtwAYZsjXOON5MTNiJ+UROEfgIdgWF1h3IBgQsBAQE?= X-IPAS-Result: =?us-ascii?q?A2FaBwC7VJtY/yI40cZdFgUBAQEDAQEBCQEBAYQyKl8BBo5?= =?us-ascii?q?MkReXQgEpgkKDNgKDLxQBAgEBAQEBAQEDXyiCWz08AQEBAQEBIwINBlgBBXkQA?= =?us-ascii?q?gEIDgouMiUCBA0BBwEBiXAOsj+IDgEHgz4BAQEHAQEBAQEjiFGCaoREEIVlBZt?= =?us-ascii?q?wAYZsjXOON5MTNiJ+UROEfgIdgWF1h3IBgQsBAQE?= Received: from ex-n22.um.umsystem.edu (HELO EX2-N22.um.umsystem.edu) ([198.209.56.34]) by um-kip8-exch-relay.um.umsystem.edu with ESMTP; 08 Feb 2017 11:27:30 -0600 Received: from EX2-T13.um.umsystem.edu (198.209.56.19) by EX2-N22.um.umsystem.edu (198.209.56.34) with Microsoft SMTP Server (version=TLS1_2, cipher=TLS_ECDHE_RSA_WITH_AES_256_CBC_SHA384_P384) id 15.1.669.32; Wed, 8 Feb 2017 11:27:30 -0600 Received: from EX2-T13.um.umsystem.edu ([198.209.56.19]) by EX2-T13.um.umsystem.edu ([198.209.56.19]) with mapi id 15.01.0669.032; Wed, 8 Feb 2017 11:27:30 -0600 From: "Montgomery-Smith, Stephen" To: Bruce Evans CC: "freebsd-numerics@freebsd.org" Subject: Re: C11 conformance of casinl-like functions. Thread-Topic: C11 conformance of casinl-like functions. Thread-Index: AQHSgfsFxO+P4JscBEe/hcN6eq2YKqFfYGeAgAATGoCAAB0GAIAACwiAgAAmToA= Date: Wed, 8 Feb 2017 17:27:30 +0000 Message-ID: <2ecd42a5-7fa3-a06e-a1cf-e82961382b46@missouri.edu> References: <20170208221449.K14261@besplex.bde.org> <0256e890-05ca-6a1b-9635-88034c544724@missouri.edu> <20170209014517.Y17120@besplex.bde.org> In-Reply-To: <20170209014517.Y17120@besplex.bde.org> Accept-Language: en-US Content-Language: en-US X-MS-Has-Attach: X-MS-TNEF-Correlator: user-agent: Mozilla/5.0 (X11; FreeBSD amd64; rv:45.0) Gecko/20100101 Thunderbird/45.5.1 x-ms-exchange-messagesentrepresentingtype: 1 x-originating-ip: [128.206.49.160] Content-Type: text/plain; charset="Windows-1252" Content-ID: <4E6071F25233484E86D14EEEBAC308F4@missouri.edu> Content-Transfer-Encoding: quoted-printable 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 Feb 2017 17:28:59 -0000 On 02/08/17 09:10, Bruce Evans wrote: > On Wed, 8 Feb 2017, Montgomery-Smith, Stephen wrote: > >> On 02/08/17 06:47, mokhi wrote: >>>> I think you mean acosl, asinl, ... >>> yeah :D >>> >>>> These were implemented quite well in 2012-2013, but not quite >>>> finished, >>>> and not committed. Only the float and double version were committed. >>>> The raw versions are still available in >>>> https://people.freebsd.org/~stephen/catrig*.c. These have rotted and >>>> require some editing. Compare with the committed parts to see most >>>> of the necessary editing. >>> Okay, I'll do. >>> Would you like that I add you on Phabricator for reviewing when my >>> editing is done? >>> >>> Thanks and best wishes, Mokhi. >> >> I would love to see these functions implemented. I spent a lot of time >> on them, and they went through a lot of testing to make sure they were >> accurate, even in very extreme edge cases. > > It was you that implemented them :-). Except for many onerous details > like testing on all supported arches and committing them. > >> Here are more details about the implementation. >> >> http://faculty.missouri.edu/~stephen/software/#catrig > > It says that the the results are good to 4 ulps except in float precision= . > I dout that that is correct. Higher precisions are just less likely to > get near the corner cases where the errors are larger. > > Corner cases are probably only near zeros and poles, but I test the > real and imaginary parts separately and this gives a dubious size for > an ulp (for the error relative to each part instead of relative to > the absolute value of the result). This gives 1-dimensional sets of > corner cases. The functions try to get a small error by this measure > and mostly succeed. > > Bruce I did some very extensive testing. In fact, I was able to find an edge=20 case in which the paper by Hull, Fairgrieve and Tang had a mistake if=20 one measured the real and imaginary parts separately for cacos/cacosh.=20 I corrected this in the Boost library:=20 https://svn.boost.org/trac/boost/ticket/7290 The algorithm in the paper by Hull, Fairgrieve and Tang would compute=20 cacos(1.00000002785941 + I*5.72464869028403e-200) as (approximately) 0 - I*0.000236048349018331 whereas it should be (approximately) 2.42520172707401e-196 - I*0.000236048349018331. (This is not only close to a zero of acos, I think it is also the end of=20 a branch point.)=