From patchwork Tue Sep 19 12:52:32 2017 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit X-Patchwork-Submitter: Szabolcs Nagy X-Patchwork-Id: 815505 Return-Path: X-Original-To: incoming@patchwork.ozlabs.org Delivered-To: patchwork-incoming@bilbo.ozlabs.org Authentication-Results: ozlabs.org; spf=pass (mailfrom) smtp.mailfrom=sourceware.org (client-ip=209.132.180.131; helo=sourceware.org; envelope-from=libc-alpha-return-84722-incoming=patchwork.ozlabs.org@sourceware.org; receiver=) Authentication-Results: ozlabs.org; dkim=pass (1024-bit key; secure) header.d=sourceware.org header.i=@sourceware.org header.b="R0+D3Aog"; dkim-atps=neutral Received: from sourceware.org (server1.sourceware.org [209.132.180.131]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (No client certificate requested) by ozlabs.org (Postfix) with ESMTPS id 3xxN8B3z4cz9s7g for ; Tue, 19 Sep 2017 22:52:58 +1000 (AEST) DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; q=dns; s=default; b= uxKKVK1kuc2gIAw4kry8adzVm57iUV90u8UfUUKUeuJK2IFsjVB7p6+vW0XQk6b+ Kjr3pMgnZIQyo7wIcH9x738sO6mnJHwdSoSpbFbJWiy+5VMCsbFIcaRUPCHbXzj7 8w/I1ljHoaIyKh9DJhunNLhlP9IdTkZx5bhrF8AF+SE= DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; s=default; bh=yWyS y/37jOGO6El9M5nQ4uFK0gI=; b=R0+D3Aog8UpvPIDcuaG+nYXQHXku3Sp11dHZ e2XLL1i8A33gd6mFJPqLJhowUBlDR2YylfuDEPYQ9o2zn5uHxAx9cP4vR1fANWcj Z9W64An4M3U0kOahFGaq4jgq+83vcJdAVlYeHKOchAA4d7hyQZO2AV5o3LefH73S rqmAFS4= Received: (qmail 17397 invoked by alias); 19 Sep 2017 12:52:51 -0000 Mailing-List: contact libc-alpha-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Unsubscribe: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: libc-alpha-owner@sourceware.org Delivered-To: mailing list libc-alpha@sourceware.org Received: (qmail 17382 invoked by uid 89); 19 Sep 2017 12:52:49 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-24.6 required=5.0 tests=AWL, BAYES_00, GIT_PATCH_0, GIT_PATCH_1, GIT_PATCH_2, GIT_PATCH_3, RCVD_IN_DNSWL_NONE, SPF_HELO_PASS, SPF_PASS autolearn=ham version=3.3.2 spammy=geoffrey, exceptional, boris, ui X-HELO: EUR01-DB5-obe.outbound.protection.outlook.com Authentication-Results: spf=none (sender IP is ) smtp.mailfrom=Szabolcs.Nagy@arm.com; Message-ID: <59C11310.5000004@arm.com> Date: Tue, 19 Sep 2017 13:52:32 +0100 From: Szabolcs Nagy User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:31.0) Gecko/20100101 Thunderbird/31.8.0 MIME-Version: 1.0 To: GNU C Library CC: nd@arm.com Subject: [PATCH 1/9 v3] Optimized generic expf and exp2f with wrappers References: <59C1123F.9080003@arm.com> In-Reply-To: <59C1123F.9080003@arm.com> X-ClientProxiedBy: DB6PR07CA0196.eurprd07.prod.outlook.com (2603:10a6:6:42::26) To DB6PR0802MB2488.eurprd08.prod.outlook.com (2603:10a6:4:a0::23) X-MS-PublicTrafficType: Email X-MS-Office365-Filtering-Correlation-Id: 8d7224db-9a39-4a72-f976-08d4ff5d4c97 X-MS-Office365-Filtering-HT: Tenant X-Microsoft-Antispam: UriScan:; BCL:0; PCL:0; RULEID:(300000500095)(300135000095)(300000501095)(300135300095)(22001)(300000502095)(300135100095)(2017030254152)(48565401081)(300000503095)(300135400095)(49563074)(201703131423075)(201703031133081)(201702281549075)(300000504095)(300135200095)(300000505095)(300135600095)(300000506095)(300135500095); SRVR:DB6PR0802MB2488; X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 3:uuukppPcB650qQZcKvoomr1k1TjexcDdfcYuxvB9cQ/3hX5vR7Y9nlpVeVuRCZvRam52guwcF9l0BZun8of+QxKdN9H74UCMYtpeIMz8QKLQx3TE+Ordh2/uNlkGUKf2Tv2+iMb/zTAmGQazJL4uA7SKB+y913V0xtm63u+w8XvSBBp10l4PZyWMEE/AwpSu7FWYfDStNEvvr0e/8J0FuCzUdxJZdF0Gez3v7wMaGm9xbJcUqhhJqkbEy6chZDPC; 25:Edx5E+ErSiZxWkkzXGqU1kNnwDDAAqd0zFa6YLyvvfLIxfPDpjSpBrKi6rUapgxuRtzom9fMB7KFxXPiOMKlefwVKuXy7PFyVHjq4WGKLip9IDv1NN4xLAmy3HrOw2JOmMO0iP6xj868tsfT6MlTJBUTy8JfqMAaW02adYylDJwnfB8BF+NM7OiAEFNMPZ6uuRtN6JZ7ci3bF6KgDZvNKt8pIbyD38UG6Gl0VPkWr4B2XUWXjhWoOPjErJIIzDlf61yy7Vn8P78BjWzt1JcMG8NadFQH4RdlYnfm54hraWqL4L/Gr/yxBlY5ghCp7DEli7cHRxIe5kagDlHyjbuY+g==; 31:gQgpfFNFWPnPUmRw9QOBYiDiicCJyVucQcgRcnZM+8cUXNF6STlxoJ5ejEwqIF45P5TmT0SOxHAau6Z1Vzs4sFdxkD6U2JDjpZwgYrmr0es7w8Jkuhiv79hlWDerjdn3sMw6gcA9DN5NSoWxXiBr3kMz5OTNQZVDdnkvd11RgCBKY808s5eSqRy5VnJVDX1t60QCtWyLkiDPj0KjEi2TF6tkUCiWhw0oqc77in7cs/M= X-MS-TrafficTypeDiagnostic: DB6PR0802MB2488: NoDisclaimer: True X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 20:fdKeeJYmecNXNLpeeblBwdTiLHsbzDBNjlPJHUOv3mX/nGRuhb1BCNKdoYkE2kDxgCmmA7+zdyMP4hJ2TAK0C5UK6fTMA4NWqouiK3YVEQseJSPAjeymG6BglCUPB+2Hvh+ckPbMJO3MbZ2wcdji9pT7KgyoDFE84c4dKaa+Yj0=; 4:Z3jsAP2Taq8OitOcQylim+ECueXyTdBzf8CCzjnGebzFl0GHhnR26gQsc1YVkwWOPjKe/cIDA8Enwz2WwjoJ+TVuh3k9dZBo8g3Zueb9687LteSVZlre3xkzQNlIY9pMMS9s4ih5+Oslp+x8xHKHUYSrochpdhbVMcVYh472JGe0g5KstXEvGeP0BUNpYEgPhlUtiqjANX/muwmBwK5bZ20IHIJjMOsUKCx5TPkQc2PXtOmq29N6EFu8Ih0eZmsAEyfcRN9BFBbZr/lDCaUrlPUIfWEPzMFELnEamvRxir6jFPeeCM82T4IyEX1iJnMTLaLCCQjhuBtBVtxWNSDgtw== X-Exchange-Antispam-Report-Test: UriScan:(180628864354917)(166708455590820); X-Microsoft-Antispam-PRVS: X-Exchange-Antispam-Report-CFA-Test: BCL:0; PCL:0; RULEID:(100000700101)(100105000095)(100000701101)(100105300095)(100000702101)(100105100095)(102415395)(6040450)(2401047)(8121501046)(5005006)(100000703101)(100105400095)(10201501046)(3002001)(93006095)(93001095)(6055026)(6041248)(20161123562025)(20161123564025)(20161123555025)(201703131423075)(201702281528075)(201703061421075)(201703061406153)(20161123560025)(20161123558100)(6072148)(201708071742011)(100000704101)(100105200095)(100000705101)(100105500095); SRVR:DB6PR0802MB2488; BCL:0; PCL:0; RULEID:(100000800101)(100110000095)(100000801101)(100110300095)(100000802101)(100110100095)(100000803101)(100110400095)(100000804101)(100110200095)(100000805101)(100110500095); SRVR:DB6PR0802MB2488; X-Forefront-PRVS: 04359FAD81 X-Forefront-Antispam-Report: SFV:NSPM; SFS:(10009020)(4630300001)(7370300001)(6009001)(6049001)(346002)(39860400002)(376002)(199003)(189002)(377424004)(83506001)(25786009)(4326008)(101416001)(59896002)(50986999)(478600001)(105586002)(106356001)(966005)(2906002)(87266999)(80316001)(189998001)(65816999)(54356999)(76176999)(84326002)(6116002)(3846002)(72206003)(2476003)(270700001)(568964002)(68736007)(4610100001)(5890100001)(6486002)(77096006)(7350300001)(2950100002)(36756003)(5660300001)(7736002)(305945005)(33656002)(66066001)(53936002)(16526017)(16576012)(65956001)(316002)(65806001)(16586007)(58126008)(64126003)(97736004)(8676002)(6306002)(6916009)(21490400003)(81166006)(81156014)(86362001)(41533002); DIR:OUT; SFP:1101; SCL:1; SRVR:DB6PR0802MB2488; H:[10.2.206.69]; FPR:; SPF:None; PTR:InfoNoRecords; A:1; MX:1; LANG:en; Received-SPF: None (protection.outlook.com: arm.com does not designate permitted sender hosts) X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 23:iY1XD+zj7jBXtxWaRAs/PaXeZQh6Z+AScOxfggmMRs4Iij2xUGRZBDM3Fp/d36M0tc/w8ymO+g2pins0R8mReWELuadSEudT3aVm2k27u0AX4HQgq/uON9TZ+TbwfdXSarP5FRtNwk+U9EGkWiibUHv1gIni6ceXcO7146kCvcHFRFhuDIqgTzCL1eXKiNDtKjgjYnzhL1dpaVJ8aqIPhszc0cedcnaAfydCs4vEawW7zAWord3Lbu7zUifG9dqipux8WaR4ERJlKtRc3ILqgOri6FFtHvjjMNkGpn6f7d6cmaDP+tJtk4mzRgBwEwRVNpcP6hlfZDaJSGlcoJ7LkpmCV81Ie54ODaePakvUKaC+0FX7sOMbk8g+jLit9pA3KwOHL9MEVgju0gIJia+XYDesfNDdRUQVmn8HfIQ3/3hT1eNfQ+XcJrIRQJPB+XBhIDFPIW1UO/lMxhKbcIkKE5BcJLgiij+NwWuzA1V6+qfil/OvEyBr+F0cUUJSRJdEW+cnuDHcNCuJCwLypydShGLo0vgW/mjfJlFjfnbwo5Xa2MiGGbcBvWwtBVADfRiN+VuhMLI70NKsN1eGWOwKZAe1IZsNBpcRwvrkIzLk3uTWweYm4iTnsOY1yKJ4N0/Hu1ZOfB1ZMilfmoDhVBQJ6ccLdlOOMsqavLZ3EQBvyqmkeBYrGZM1sxGqyk5yOPiJAWiyob8Ybnv6Kvunh2wdpWHQwJEjlQWGoZM4ReNGa2DunyzXzE8HCIM2Bn0kATIsy2RKiOcKD8doPj0x/2UT6WFeHnwIhCjqVgoVvjFVEOLjUIdXJOfQ3t4Jxx2Mxbr3DKf9JrWiRTikTbDdeZyJ6R9qKjffIDAeD6s6MQYxtiqDBWikHWn/1AMCFUBZElvKr0KI+hzA1hb6b7eA1cWHWTC3216FS3x381TeVysIoveGqnmpkOTq9gpLy6Dz65Ei6P4yH3tgJttvfnLHTVxG6Whob7Gy0MMophT0gtENVLkO84H9rjF+DKAhQxRnnvnmwh437nid4Mby4zAlCy/KJ3eEUjK/IqMVJ7xJrvG+ekgn7RTEqqoZrRpkeFRcTAa3IfXGrjwrWRP5zKTo5kDF+M/YIv4sbL1jw6/XvXRAstIc+DCc5SQqIGWkonv8ydpIXFBn4GlAbFwbek7vIZVl7usQMakaxAyyNdjJVV7q0yZg7ZihJgcB1sV2IkamQW/j85WVrVnLX5C6fVvCBhscPAjV/ZQLgWglkXeTDYKKYkBK8YGjMaId87QD5iTJKS2Pjqfl+FJi26A8U05ZvtEXOwkqX1NZOFzGNINDv6HQbQM9uWqI5sVvYJwJt7AxctuW2RiBMN2XpmVB4L25bMv7B96e/hQ+ySibKDlpeubrMZsNiy9mFGHWmOpT6+fwKUQM/7YDMTRovvcoAwLoYlONR0OMDXpapEwmuKgVfTtWu8FahvssNJdFSvr7FDaYvVwm4U1QxTDzCDtU/Rw/erxvZp8r+ORc/kPy7rUQMhsgf5KzDUwTSJbzU5hEyUyq3QVZ X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 6:jo8O7IYKwYMTKnyfF1QT3bPC8uM/83EMj5e/Geyt8Ec505CiuNoV6kDDsXYm6/uULMPXavhN3ABD+eP1k7JbAMNGzMEgUrSuepJRJsQ8KeMYKX2/gtfx9tKUquDEpkwCqyI6PHkrgngNnJxQhtvs7I54vd25SpidP8S73UCN159teXCQBdPCMl0gA+m0CbIAx6crnuZtdwjD7G/M1KXMmFcZLOY9aD/+hNiy9OPNySG/Odl19PVo7Tm3Ok8Av3V1hv7piJVE0sChGxtNQLmxr7eDzz8ZiA/k4qkTdva8CscvwCITlvffxtw0NEp3DUGk1VMBY08fnJcHbfGCa6O2MA==; 5:eBiDil352u+kLp1ZtFkgvn+bxqyamKOX45bnd87Fv7zCLgfb6jakLoHyQQDEbj6ctCX1eyiNV6ysMWSzBWUjyIMmlzNu7+smiKvHjEIN4O7MhHKrbCMSqiMWDKyAjn34bOjqK+WLyrD5SqYtxxlVsw==; 24:vaAnS/azoFQ450h4HlsrI5orehqzHPRiTA1F+TUd2BinWeU/gHkXg+8rhFKHbZRi7LDJbj0jGlkCQT4w/Kf6hV6+6CgVZTR0Izs5z2ufzxw=; 7:K29Qh4cC8SYrZ6fE2D977NwX107Rkisgqt+6qHRaTsoQojxA45VHjaIHDF4+y8QOisBtTC/Tn+DlaN9+fuEjvxbx1Q2531VTTbIZSWrJ0qdoQcvUkh0eHIGwwvRT+IFQzLaoZseMEXDh7YPD1ztOkfiL5Q/tErxeOpkGh2G1qOw0FFbt0EZ4Xyjny5J+lESnbquUeaNBLMvYL3gHz8XVxXa9cMm/PcEcMG55hdMj5Fc= SpamDiagnosticOutput: 1:99 SpamDiagnosticMetadata: NSPM X-OriginatorOrg: arm.com X-MS-Exchange-CrossTenant-OriginalArrivalTime: 19 Sep 2017 12:52:34.5652 (UTC) X-MS-Exchange-CrossTenant-FromEntityHeader: Hosted X-MS-Exchange-Transport-CrossTenantHeadersStamped: DB6PR0802MB2488 Based on new expf and exp2f code from https://github.com/ARM-software/optimized-routines/ with the new expf benchmark (with wrapper, aarch64): reciprocal-throughput: 2.3x faster latency: 1.7x faster with the new expf benchmark (without wrapper, aarch64): reciprocal-throughput: 3.3x faster latency: 1.7x faster with naive ubenchmark (without wrapper, aarch64): reciprocal-throughput: 3.7x faster latency: 1.3x faster libm.so size on aarch64: .text size: -152 bytes .rodata size: -1740 bytes expf/exp2f worst case nearest rounding error: 0.502 ulp worst case non-nearest rounding error: 1 ulp Error checks are inline and errno setting is in separate tail called functions, but the wrappers are kept in this patch to handle the _LIB_VERSION==_SVID_ case. (So e.g. errno is set twice for expf calls and once __expf_finite calls now on targets where the new code is used.) Double precision arithmetics is used which is expected to be faster on most targets (including soft-float) than using single precision and it is easier to get good precision result with it. Const data is kept in a separate translation unit which complicates maintenance a bit, but is expected to give good code for literal loads on most targets and allows sharing data across expf, exp2f and powf. (This data is disabled on i386 and ia64 which have their own expf, exp2f and powf code.) Some configuration is in a new math_config.h the settings may need further discussion. Some details may need target specific tweaks: - best convert and round to int operation in the arg reduction may be different across targets. - code was optimized on fma target, optimal polynomial eval may be different without fma. - gcc does not always generate good code for fp bit representation access via unions or it may be inherently slow on some target. The libm-test-ulps will need adjustment because.. - The argument reduction ideally uses nearest rounded rint, but that is not efficient on most targets, so the polynomial can get evaluated on a wider interval in non-nearest rounding mode making 1 ulp errors common in that case. - The polynomial is evaluated such that it has 1 ulp error on negative tiny inputs with upward rounding, but in exchange the evaluation is better pipelined. v3: - Add sysdeps/m68k/m680x0/fpu/{math_errf/e_exp2f_data}.c 2017-09-19 Szabolcs Nagy * math/Makefile (type-float-routines): Add math_errf and e_exp2f_data. * sysdeps/aarch64/fpu/math_private.h (TOINT_INTRINSICS): Define. (roundtoint, converttoint): Likewise. * sysdeps/ieee754/flt-32/e_expf.c: New implementation. * sysdeps/ieee754/flt-32/e_exp2f.c: New implementation. * sysdeps/ieee754/flt-32/e_exp2f_data.c: New file. * sysdeps/ieee754/flt-32/math_config.h: New file. * sysdeps/ieee754/flt-32/math_errf.c: New file. * sysdeps/ieee754/flt-32/t_exp2f.h: Remove. * sysdeps/i386/fpu/e_exp2f_data.c: New file. * sysdeps/i386/fpu/math_errf.c: New file. * sysdeps/ia64/fpu/e_exp2f_data.c: New file. * sysdeps/ia64/fpu/math_errf.c: New file. * sysdeps/m68k/m680x0/fpu/e_exp2f_data.c: New file. * sysdeps/m68k/m680x0/fpu/math_errf.c: New file. diff --git a/math/Makefile b/math/Makefile index 0601f3ac43..04586156f8 100644 --- a/math/Makefile +++ b/math/Makefile @@ -115,7 +115,7 @@ type-double-routines := branred doasin dosincos halfulp mpa mpatan2 \ # float support type-float-suffix := f -type-float-routines := k_rem_pio2f +type-float-routines := k_rem_pio2f math_errf e_exp2f_data # _Float128 support type-float128-suffix := f128 diff --git a/sysdeps/aarch64/fpu/math_private.h b/sysdeps/aarch64/fpu/math_private.h index 807111ea5a..775237ffb2 100644 --- a/sysdeps/aarch64/fpu/math_private.h +++ b/sysdeps/aarch64/fpu/math_private.h @@ -319,6 +319,26 @@ libc_feresetround_noex_aarch64_ctx (struct rm_ctx *ctx) #define libc_feresetround_noexf_ctx libc_feresetround_noex_aarch64_ctx #define libc_feresetround_noexl_ctx libc_feresetround_noex_aarch64_ctx +/* Hack: only include the large arm_neon.h when needed. */ +#ifdef _MATH_CONFIG_H +#include + +/* ACLE intrinsics for frintn and fcvtns instructions. */ +#define TOINT_INTRINSICS 1 + +static inline double_t +roundtoint (double_t x) +{ + return vget_lane_f64 (vrndn_f64 (vld1_f64 (&x)), 0); +} + +static inline uint64_t +converttoint (double_t x) +{ + return vcvtnd_s64_f64 (x); +} +#endif + #include_next #endif diff --git a/sysdeps/i386/fpu/e_exp2f_data.c b/sysdeps/i386/fpu/e_exp2f_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/i386/fpu/e_exp2f_data.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/i386/fpu/math_errf.c b/sysdeps/i386/fpu/math_errf.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/i386/fpu/math_errf.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ia64/fpu/e_exp2f_data.c b/sysdeps/ia64/fpu/e_exp2f_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/ia64/fpu/e_exp2f_data.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ia64/fpu/math_errf.c b/sysdeps/ia64/fpu/math_errf.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/ia64/fpu/math_errf.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ieee754/flt-32/e_exp2f.c b/sysdeps/ieee754/flt-32/e_exp2f.c index 567d3ff6d0..be53e5c699 100644 --- a/sysdeps/ieee754/flt-32/e_exp2f.c +++ b/sysdeps/ieee754/flt-32/e_exp2f.c @@ -1,7 +1,6 @@ -/* Single-precision floating point 2^x. - Copyright (C) 1997-2017 Free Software Foundation, Inc. +/* Single-precision 2^x function. + Copyright (C) 2017 Free Software Foundation, Inc. This file is part of the GNU C Library. - Contributed by Geoffrey Keating The GNU C Library is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser General Public @@ -17,116 +16,73 @@ License along with the GNU C Library; if not, see . */ -/* The basic design here is from - Shmuel Gal and Boris Bachelis, "An Accurate Elementary Mathematical - Library for the IEEE Floating Point Standard", ACM Trans. Math. Soft., - 17 (1), March 1991, pp. 26-45. - It has been slightly modified to compute 2^x instead of e^x, and for - single-precision. - */ -#ifndef _GNU_SOURCE -# define _GNU_SOURCE -#endif -#include -#include -#include #include -#include -#include -#include - -#include "t_exp2f.h" - -static const float TWOM100 = 7.88860905e-31; -static const float TWO127 = 1.7014118346e+38; +#include +#include "math_config.h" + +/* +EXP2F_TABLE_BITS = 5 +EXP2F_POLY_ORDER = 3 + +ULP error: 0.502 (nearest rounding.) +Relative error: 1.69 * 2^-34 in [-1/64, 1/64] (before rounding.) +Wrong count: 168353 (all nearest rounding wrong results with fma.) +Non-nearest ULP error: 1 (rounded ULP error) +*/ + +#define N (1 << EXP2F_TABLE_BITS) +#define T __exp2f_data.tab +#define C __exp2f_data.poly +#define SHIFT __exp2f_data.shift_scaled + +static inline uint32_t +top12 (float x) +{ + return asuint (x) >> 20; +} float __ieee754_exp2f (float x) { - static const float himark = (float) FLT_MAX_EXP; - static const float lomark = (float) (FLT_MIN_EXP - FLT_MANT_DIG - 1); - - /* Check for usual case. */ - if (isless (x, himark) && isgreaterequal (x, lomark)) + uint32_t abstop; + uint64_t ki, t; + /* double_t for better performance on targets with FLT_EVAL_METHOD==2. */ + double_t kd, xd, z, r, r2, y, s; + + xd = (double_t) x; + abstop = top12 (x) & 0x7ff; + if (__glibc_unlikely (abstop >= top12 (128.0f))) { - static const float THREEp14 = 49152.0; - int tval, unsafe; - float rx, x22, result; - union ieee754_float ex2_u, scale_u; - - if (fabsf (x) < FLT_EPSILON / 4.0f) - return 1.0f + x; - - { - SET_RESTORE_ROUND_NOEXF (FE_TONEAREST); - - /* 1. Argument reduction. - Choose integers ex, -128 <= t < 128, and some real - -1/512 <= x1 <= 1/512 so that - x = ex + t/512 + x1. - - First, calculate rx = ex + t/256. */ - rx = x + THREEp14; - rx -= THREEp14; - x -= rx; /* Compute x=x1. */ - /* Compute tval = (ex*256 + t)+128. - Now, t = (tval mod 256)-128 and ex=tval/256 [that's mod, NOT %; - and /-round-to-nearest not the usual c integer /]. */ - tval = (int) (rx * 256.0f + 128.0f); - - /* 2. Adjust for accurate table entry. - Find e so that - x = ex + t/256 + e + x2 - where -7e-4 < e < 7e-4, and - (float)(2^(t/256+e)) - is accurate to one part in 2^-64. */ - - /* 'tval & 255' is the same as 'tval%256' except that it's always - positive. - Compute x = x2. */ - x -= __exp2f_deltatable[tval & 255]; - - /* 3. Compute ex2 = 2^(t/255+e+ex). */ - ex2_u.f = __exp2f_atable[tval & 255]; - tval >>= 8; - /* x2 is an integer multiple of 2^-30; avoid intermediate - underflow from the calculation of x22 * x. */ - unsafe = abs(tval) >= -FLT_MIN_EXP - 32; - ex2_u.ieee.exponent += tval >> unsafe; - scale_u.f = 1.0; - scale_u.ieee.exponent += tval - (tval >> unsafe); - - /* 4. Approximate 2^x2 - 1, using a second-degree polynomial, - with maximum error in [-2^-9 - 2^-14, 2^-9 + 2^-14] - less than 1.3e-10. */ - - x22 = (.24022656679f * x + .69314736128f) * ex2_u.f; - } - - /* 5. Return (2^x2-1) * 2^(t/512+e+ex) + 2^(t/512+e+ex). */ - result = x22 * x + ex2_u.f; - - if (!unsafe) - return result; - else - { - result *= scale_u.f; - math_check_force_underflow_nonneg (result); - return result; - } - } - /* Exceptional cases: */ - else if (isless (x, himark)) - { - if (isinf (x)) - /* e^-inf == 0, with no error. */ - return 0; - else - /* Underflow */ - return TWOM100 * TWOM100; + /* |x| >= 128 or x is nan. */ + if (asuint (x) == asuint (-INFINITY)) + return 0.0f; + if (abstop >= top12 (INFINITY)) + return x + x; + if (x > 0.0f) + return __math_oflowf (0); + if (x <= -150.0f) + return __math_uflowf (0); +#if WANT_ERRNO_UFLOW + if (x < -149.0f) + return __math_may_uflowf (0); +#endif } - else - /* Return x, if x is a NaN or Inf; or overflow, otherwise. */ - return TWO127*x; + + /* x = k/N + r with r in [-1/(2N), 1/(2N)] and int k. */ + kd = (double) (xd + SHIFT); /* Rounding to double precision is required. */ + ki = asuint64 (kd); + kd -= SHIFT; /* k/N for int k. */ + r = xd - kd; + + /* exp2(x) = 2^(k/N) * 2^r ~= s * (C0*r^3 + C1*r^2 + C2*r + 1) */ + t = T[ki % N]; + t += ki << (52 - EXP2F_TABLE_BITS); + s = asdouble (t); + z = C[0] * r + C[1]; + r2 = r * r; + y = C[2] * r + 1; + y = z * r2 + y; + y = y * s; + return (float) y; } strong_alias (__ieee754_exp2f, __exp2f_finite) diff --git a/sysdeps/ieee754/flt-32/e_exp2f_data.c b/sysdeps/ieee754/flt-32/e_exp2f_data.c new file mode 100644 index 0000000000..390dcae333 --- /dev/null +++ b/sysdeps/ieee754/flt-32/e_exp2f_data.c @@ -0,0 +1,44 @@ +/* Shared data between expf, exp2f and powf. + Copyright (C) 2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#include "math_config.h" + +#define N (1 << EXP2F_TABLE_BITS) + +const struct exp2f_data __exp2f_data = { + /* tab[i] = uint(2^(i/N)) - (i << 52-BITS) + used for computing 2^(k/N) for an int |k| < 150 N as + double(tab[k%N] + (k << 52-BITS)) */ + .tab = { +0x3ff0000000000000, 0x3fefd9b0d3158574, 0x3fefb5586cf9890f, 0x3fef9301d0125b51, +0x3fef72b83c7d517b, 0x3fef54873168b9aa, 0x3fef387a6e756238, 0x3fef1e9df51fdee1, +0x3fef06fe0a31b715, 0x3feef1a7373aa9cb, 0x3feedea64c123422, 0x3feece086061892d, +0x3feebfdad5362a27, 0x3feeb42b569d4f82, 0x3feeab07dd485429, 0x3feea47eb03a5585, +0x3feea09e667f3bcd, 0x3fee9f75e8ec5f74, 0x3feea11473eb0187, 0x3feea589994cce13, +0x3feeace5422aa0db, 0x3feeb737b0cdc5e5, 0x3feec49182a3f090, 0x3feed503b23e255d, +0x3feee89f995ad3ad, 0x3feeff76f2fb5e47, 0x3fef199bdd85529c, 0x3fef3720dcef9069, +0x3fef5818dcfba487, 0x3fef7c97337b9b5f, 0x3fefa4afa2a490da, 0x3fefd0765b6e4540, + }, + .shift_scaled = 0x1.8p+52 / N, + .poly = { 0x1.c6af84b912394p-5, 0x1.ebfce50fac4f3p-3, 0x1.62e42ff0c52d6p-1 }, + .shift = 0x1.8p+52, + .invln2_scaled = 0x1.71547652b82fep+0 * N, + .poly_scaled = { +0x1.c6af84b912394p-5/N/N/N, 0x1.ebfce50fac4f3p-3/N/N, 0x1.62e42ff0c52d6p-1/N + }, +}; diff --git a/sysdeps/ieee754/flt-32/e_expf.c b/sysdeps/ieee754/flt-32/e_expf.c index 782072f213..b8244a82cf 100644 --- a/sysdeps/ieee754/flt-32/e_expf.c +++ b/sysdeps/ieee754/flt-32/e_expf.c @@ -1,7 +1,6 @@ -/* Single-precision floating point e^x. - Copyright (C) 1997-2017 Free Software Foundation, Inc. +/* Single-precision e^x function. + Copyright (C) 2017 Free Software Foundation, Inc. This file is part of the GNU C Library. - Contributed by Geoffrey Keating The GNU C Library is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser General Public @@ -17,117 +16,87 @@ License along with the GNU C Library; if not, see . */ -/* How this works: - - The input value, x, is written as - - x = n * ln(2) + t/512 + delta[t] + x; - - where: - - n is an integer, 127 >= n >= -150; - - t is an integer, 177 >= t >= -177 - - delta is based on a table entry, delta[t] < 2^-28 - - x is whatever is left, |x| < 2^-10 - - Then e^x is approximated as - - e^x = 2^n ( e^(t/512 + delta[t]) - + ( e^(t/512 + delta[t]) - * ( p(x + delta[t] + n * ln(2)) - delta ) ) ) - - where - - p(x) is a polynomial approximating e(x)-1; - - e^(t/512 + delta[t]) is obtained from a table. - - The table used is the same one as for the double precision version; - since we have the table, we might as well use it. - - It turns out to be faster to do calculations in double precision than - to perform an 'accurate table method' expf, because of the range reduction - overhead (compare exp2f). - */ -#include -#include #include -#include -#include -#include - -extern const float __exp_deltatable[178]; -extern const double __exp_atable[355] /* __attribute__((mode(DF))) */; - -static const float TWOM100 = 7.88860905e-31; -static const float TWO127 = 1.7014118346e+38; +#include +#include "math_config.h" + +/* +EXP2F_TABLE_BITS = 5 +EXP2F_POLY_ORDER = 3 + +ULP error: 0.502 (nearest rounding.) +Relative error: 1.69 * 2^-34 in [-ln2/64, ln2/64] (before rounding.) +Wrong count: 170635 (all nearest rounding wrong results with fma.) +Non-nearest ULP error: 1 (rounded ULP error) +*/ + +#define N (1 << EXP2F_TABLE_BITS) +#define InvLn2N __exp2f_data.invln2_scaled +#define T __exp2f_data.tab +#define C __exp2f_data.poly_scaled + +static inline uint32_t +top12 (float x) +{ + return asuint (x) >> 20; +} float __ieee754_expf (float x) { - static const float himark = 88.72283935546875; - static const float lomark = -103.972084045410; - /* Check for usual case. */ - if (isless (x, himark) && isgreater (x, lomark)) + uint32_t abstop; + uint64_t ki, t; + /* double_t for better performance on targets with FLT_EVAL_METHOD==2. */ + double_t kd, xd, z, r, r2, y, s; + + xd = (double_t) x; + abstop = top12 (x) & 0x7ff; + if (__glibc_unlikely (abstop >= top12 (88.0f))) { - static const float THREEp42 = 13194139533312.0; - static const float THREEp22 = 12582912.0; - /* 1/ln(2). */ -#undef M_1_LN2 - static const float M_1_LN2 = 1.44269502163f; - /* ln(2) */ -#undef M_LN2 - static const double M_LN2 = .6931471805599452862; - - int tval; - double x22, t, result, dx; - float n, delta; - union ieee754_double ex2_u; - - { - SET_RESTORE_ROUND_NOEXF (FE_TONEAREST); - - /* Calculate n. */ - n = x * M_1_LN2 + THREEp22; - n -= THREEp22; - dx = x - n*M_LN2; - - /* Calculate t/512. */ - t = dx + THREEp42; - t -= THREEp42; - dx -= t; - - /* Compute tval = t. */ - tval = (int) (t * 512.0); - - if (t >= 0) - delta = - __exp_deltatable[tval]; - else - delta = __exp_deltatable[-tval]; - - /* Compute ex2 = 2^n e^(t/512+delta[t]). */ - ex2_u.d = __exp_atable[tval+177]; - ex2_u.ieee.exponent += (int) n; - - /* Approximate e^(dx+delta) - 1, using a second-degree polynomial, - with maximum error in [-2^-10-2^-28,2^-10+2^-28] - less than 5e-11. */ - x22 = (0.5000000496709180453 * dx + 1.0000001192102037084) * dx + delta; - } - - /* Return result. */ - result = x22 * ex2_u.d + ex2_u.d; - return (float) result; + /* |x| >= 88 or x is nan. */ + if (asuint (x) == asuint (-INFINITY)) + return 0.0f; + if (abstop >= top12 (INFINITY)) + return x + x; + if (x > 0x1.62e42ep6f) /* x > log(0x1p128) ~= 88.72 */ + return __math_oflowf (0); + if (x < -0x1.9fe368p6f) /* x < log(0x1p-150) ~= -103.97 */ + return __math_uflowf (0); +#if WANT_ERRNO_UFLOW + if (x < -0x1.9d1d9ep6f) /* x < log(0x1p-149) ~= -103.28 */ + return __math_may_uflowf (0); +#endif } - /* Exceptional cases: */ - else if (isless (x, himark)) - { - if (isinf (x)) - /* e^-inf == 0, with no error. */ - return 0; - else - /* Underflow */ - return TWOM100 * TWOM100; - } - else - /* Return x, if x is a NaN or Inf; or overflow, otherwise. */ - return TWO127*x; + + /* x*N/Ln2 = k + r with r in [-1/2, 1/2] and int k. */ + z = InvLn2N * xd; + + /* Round and convert z to int, the result is in [-150*N, 128*N] and + ideally ties-to-even rule is used, otherwise the magnitude of r + can be bigger which gives larger approximation error. */ +#if TOINT_INTRINSICS + kd = roundtoint (z); + ki = converttoint (z); +#elif TOINT_RINT + kd = rint (z); + ki = (long) kd; +#elif TOINT_SHIFT +# define SHIFT __exp2f_data.shift + kd = (double) (z + SHIFT); /* Rounding to double precision is required. */ + ki = asuint64 (kd); + kd -= SHIFT; +#endif + r = z - kd; + + /* exp(x) = 2^(k/N) * 2^(r/N) ~= s * (C0*r^3 + C1*r^2 + C2*r + 1) */ + t = T[ki % N]; + t += ki << (52 - EXP2F_TABLE_BITS); + s = asdouble (t); + z = C[0] * r + C[1]; + r2 = r * r; + y = C[2] * r + 1; + y = z * r2 + y; + y = y * s; + return (float) y; } strong_alias (__ieee754_expf, __expf_finite) diff --git a/sysdeps/ieee754/flt-32/math_config.h b/sysdeps/ieee754/flt-32/math_config.h new file mode 100644 index 0000000000..1759c447c8 --- /dev/null +++ b/sysdeps/ieee754/flt-32/math_config.h @@ -0,0 +1,114 @@ +/* Configuration for math routines. + Copyright (C) 2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#ifndef _MATH_CONFIG_H +#define _MATH_CONFIG_H + +#include +#include +#include + +#ifndef WANT_ROUNDING +/* Correct special case results in non-nearest rounding modes. */ +# define WANT_ROUNDING 1 +#endif +#ifndef WANT_ERRNO +/* Set errno according to ISO C with (math_errhandling & MATH_ERRNO) != 0. */ +# define WANT_ERRNO 1 +#endif +#ifndef WANT_ERRNO_UFLOW +/* Set errno to ERANGE if result underflows to 0 (in all rounding modes). */ +# define WANT_ERRNO_UFLOW (WANT_ROUNDING && WANT_ERRNO) +#endif + +#ifndef TOINT_INTRINSICS +# define TOINT_INTRINSICS 0 +#endif +#ifndef TOINT_RINT +# define TOINT_RINT 0 +#endif +#ifndef TOINT_SHIFT +# define TOINT_SHIFT 1 +#endif + +static inline uint32_t +asuint (float f) +{ + union + { + float f; + uint32_t i; + } u = {f}; + return u.i; +} + +static inline float +asfloat (uint32_t i) +{ + union + { + uint32_t i; + float f; + } u = {i}; + return u.f; +} + +static inline uint64_t +asuint64 (double f) +{ + union + { + double f; + uint64_t i; + } u = {f}; + return u.i; +} + +static inline double +asdouble (uint64_t i) +{ + union + { + uint64_t i; + double f; + } u = {i}; + return u.f; +} + +#define NOINLINE __attribute__ ((noinline)) + +attribute_hidden float __math_oflowf (unsigned long); +attribute_hidden float __math_uflowf (unsigned long); +attribute_hidden float __math_may_uflowf (unsigned long); +attribute_hidden float __math_divzerof (unsigned long); +attribute_hidden float __math_invalidf (float); + +/* Shared between expf, exp2f and powf. */ +#define EXP2F_TABLE_BITS 5 +#define EXP2F_POLY_ORDER 3 +extern const struct exp2f_data +{ + uint64_t tab[1 << EXP2F_TABLE_BITS]; + double shift_scaled; + double poly[EXP2F_POLY_ORDER]; + double shift; + double invln2_scaled; + double poly_scaled[EXP2F_POLY_ORDER]; +} __exp2f_data; + +#endif diff --git a/sysdeps/ieee754/flt-32/math_errf.c b/sysdeps/ieee754/flt-32/math_errf.c new file mode 100644 index 0000000000..ab546e24cf --- /dev/null +++ b/sysdeps/ieee754/flt-32/math_errf.c @@ -0,0 +1,76 @@ +/* Single-precision math error handling. + Copyright (C) 2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#include "math_config.h" + +#if WANT_ERRNO +# include +/* NOINLINE reduces code size. */ +NOINLINE static float +with_errnof (float y, int e) +{ + errno = e; + return y; +} +#else +# define with_errnof(x, e) (x) +#endif + +/* NOINLINE prevents fenv semantics breaking optimizations. */ +NOINLINE static float +xflowf (unsigned long sign, float y) +{ + y = (sign ? -y : y) * y; + return with_errnof (y, ERANGE); +} + +attribute_hidden float +__math_uflowf (unsigned long sign) +{ + return xflowf (sign, 0x1p-95f); +} + +#if WANT_ERRNO_UFLOW +/* Underflows to zero in some non-nearest rounding mode, setting errno + is valid even if the result is non-zero, but in the subnormal range. */ +attribute_hidden float +__math_may_uflowf (unsigned long sign) +{ + return xflowf (sign, 0x1.4p-75f); +} +#endif + +attribute_hidden float +__math_oflowf (unsigned long sign) +{ + return xflowf (sign, 0x1p97f); +} + +attribute_hidden float +__math_divzerof (unsigned long sign) +{ + float y = 0; + return with_errnof ((sign ? -1 : 1) / y, ERANGE); +} + +attribute_hidden float +__math_invalidf (float x) +{ + float y = (x - x) / (x - x); + return isnan (x) ? y : with_errnof (y, EDOM); +} diff --git a/sysdeps/ieee754/flt-32/t_exp2f.h b/sysdeps/ieee754/flt-32/t_exp2f.h deleted file mode 100644 index aecabcc372..0000000000 --- a/sysdeps/ieee754/flt-32/t_exp2f.h +++ /dev/null @@ -1,351 +0,0 @@ -/* Accurate tables for exp2f(). - Copyright (C) 1998-2017 Free Software Foundation, Inc. - This file is part of the GNU C Library. - Contributed by Geoffrey Keating - - The GNU C Library is free software; you can redistribute it and/or - modify it under the terms of the GNU Lesser General Public - License as published by the Free Software Foundation; either - version 2.1 of the License, or (at your option) any later version. - - The GNU C Library is distributed in the hope that it will be useful, - but WITHOUT ANY WARRANTY; without even the implied warranty of - MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU - Lesser General Public License for more details. - - You should have received a copy of the GNU Lesser General Public - License along with the GNU C Library; if not, see - . */ - -/* This table has the property that, for all integers -128 <= i <= 127, - exp(i/256.0 + __exp2f_deltatable[i-128]) == __exp2f_atable[i+128] + r - for some -2^-35 < r < 2^-35 (abs(r) < 2^-36 if i <= 0); and that - __exp2f_deltatable[i+128] == t * 2^-30 - for integer t so that abs(t) <= 43447 * 2^0. */ - -#define W30 (9.31322575e-10) -static const float __exp2f_deltatable[256] = { - -810*W30, 283*W30, -1514*W30, 1304*W30, - -1148*W30, -98*W30, -744*W30, -156*W30, - -419*W30, -155*W30, 474*W30, 167*W30, - -1984*W30, -826*W30, 692*W30, 781*W30, - -578*W30, -411*W30, -129*W30, -1500*W30, - 654*W30, -141*W30, -816*W30, -53*W30, - 148*W30, 493*W30, -2214*W30, 760*W30, - 260*W30, 750*W30, -1300*W30, 1424*W30, - -1445*W30, -339*W30, -680*W30, -349*W30, - -922*W30, 531*W30, 193*W30, -2892*W30, - 290*W30, -2145*W30, -276*W30, 485*W30, - -695*W30, 215*W30, -7093*W30, 412*W30, - -4596*W30, 367*W30, 592*W30, -615*W30, - -97*W30, -1066*W30, 972*W30, -226*W30, - -625*W30, -374*W30, -5647*W30, -180*W30, - 20349*W30, -447*W30, 111*W30, -4164*W30, - -87*W30, -21*W30, -251*W30, 66*W30, - -517*W30, 2093*W30, -263*W30, 182*W30, - -601*W30, 475*W30, -483*W30, -1251*W30, - -373*W30, 1471*W30, -92*W30, -215*W30, - -97*W30, -190*W30, 0*W30, -290*W30, - -2647*W30, 1940*W30, -582*W30, 28*W30, - 833*W30, 1493*W30, 34*W30, 321*W30, - 3327*W30, -35*W30, 177*W30, -135*W30, - -796*W30, -428*W30, 129*W30, 9332*W30, - -12*W30, -69*W30, -1743*W30, 6508*W30, - -60*W30, 359*W30, 43447*W30, 15*W30, - -23*W30, -305*W30, -375*W30, -652*W30, - 667*W30, 269*W30, -1575*W30, 185*W30, - -329*W30, 200*W30, 6002*W30, 163*W30, - -647*W30, 19*W30, -603*W30, -755*W30, - 742*W30, -438*W30, 3587*W30, 2560*W30, - 0*W30, -520*W30, -241*W30, -299*W30, - -1270*W30, -991*W30, -1138*W30, 255*W30, - -1192*W30, 1722*W30, 1023*W30, 3700*W30, - -1388*W30, -1551*W30, -2549*W30, 27*W30, - 282*W30, 673*W30, 113*W30, 1561*W30, - 72*W30, 873*W30, 87*W30, -395*W30, - -433*W30, 629*W30, 3440*W30, -284*W30, - -592*W30, -103*W30, -46*W30, -3844*W30, - 1712*W30, 303*W30, 1555*W30, -631*W30, - -1400*W30, -961*W30, -854*W30, -276*W30, - 407*W30, 833*W30, -345*W30, -1501*W30, - 121*W30, -1581*W30, 400*W30, 150*W30, - 1224*W30, -139*W30, -563*W30, 879*W30, - 933*W30, 2939*W30, 788*W30, 211*W30, - 530*W30, -192*W30, 706*W30, -13347*W30, - 1065*W30, 3*W30, 111*W30, -208*W30, - -360*W30, -532*W30, -291*W30, 483*W30, - 987*W30, -33*W30, -1373*W30, -166*W30, - -1174*W30, -3955*W30, 1601*W30, -280*W30, - 1405*W30, 600*W30, -1659*W30, -23*W30, - 390*W30, 449*W30, 570*W30, -13143*W30, - -9*W30, -1646*W30, 1201*W30, 294*W30, - 2181*W30, -1173*W30, 1388*W30, -4504*W30, - 190*W30, -2304*W30, 211*W30, 239*W30, - 48*W30, -817*W30, 1018*W30, 1828*W30, - -663*W30, 1408*W30, 408*W30, -36*W30, - 1295*W30, -230*W30, 1341*W30, 9*W30, - 40*W30, 705*W30, 186*W30, 376*W30, - 557*W30, 5866*W30, 363*W30, -1558*W30, - 718*W30, 669*W30, 1369*W30, -2972*W30, - -468*W30, -121*W30, -219*W30, 667*W30, - 29954*W30, 366*W30, 48*W30, -203*W30 -}; - -static const float __exp2f_atable[256] /* __attribute__((mode(SF))) */ = { - 0.707106411447, /* 0x0.b504ecfff */ - 0.709024071690, /* 0x0.b58299fff */ - 0.710945606239, /* 0x0.b60088000 */ - 0.712874472142, /* 0x0.b67ef1000 */ - 0.714806139464, /* 0x0.b6fd88fff */ - 0.716744661340, /* 0x0.b77c94000 */ - 0.718687653549, /* 0x0.b7fbea000 */ - 0.720636486992, /* 0x0.b87ba1fff */ - 0.722590208040, /* 0x0.b8fbabfff */ - 0.724549472323, /* 0x0.b97c12fff */ - 0.726514220228, /* 0x0.b9fcd5fff */ - 0.728483855735, /* 0x0.ba7deb000 */ - 0.730457961549, /* 0x0.baff4afff */ - 0.732438981522, /* 0x0.bb811efff */ - 0.734425544748, /* 0x0.bc0350000 */ - 0.736416816713, /* 0x0.bc85d0000 */ - 0.738412797450, /* 0x0.bd089efff */ - 0.740414917465, /* 0x0.bd8bd4fff */ - 0.742422521111, /* 0x0.be0f66fff */ - 0.744434773914, /* 0x0.be9346fff */ - 0.746454179287, /* 0x0.bf179f000 */ - 0.748477637755, /* 0x0.bf9c3afff */ - 0.750506639473, /* 0x0.c02133fff */ - 0.752541840064, /* 0x0.c0a694fff */ - 0.754582285889, /* 0x0.c12c4e000 */ - 0.756628334525, /* 0x0.c1b265000 */ - 0.758678436269, /* 0x0.c238bffff */ - 0.760736882681, /* 0x0.c2bfa6fff */ - 0.762799203401, /* 0x0.c346cf000 */ - 0.764867603790, /* 0x0.c3ce5d000 */ - 0.766940355298, /* 0x0.c45633fff */ - 0.769021093841, /* 0x0.c4de90fff */ - 0.771104693409, /* 0x0.c5671dfff */ - 0.773195922364, /* 0x0.c5f02afff */ - 0.775292098512, /* 0x0.c6798afff */ - 0.777394294745, /* 0x0.c70350000 */ - 0.779501736166, /* 0x0.c78d6d000 */ - 0.781615912910, /* 0x0.c817fafff */ - 0.783734917628, /* 0x0.c8a2d9fff */ - 0.785858273516, /* 0x0.c92e02000 */ - 0.787990570071, /* 0x0.c9b9c0000 */ - 0.790125787245, /* 0x0.ca45aefff */ - 0.792268991467, /* 0x0.cad223fff */ - 0.794417440881, /* 0x0.cb5ef0fff */ - 0.796570718287, /* 0x0.cbec0efff */ - 0.798730909811, /* 0x0.cc79a0fff */ - 0.800892710672, /* 0x0.cd074dfff */ - 0.803068041795, /* 0x0.cd95ddfff */ - 0.805242776881, /* 0x0.ce2464000 */ - 0.807428598393, /* 0x0.ceb3a3fff */ - 0.809617877002, /* 0x0.cf431dfff */ - 0.811812341211, /* 0x0.cfd2eefff */ - 0.814013659956, /* 0x0.d06333000 */ - 0.816220164311, /* 0x0.d0f3ce000 */ - 0.818434238424, /* 0x0.d184e7fff */ - 0.820652604094, /* 0x0.d21649fff */ - 0.822877407074, /* 0x0.d2a818000 */ - 0.825108587751, /* 0x0.d33a51000 */ - 0.827342867839, /* 0x0.d3ccbdfff */ - 0.829588949684, /* 0x0.d45ff1000 */ - 0.831849217401, /* 0x0.d4f411fff */ - 0.834093391880, /* 0x0.d58724fff */ - 0.836355149750, /* 0x0.d61b5f000 */ - 0.838620424257, /* 0x0.d6afd3fff */ - 0.840896368027, /* 0x0.d744fc000 */ - 0.843176305293, /* 0x0.d7da66fff */ - 0.845462262643, /* 0x0.d87037000 */ - 0.847754716864, /* 0x0.d90673fff */ - 0.850052893157, /* 0x0.d99d10fff */ - 0.852359056469, /* 0x0.da3433fff */ - 0.854668736446, /* 0x0.dacb91fff */ - 0.856986224651, /* 0x0.db6373000 */ - 0.859309315673, /* 0x0.dbfbb1fff */ - 0.861639738080, /* 0x0.dc946bfff */ - 0.863975346095, /* 0x0.dd2d7d000 */ - 0.866317391394, /* 0x0.ddc6f9fff */ - 0.868666708472, /* 0x0.de60f1000 */ - 0.871022939695, /* 0x0.defb5c000 */ - 0.873383641229, /* 0x0.df9611fff */ - 0.875751554968, /* 0x0.e03141000 */ - 0.878126025200, /* 0x0.e0ccde000 */ - 0.880506813521, /* 0x0.e168e4fff */ - 0.882894217966, /* 0x0.e2055afff */ - 0.885287821299, /* 0x0.e2a239000 */ - 0.887686729423, /* 0x0.e33f6ffff */ - 0.890096127973, /* 0x0.e3dd56fff */ - 0.892507970338, /* 0x0.e47b67000 */ - 0.894928157336, /* 0x0.e51a03000 */ - 0.897355020043, /* 0x0.e5b90efff */ - 0.899788379682, /* 0x0.e65888000 */ - 0.902227103705, /* 0x0.e6f85afff */ - 0.904673457151, /* 0x0.e798ae000 */ - 0.907128036008, /* 0x0.e8398afff */ - 0.909585535528, /* 0x0.e8da99000 */ - 0.912051796915, /* 0x0.e97c3a000 */ - 0.914524436003, /* 0x0.ea1e46000 */ - 0.917003571999, /* 0x0.eac0bf000 */ - 0.919490039339, /* 0x0.eb63b2fff */ - 0.921983361257, /* 0x0.ec071a000 */ - 0.924488604054, /* 0x0.ecab48fff */ - 0.926989555360, /* 0x0.ed4f30000 */ - 0.929502844812, /* 0x0.edf3e6000 */ - 0.932021975503, /* 0x0.ee98fdfff */ - 0.934553921208, /* 0x0.ef3eecfff */ - 0.937083780759, /* 0x0.efe4b8fff */ - 0.939624726786, /* 0x0.f08b3f000 */ - 0.942198514924, /* 0x0.f133ebfff */ - 0.944726586343, /* 0x0.f1d99a000 */ - 0.947287976728, /* 0x0.f28176fff */ - 0.949856162070, /* 0x0.f329c5fff */ - 0.952431440345, /* 0x0.f3d28bfff */ - 0.955013573175, /* 0x0.f47bc5000 */ - 0.957603693021, /* 0x0.f52584000 */ - 0.960199773321, /* 0x0.f5cfa7000 */ - 0.962801992906, /* 0x0.f67a31000 */ - 0.965413510788, /* 0x0.f72556fff */ - 0.968030691152, /* 0x0.f7d0dc000 */ - 0.970655620084, /* 0x0.f87ce2fff */ - 0.973290979849, /* 0x0.f92998fff */ - 0.975926160805, /* 0x0.f9d64bfff */ - 0.978571653370, /* 0x0.fa83ac000 */ - 0.981225252139, /* 0x0.fb3193fff */ - 0.983885228626, /* 0x0.fbdfe6fff */ - 0.986552715296, /* 0x0.fc8eb7fff */ - 0.989228487027, /* 0x0.fd3e14000 */ - 0.991909801964, /* 0x0.fdedcd000 */ - 0.994601726545, /* 0x0.fe9e38000 */ - 0.997297704209, /* 0x0.ff4ee6fff */ - 1.000000000000, /* 0x1.000000000 */ - 1.002710938457, /* 0x1.00b1aa000 */ - 1.005429744692, /* 0x1.0163d7ffe */ - 1.008155703526, /* 0x1.02167dffe */ - 1.010888457284, /* 0x1.02c995fff */ - 1.013629436498, /* 0x1.037d38000 */ - 1.016377568250, /* 0x1.043152000 */ - 1.019134163841, /* 0x1.04e5f9ffe */ - 1.021896362316, /* 0x1.059b00000 */ - 1.024668931945, /* 0x1.0650b3ffe */ - 1.027446627635, /* 0x1.0706be001 */ - 1.030234098408, /* 0x1.07bd6bffe */ - 1.033023953416, /* 0x1.087441ffe */ - 1.035824656494, /* 0x1.092bce000 */ - 1.038632392900, /* 0x1.09e3d0001 */ - 1.041450142840, /* 0x1.0a9c79ffe */ - 1.044273972530, /* 0x1.0b558a001 */ - 1.047105550795, /* 0x1.0c0f1c001 */ - 1.049944162390, /* 0x1.0cc924001 */ - 1.052791833895, /* 0x1.0d83c4001 */ - 1.055645227426, /* 0x1.0e3ec3fff */ - 1.058507919326, /* 0x1.0efa60001 */ - 1.061377286898, /* 0x1.0fb66bfff */ - 1.064254641510, /* 0x1.1072fdffe */ - 1.067140102389, /* 0x1.113018000 */ - 1.070034146304, /* 0x1.11edc1fff */ - 1.072937250162, /* 0x1.12ac04001 */ - 1.075843691823, /* 0x1.136a7dfff */ - 1.078760385496, /* 0x1.1429a3ffe */ - 1.081685543070, /* 0x1.14e958000 */ - 1.084618330005, /* 0x1.15a98c000 */ - 1.087556362176, /* 0x1.166a18001 */ - 1.090508937863, /* 0x1.172b98001 */ - 1.093464612954, /* 0x1.17ed4bfff */ - 1.096430182434, /* 0x1.18afa5ffe */ - 1.099401354802, /* 0x1.19725e000 */ - 1.102381587017, /* 0x1.1a35adfff */ - 1.105370759965, /* 0x1.1af994000 */ - 1.108367800686, /* 0x1.1bbdfdffe */ - 1.111373305331, /* 0x1.1c82f6000 */ - 1.114387035385, /* 0x1.1d4878001 */ - 1.117408752440, /* 0x1.1e0e7ffff */ - 1.120437502874, /* 0x1.1ed4fe000 */ - 1.123474478729, /* 0x1.1f9c06000 */ - 1.126521706601, /* 0x1.2063ba001 */ - 1.129574775716, /* 0x1.212bd0001 */ - 1.132638812065, /* 0x1.21f49e000 */ - 1.135709524130, /* 0x1.22bddbffe */ - 1.138789534565, /* 0x1.2387b5fff */ - 1.141876101508, /* 0x1.2451fe000 */ - 1.144971728301, /* 0x1.251cddffe */ - 1.148077130296, /* 0x1.25e861ffe */ - 1.151189923305, /* 0x1.26b462001 */ - 1.154312610610, /* 0x1.278107ffe */ - 1.157440662410, /* 0x1.284e08001 */ - 1.160578370109, /* 0x1.291baa001 */ - 1.163725256932, /* 0x1.29e9e6000 */ - 1.166879892324, /* 0x1.2ab8a3ffe */ - 1.170044302935, /* 0x1.2b8805fff */ - 1.173205971694, /* 0x1.2c5739ffe */ - 1.176397800428, /* 0x1.2d2867ffe */ - 1.179586529747, /* 0x1.2df962001 */ - 1.182784795737, /* 0x1.2ecafbffe */ - 1.185991406414, /* 0x1.2f9d21ffe */ - 1.189206838636, /* 0x1.306fdc001 */ - 1.192430973067, /* 0x1.314328000 */ - 1.195664167430, /* 0x1.32170c001 */ - 1.198906540890, /* 0x1.32eb8a001 */ - 1.202157497408, /* 0x1.33c098000 */ - 1.205416083326, /* 0x1.349625fff */ - 1.208683252332, /* 0x1.356c43fff */ - 1.211961269402, /* 0x1.364318001 */ - 1.215246438983, /* 0x1.371a64000 */ - 1.218539118740, /* 0x1.37f22dffe */ - 1.221847295770, /* 0x1.38cafc000 */ - 1.225158572187, /* 0x1.39a3fdfff */ - 1.228481650325, /* 0x1.3a7dc5ffe */ - 1.231811761846, /* 0x1.3b5803fff */ - 1.235149741144, /* 0x1.3c32c5ffe */ - 1.238499879811, /* 0x1.3d0e53ffe */ - 1.241858124726, /* 0x1.3dea69fff */ - 1.245225191102, /* 0x1.3ec713fff */ - 1.248601436624, /* 0x1.3fa458000 */ - 1.251975655584, /* 0x1.40817a001 */ - 1.255380749731, /* 0x1.4160a2001 */ - 1.258783102010, /* 0x1.423f9bffe */ - 1.262198328973, /* 0x1.431f6e000 */ - 1.265619754780, /* 0x1.43ffa7fff */ - 1.269052743928, /* 0x1.44e0a4001 */ - 1.272490739830, /* 0x1.45c1f4000 */ - 1.275942921659, /* 0x1.46a432001 */ - 1.279397487615, /* 0x1.478697ffe */ - 1.282870173427, /* 0x1.486a2dffe */ - 1.286346316319, /* 0x1.494dfdffe */ - 1.289836049094, /* 0x1.4a32b2001 */ - 1.293333172770, /* 0x1.4b17e1ffe */ - 1.296839594835, /* 0x1.4bfdadfff */ - 1.300354957560, /* 0x1.4ce40fffe */ - 1.303882122055, /* 0x1.4dcb38001 */ - 1.307417988757, /* 0x1.4eb2f1ffe */ - 1.310960650439, /* 0x1.4f9b1dfff */ - 1.314516782746, /* 0x1.50842bfff */ - 1.318079948424, /* 0x1.516daffff */ - 1.321653246888, /* 0x1.5257de000 */ - 1.325237751030, /* 0x1.5342c8001 */ - 1.328829526907, /* 0x1.542e2c000 */ - 1.332433700535, /* 0x1.551a5fffe */ - 1.336045145966, /* 0x1.56070dffe */ - 1.339667558645, /* 0x1.56f473ffe */ - 1.343300342533, /* 0x1.57e287ffe */ - 1.346941947961, /* 0x1.58d130001 */ - 1.350594043714, /* 0x1.59c087ffe */ - 1.354256033883, /* 0x1.5ab085fff */ - 1.357932448365, /* 0x1.5ba175ffe */ - 1.361609339707, /* 0x1.5c926dfff */ - 1.365299344044, /* 0x1.5d8441ffe */ - 1.369003057507, /* 0x1.5e76fc001 */ - 1.372714757920, /* 0x1.5f6a3c000 */ - 1.376437187179, /* 0x1.605e2fffe */ - 1.380165219333, /* 0x1.615282001 */ - 1.383909463864, /* 0x1.6247e3ffe */ - 1.387661933907, /* 0x1.633dd0000 */ - 1.391424179060, /* 0x1.64345fffe */ - 1.395197510706, /* 0x1.652ba9fff */ - 1.399006724329, /* 0x1.66254dffe */ - 1.402773022651, /* 0x1.671c22000 */ - 1.406576037403, /* 0x1.68155dfff */ - 1.410389423392, /* 0x1.690f48001 */ -}; diff --git a/sysdeps/m68k/m680x0/fpu/e_exp2f_data.c b/sysdeps/m68k/m680x0/fpu/e_exp2f_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/e_exp2f_data.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/m68k/m680x0/fpu/math_errf.c b/sysdeps/m68k/m680x0/fpu/math_errf.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/math_errf.c @@ -0,0 +1 @@ +/* Not needed. */ From patchwork Tue Sep 19 12:56:23 2017 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit X-Patchwork-Submitter: Szabolcs Nagy X-Patchwork-Id: 815508 Return-Path: X-Original-To: incoming@patchwork.ozlabs.org Delivered-To: patchwork-incoming@bilbo.ozlabs.org Authentication-Results: ozlabs.org; spf=pass (mailfrom) smtp.mailfrom=sourceware.org (client-ip=209.132.180.131; helo=sourceware.org; envelope-from=libc-alpha-return-84724-incoming=patchwork.ozlabs.org@sourceware.org; receiver=) Authentication-Results: ozlabs.org; dkim=pass (1024-bit key; secure) header.d=sourceware.org header.i=@sourceware.org header.b="XIgv6kaa"; dkim-atps=neutral Received: from sourceware.org (server1.sourceware.org [209.132.180.131]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (No client certificate requested) by ozlabs.org (Postfix) with ESMTPS id 3xxNDS1TPwz9s7g for ; Tue, 19 Sep 2017 22:56:40 +1000 (AEST) DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; q=dns; s=default; b= YRdzKL7h/hrNH/c/1RC2okTiJ3qJfSqy3srrP+6y8LwHDwlSVlpafgY/7xCQ2WvV KixE9llEJcD+XCTrmyHULERW4RdE2DQzeoEeu5u6X5coGN9qitOFpmIN4ZZz+AQ+ o72oQ2TLX7RPqrBwxOWCnZ2HvnoOgDjAv8n+V18s1xA= DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; s=default; bh=A8Ff Rj8uvpEXH9QHk781JtUHaOg=; b=XIgv6kaazcU2m6A3r3WINKMvk+VVHb8IRBUh A9Zo7QoVveHp5VjGxZ/ksXldRHuYfy2CjxpIf5fECXiqXWbYpFb+tCMqHKXwN/VN s0WY2SrJE4FteQ8SdDF1466skuV4mZ8+X8DOABhxya8aSqqkeZYD+oDQz1ZExd77 RanRJ0o= Received: (qmail 22793 invoked by alias); 19 Sep 2017 12:56:33 -0000 Mailing-List: contact libc-alpha-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Unsubscribe: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: libc-alpha-owner@sourceware.org Delivered-To: mailing list libc-alpha@sourceware.org Received: (qmail 22784 invoked by uid 89); 19 Sep 2017 12:56:32 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-24.6 required=5.0 tests=AWL, BAYES_00, GIT_PATCH_0, GIT_PATCH_1, GIT_PATCH_2, GIT_PATCH_3, RCVD_IN_DNSWL_NONE, SPF_HELO_PASS, SPF_PASS autolearn=ham version=3.3.2 spammy= X-HELO: EUR03-DB5-obe.outbound.protection.outlook.com Message-ID: <59C113F7.1030404@arm.com> Date: Tue, 19 Sep 2017 13:56:23 +0100 From: Szabolcs Nagy User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:31.0) Gecko/20100101 Thunderbird/31.8.0 MIME-Version: 1.0 To: GNU C Library CC: nd@arm.com Subject: [PATCH 3/9 v2] Do not wrap expf and exp2f References: <59C1123F.9080003@arm.com> In-Reply-To: <59C1123F.9080003@arm.com> X-ClientProxiedBy: DB6PR06CA0012.eurprd06.prod.outlook.com (2603:10a6:6:1::25) To VI1PR0802MB2494.eurprd08.prod.outlook.com (2603:10a6:800:b6::22) X-MS-PublicTrafficType: Email X-MS-Office365-Filtering-Correlation-Id: 0ae05b96-c7f6-4185-86a9-08d4ff5dd6c9 X-MS-Office365-Filtering-HT: Tenant X-Microsoft-Antispam: UriScan:; BCL:0; PCL:0; RULEID:(300000500095)(300135000095)(300000501095)(300135300095)(22001)(300000502095)(300135100095)(2017030254152)(300000503095)(300135400095)(48565401081)(49563074)(201703131423075)(201703031133081)(201702281549075)(300000504095)(300135200095)(300000505095)(300135600095)(300000506095)(300135500095); SRVR:VI1PR0802MB2494; X-Microsoft-Exchange-Diagnostics: 1; VI1PR0802MB2494; 3:BN7pVsNQ2gIgTsIQmn/UUNk2zlRTzMC9+ygDKA8IeNffUjv4xlXgQ5XibcaxSz/2M6AZa6XaeTBJE9mex2PgFHPwUCN0WW7Yvfd0HXO65XxwGbexzZYFVlLuO1Csx3Dqmdzqi8sfpYjokewY6CAIZhyOBHgs+GpQ7WiVA+RUF+8NDWmhdKzvnz+6VgGgCoB+ipJyF9+G/YiqWn8lUOs2CZyNaXRw6hI/Lv0qZvu3fFA77tg87oGa74REd3Ez/avR; 25:V1a7F1MOofwvyWyN/ew35T8vjuALpuuDv5q+f4IpZwfslj05a7UmmtV49dZaOpXpruE5o683HfOXaBMacEhEBTOKcApFubhSwXSbV06YbuRNgo5cLUGmoZzWE1Lkd/sqtFfNs9nYicsBkw2eDrUjKxux17DLyy5yB6Hk29RGgTFD2ADVlaQQVwl+9x1hZg19jgVfr/njAXDx6EY5d3Tg2dBd0TST7hj4pPgYXbIXuk8CSyPyEt5tfM/fYT+LfH+UNc25MojQE5Cn37AlEUOCSkvU6NYAXa5bW5I/CNKoYuGM2l0S192UPPTd9XqaaRleQqw/U4F/u/ofbrRKSdOLrA==; 31:dOXohLNChHrLXG9+rJaW5vGX/o4sS9mcVVxa53PxJYHRjvTnoM2uYB40xCe800HSknu8RIbhk2d/0stn9xHST8/ijbMrcHrR4rwKYZz5HlEv14QoE52JMA0I9/hXPQT/cUsVumQdxVG6qxSOlUw8Q2LwTgr6GI+lSxDy4NLnk+vs6c5FSfipjzDyBCAH+QNT+/gW95ZlPosuXzQs48NiPq5GXDllYQoDMQXSBfLdYEs= X-MS-TrafficTypeDiagnostic: VI1PR0802MB2494: Authentication-Results: spf=none (sender IP is ) smtp.mailfrom=Szabolcs.Nagy@arm.com; NoDisclaimer: True X-Microsoft-Exchange-Diagnostics: 1; VI1PR0802MB2494; 20:hxzQZSeMWpIJzEnNHsEW9yzsxmfjDsNcYB8MdfzK4z0h+Pm01Lvp3jgTYK6jhTuBT662BVUeaJu/A7c/yr59j3WM8cVgmge931hdvC26q2/a6Cmms0Cc7MMkrrKfVsYh8tr3RuGqAKF5pgUwMMaZAKftBIKFNCWpOz883Z1ufDY=; 4:J+HjONEtwNaC0Mdi2XIfeMkfZRRNwf6/hFwtfgbrM8MAYO068ouL8yLTiqXdm8H4WArOg9pD5SWAV4UpLRO7YjPsrMbfWLwZh71APCBVqRGNsmI5ax4B4wWMT1fObw2GNWxkuKOHDklj44zUNtkA8H5AjbZvdo0J2azVu96FTyIn6FkC2UG5kdZ0qhu26fc83OIpEXTtLysU04BijGOYgg37ESC0oJ6ck1w4Wl0+GewiXXlDPVBRqde+ZEa6qbbVVo/bxpwv229Q0Ei++GIhU/MET4XWR/GQsEWkF3uZ5HM= X-Exchange-Antispam-Report-Test: UriScan:(180628864354917); X-Microsoft-Antispam-PRVS: X-Exchange-Antispam-Report-CFA-Test: BCL:0; PCL:0; RULEID:(100000700101)(100105000095)(100000701101)(100105300095)(100000702101)(100105100095)(102415395)(6040450)(2401047)(5005006)(8121501046)(3002001)(93006095)(93001095)(10201501046)(100000703101)(100105400095)(6055026)(6041248)(20161123564025)(20161123555025)(20161123558100)(20161123562025)(20161123560025)(201703131423075)(201702281528075)(201703061421075)(201703061406153)(6072148)(201708071742011)(100000704101)(100105200095)(100000705101)(100105500095); SRVR:VI1PR0802MB2494; BCL:0; PCL:0; RULEID:(100000800101)(100110000095)(100000801101)(100110300095)(100000802101)(100110100095)(100000803101)(100110400095)(100000804101)(100110200095)(100000805101)(100110500095); SRVR:VI1PR0802MB2494; X-Forefront-PRVS: 04359FAD81 X-Forefront-Antispam-Report: SFV:NSPM; SFS:(10009020)(4630300001)(7370300001)(6009001)(6049001)(376002)(39860400002)(346002)(377424004)(199003)(189002)(72206003)(5660300001)(5890100001)(8676002)(68736007)(86362001)(21490400003)(77096006)(59896002)(106356001)(6486002)(6666003)(36756003)(80316001)(189998001)(7350300001)(87266999)(76176999)(65816999)(54356999)(81166006)(81156014)(50986999)(3846002)(101416001)(6116002)(53936002)(16586007)(305945005)(105586002)(16576012)(7736002)(2476003)(25786009)(2906002)(478600001)(84326002)(16526017)(83506001)(4610100001)(6916009)(316002)(568964002)(2950100002)(270700001)(64126003)(33656002)(66066001)(97736004)(58126008)(65956001)(65806001)(4326008); DIR:OUT; SFP:1101; SCL:1; SRVR:VI1PR0802MB2494; H:[10.2.206.69]; FPR:; SPF:None; PTR:InfoNoRecords; A:1; MX:1; LANG:en; Received-SPF: None (protection.outlook.com: arm.com does not designate permitted sender hosts) X-Microsoft-Exchange-Diagnostics: 1; VI1PR0802MB2494; 23:sxR6wV3dFjyIxyI+D2tV21ki5oc+J7GR7bbK+wlWX5khsIhhIQpyFPxjfiKTwiIbymkB3cgBhml8d9wcFLhq5HgqqzsABMnAiDr+mbODzmY/+06Cxhpp/rrGrzuhUImVTI20YXuHW1qn12jima6LE4f9BrpUrKFEspBouyhSzgejfnjfBxcCYmBLK5R4pefbCYJqmFAb+AoikcNu43zN/YOO1ko+iL/NaP3+OJVI2j1SjpF9PIXvQPrBJ5GkKFd4sFQ2J9jfS+ghIaaIbxuhiwh421KapC7Cwu2pgBEc2SNCRR0LMq213LNKVrjMH20kZqWXe4nkswPIKoflhni5CH0a0Or7YnXFww686ImjXeG8lV2SZXVdEpNZg9yzbCyYSodrViLjP3DrSTMzcRi6lUm3gfWsUuHqHWfoih2WkdK1+4NgJQrjWArFFvCz9eSgiDpeve/PvcZl+NPQWrfIfHer+pU69bVk0uOSQPg8sv0EulXqfJkJ5yhkMfQKHMHP3r/A9luJzAsTO9uqGOQbawXkIvMJjhh/Y/jxGPEgVWn1gXzUBWBFUjg8g4US6+RhGFIzhNebzZ091hhCxDWq72gBdrANsMUpH5yE+7Amkm0kdQSqfx3rBKZosCT1PZpkPytj8h180EVsRU/46oxGHkdHQLT53SgGDInjDFOgv5ZDRBScnb7fyrvvQ0Y0X0A3AR1kbaNJs80+kA2DyKqBHvHl7JHLETYw57ZPP2/YoXVNolZfN2KHxf1HFRbScVm4qULcDgvG9qnhRj6U3MijVKLcEDyd1WpXS/Wl8YRokdOyo3BHWBkb9jpYzsajHmWH6JDBndPTEQ1BJiCybvBpnPFL0/cDEt6fePa7ydcIzm9Quc0yc8JswcNulQP/QVgqckE7G1IuI30gfh6TJu6RIWDbF0nDUHlSzlXFPBny4+ClSEREobo/vEF5D6mQbypQyP9EJaMW5JnoRKOR9cpGwV5QRx6rf3H+L79kv1VDn4a8Sni17J8goPeVFK6SiQ0tTZY6wrTDWY+fbdCoqYfoAuoTI7wr2TguGXY5Yg3EV0Xppqb/4DNX4bDbTnZAPiiPb4GTg95MzkDnDSNCGNYX3iXR/3Gc87XC0B4s8yascTdIGVzATwPXdSooBSUsvfYV7yWUun9HYkAPEtXxBRjlnefCBWAEQZAMIwmqo7nHTsIjbilFbkxuIrgqWCSrBCn5WzSzXceBdmIjf//irMs4M+WUu0IXWgaLVVaPjG9vnDNzMAHLnlyuIh1qHY1uMtfkjaSmKAueFxqh8r2y/jPZi8cjeDuCUT6ldYL9OvNtxFlEQxSWVd1SFMbhUxoaVVTkyj+pQYoqFWgxY7V/e2OqkahzZZwJThgO3E2kJtsqZDAOAhAogQHiuYio12XCk/ZHv/gEeq3twH5/5WTnWavI8zNgH0S54tS4owStbhKlTSVVwf+K7hvR8iskF9wBDtcjsym1hmkd3Ig/IZd4aTaL8g== X-Microsoft-Exchange-Diagnostics: 1; VI1PR0802MB2494; 6:7DskRqcHDxca1BAfQ/qxUFH103yiKvMl3emdP+SmnpKxj0Fg1PdCS9KuO/mEZ8sEguBgWUUlRNIfgeIqMZQKaKhDNpYVIIwYshIpPjOVCIbZGwzaQdlX1BGSHobzjfghZaZYJ320G5YD/E5Epjf9MZ3K2UxUdsvVjJA8pirTjo1UW4JWzH2uvro0p1fU+HjPk8h1tevVjr8ZAv6aCtNNHkWrVe7+slo7Zh0qMhnqJumCIKVtUePhP1fCACXcnxAaNXj9jT12dkrJ+QBW09ZX4FmwquvvlI3GsSOyNiuKGowhklep48Ci6RaA6SITrfzxZpVJx9cPx3UzNWsVk+J6jA==; 5:w4ggKz6jIPAsEnFLE7rxfZ3n+eDWDgYkJz9sBm8Oz9pW8I8tnf4w1XNbBqOb0MUP1C6zGCWoWAyTNHDaNdH0B6bDS6NI581/93Uw1tkCjRJCMwuXBU8t4qGB7xM6rpFcMFLRtQI0FrtD5srEsL4Seg==; 24:h0faWY9f1cWK8Px7aU9aD/qiTVxx72kNciU981TPVV9nyYfmveDwwHynqvOGpwPbQEh6+R2fVXXwv985PNa9ZutdxWlgObM9JJWERyU9nII=; 7:2cdWs8mZkNbzg4imoRGHyxIpvA2ysIw4nMFhyVuTeYSb2LuYzTePlsN9H+RRvi00Rs5pQeXEopKu7L7q3f4Lc9/z1zUqubzcnvvusgDT0mrmA9nZO9g5kdvQdKTiEhOGguQGi3YSP7nBrIz2qmC5SjG54KR56OI0gopRX/W60UTWetir8qWOMzxC7bf9ue9spWM7OtXdESqHczY4Fl/msgO6txR5ih2qwCH2ORhbtTc= SpamDiagnosticOutput: 1:99 SpamDiagnosticMetadata: NSPM X-OriginatorOrg: arm.com X-MS-Exchange-CrossTenant-OriginalArrivalTime: 19 Sep 2017 12:56:26.2749 (UTC) X-MS-Exchange-CrossTenant-FromEntityHeader: Hosted X-MS-Exchange-Transport-CrossTenantHeadersStamped: VI1PR0802MB2494 The new generic expf and exp2f code don't need wrappers any more, they set errno inline, so only use the wrappers on targets that need it. (The top level wrapper code is included into sysdeps dirs where the wrapper is needed and empty w_exp*f.c is used to suppress the wrapper). A powerpc64 expf implementation includes the expf c code directly which needed some changes. v2: - Fix sysdeps/powerpc64/fpu/multiarch/e_expf-ppc64.c - Use "sysdeps/../math/" hack in sysdeps/*/w_expf.c - Add sysdeps/m68k/m680x0/fpu/w_exp2f.c 2017-09-19 Szabolcs Nagy * sysdeps/ieee754/flt-32/e_exp2f.c (__exp2f): Define without wrapper. * sysdeps/ieee754/flt-32/e_expf.c (__expf): Likewise * sysdeps/ieee754/flt-32/w_exp2f.c: New file. * sysdeps/ieee754/flt-32/w_expf.c: New file. * sysdeps/powerpc/powerpc64/fpu/multiarch/e_expf-ppc64.c: Update for the new expf code. * sysdeps/powerpc/powerpc64/fpu/multiarch/w_expf.c: New file. * sysdeps/m68k/m680x0/fpu/w_exp2f.c: New file. * sysdeps/m68k/m680x0/fpu/w_expf.c: New file. * sysdeps/i386/fpu/w_exp2f.c: New file. * sysdeps/i386/fpu/w_expf.c: New file. * sysdeps/i386/i686/fpu/multiarch/w_expf.c: New file. * sysdeps/x86_64/fpu/w_expf.c: New file. diff --git a/sysdeps/i386/fpu/w_exp2f.c b/sysdeps/i386/fpu/w_exp2f.c new file mode 100644 index 0000000000..583065d12a --- /dev/null +++ b/sysdeps/i386/fpu/w_exp2f.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/i386/fpu/w_expf.c b/sysdeps/i386/fpu/w_expf.c new file mode 100644 index 0000000000..b5fe164520 --- /dev/null +++ b/sysdeps/i386/fpu/w_expf.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/i386/i686/fpu/multiarch/w_expf.c b/sysdeps/i386/i686/fpu/multiarch/w_expf.c new file mode 100644 index 0000000000..b5fe164520 --- /dev/null +++ b/sysdeps/i386/i686/fpu/multiarch/w_expf.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/ieee754/flt-32/e_exp2f.c b/sysdeps/ieee754/flt-32/e_exp2f.c index be53e5c699..7a3c5fdcff 100644 --- a/sysdeps/ieee754/flt-32/e_exp2f.c +++ b/sysdeps/ieee754/flt-32/e_exp2f.c @@ -18,6 +18,7 @@ #include #include +#include #include "math_config.h" /* @@ -42,7 +43,7 @@ top12 (float x) } float -__ieee754_exp2f (float x) +__exp2f (float x) { uint32_t abstop; uint64_t ki, t; @@ -85,4 +86,6 @@ __ieee754_exp2f (float x) y = y * s; return (float) y; } -strong_alias (__ieee754_exp2f, __exp2f_finite) +strong_alias (__exp2f, __ieee754_exp2f) +strong_alias (__exp2f, __exp2f_finite) +versioned_symbol (libm, __exp2f, exp2f, GLIBC_2_27); diff --git a/sysdeps/ieee754/flt-32/e_expf.c b/sysdeps/ieee754/flt-32/e_expf.c index b8244a82cf..8a38c2a0c6 100644 --- a/sysdeps/ieee754/flt-32/e_expf.c +++ b/sysdeps/ieee754/flt-32/e_expf.c @@ -18,6 +18,7 @@ #include #include +#include #include "math_config.h" /* @@ -42,7 +43,7 @@ top12 (float x) } float -__ieee754_expf (float x) +__expf (float x) { uint32_t abstop; uint64_t ki, t; @@ -99,4 +100,7 @@ __ieee754_expf (float x) y = y * s; return (float) y; } -strong_alias (__ieee754_expf, __expf_finite) +hidden_def (__expf) +strong_alias (__expf, __ieee754_expf) +strong_alias (__expf, __expf_finite) +versioned_symbol (libm, __expf, expf, GLIBC_2_27); diff --git a/sysdeps/ieee754/flt-32/w_exp2f.c b/sysdeps/ieee754/flt-32/w_exp2f.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/ieee754/flt-32/w_exp2f.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ieee754/flt-32/w_expf.c b/sysdeps/ieee754/flt-32/w_expf.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/ieee754/flt-32/w_expf.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/m68k/m680x0/fpu/w_exp2f.c b/sysdeps/m68k/m680x0/fpu/w_exp2f.c new file mode 100644 index 0000000000..583065d12a --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/w_exp2f.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/m68k/m680x0/fpu/w_expf.c b/sysdeps/m68k/m680x0/fpu/w_expf.c new file mode 100644 index 0000000000..b5fe164520 --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/w_expf.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/powerpc/powerpc64/fpu/multiarch/e_expf-ppc64.c b/sysdeps/powerpc/powerpc64/fpu/multiarch/e_expf-ppc64.c index b236290ea2..ee5667c0d5 100644 --- a/sysdeps/powerpc/powerpc64/fpu/multiarch/e_expf-ppc64.c +++ b/sysdeps/powerpc/powerpc64/fpu/multiarch/e_expf-ppc64.c @@ -16,9 +16,16 @@ License along with the GNU C Library; if not, see . */ +#include #undef strong_alias #define strong_alias(a, b) +#undef hidden_def +#define hidden_def(a) +#undef libm_hidden_proto +#define libm_hidden_proto(a) +#undef versioned_symbol +#define versioned_symbol(a, b, c, d) -#define __ieee754_expf __ieee754_expf_ppc64 +#define __expf __ieee754_expf_ppc64 #include diff --git a/sysdeps/powerpc/powerpc64/fpu/multiarch/w_expf.c b/sysdeps/powerpc/powerpc64/fpu/multiarch/w_expf.c new file mode 100644 index 0000000000..b5fe164520 --- /dev/null +++ b/sysdeps/powerpc/powerpc64/fpu/multiarch/w_expf.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/x86_64/fpu/w_expf.c b/sysdeps/x86_64/fpu/w_expf.c new file mode 100644 index 0000000000..b5fe164520 --- /dev/null +++ b/sysdeps/x86_64/fpu/w_expf.c @@ -0,0 +1 @@ +#include From patchwork Tue Sep 19 12:58:19 2017 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit X-Patchwork-Submitter: Szabolcs Nagy X-Patchwork-Id: 815511 Return-Path: X-Original-To: incoming@patchwork.ozlabs.org Delivered-To: patchwork-incoming@bilbo.ozlabs.org Authentication-Results: ozlabs.org; spf=pass (mailfrom) smtp.mailfrom=sourceware.org (client-ip=209.132.180.131; helo=sourceware.org; envelope-from=libc-alpha-return-84725-incoming=patchwork.ozlabs.org@sourceware.org; receiver=) Authentication-Results: ozlabs.org; dkim=pass (1024-bit key; secure) header.d=sourceware.org header.i=@sourceware.org header.b="vr7J6o/D"; dkim-atps=neutral Received: from sourceware.org (server1.sourceware.org [209.132.180.131]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (No client certificate requested) by ozlabs.org (Postfix) with ESMTPS id 3xxNGg5qhxz9s7g for ; Tue, 19 Sep 2017 22:58:35 +1000 (AEST) DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; q=dns; s=default; b= tk+62UUeBEw4mqJlAw2KAOIyitTy7YxnZLjAvwfpN8twxYv4I8/SkGyTIskrS6/i oq3lZJW4KIGNt8FJOPlRXjLp/nLPWAnL3SmvhX1zwR7Ub4Ioy1n9FW7Oq+J+P7zg FBMXI48NnfKAHrKSRH+zILbeVFjtAKD1PjnU8GwFefQ= DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; s=default; bh=1enC k1ytJYw+MqrUHadjvlAk/lU=; b=vr7J6o/Dw5ZesmdLdJsxf8obA5sDhq9+BrUs oEN9Of8HFGES0JYUZAL1ZuJMFhoLvVsgKk5LgJpLO/PiWGwWkeaTISaV7IzNAMqS HqnAszv0cjc5Cm36VS7dpFxRKJE9WOi3cNCHy+Jf3JjCwhQv5kuFjmNVDaMd5FkP EBJHBoY= Received: (qmail 25003 invoked by alias); 19 Sep 2017 12:58:29 -0000 Mailing-List: contact libc-alpha-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Unsubscribe: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: libc-alpha-owner@sourceware.org Delivered-To: mailing list libc-alpha@sourceware.org Received: (qmail 24978 invoked by uid 89); 19 Sep 2017 12:58:27 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-23.9 required=5.0 tests=AWL, BAYES_00, GIT_PATCH_0, GIT_PATCH_1, GIT_PATCH_2, GIT_PATCH_3, KAM_ASCII_DIVIDERS, RCVD_IN_DNSWL_NONE, SPF_HELO_PASS, SPF_PASS autolearn=ham version=3.3.2 spammy= X-HELO: EUR01-HE1-obe.outbound.protection.outlook.com Authentication-Results: spf=none (sender IP is ) smtp.mailfrom=Szabolcs.Nagy@arm.com; Message-ID: <59C1146B.5090606@arm.com> Date: Tue, 19 Sep 2017 13:58:19 +0100 From: Szabolcs Nagy User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:31.0) Gecko/20100101 Thunderbird/31.8.0 MIME-Version: 1.0 To: GNU C Library CC: nd@arm.com Subject: [PATCH 4/9] New generic logf References: <59C1123F.9080003@arm.com> In-Reply-To: <59C1123F.9080003@arm.com> X-ClientProxiedBy: DB6PR06CA0022.eurprd06.prod.outlook.com (2603:10a6:6:1::35) To DB6PR0802MB2488.eurprd08.prod.outlook.com (2603:10a6:4:a0::23) X-MS-PublicTrafficType: Email X-MS-Office365-Filtering-Correlation-Id: 14a721ce-fdac-4d41-e5f2-08d4ff5e1ae4 X-MS-Office365-Filtering-HT: Tenant X-Microsoft-Antispam: UriScan:; BCL:0; PCL:0; RULEID:(300000500095)(300135000095)(300000501095)(300135300095)(22001)(300000502095)(300135100095)(2017030254152)(48565401081)(300000503095)(300135400095)(49563074)(201703131423075)(201703031133081)(201702281549075)(300000504095)(300135200095)(300000505095)(300135600095)(300000506095)(300135500095); SRVR:DB6PR0802MB2488; X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 3:Otpc3RTrK31XnvTN71MdhzIVK22z/+C3VOer33AyGdOKy3fIK1e4AxXFvHHzCWngo/dh4NIr8RJ9i7yLda2dUrvG3lcmpUjL4+O/Fy33VfThmxFOFl1Sgnce8Sp1zayKiKQUg9E2F5u8DgTVCJPKcXLb1zYjwzgxajSBoGD93okYy/eIQJhSLt42izDIXGfg6j0afu9xy0q0ObzysLewTAFGCEy6uvcmLpzyHuNQhsRlnD0hjriChpSnMd/W4nGw; 25:Saei9apP7FxfnwyoobdBvd0H1WSNTWxdKVgZ4OS7P2GAAwnYknaORj3TL0GpREQtf4Kx2aGehFtGh1X3Gpv54zPxh1l8h5/57ReeOB/PoAbPoB/2xYqJbZ19QSJDS1lvvjyeGVjVtrEaELOM4rCVJ1Zpt/1z65JsAsGQTYpLxOlcbSghRgEhSsGoCH5Z7YP19QegJwnN7BCh0Fty56IS7AphZtazf5KjE3DS2XFUawRRDjowEknto36uL1okaJzjrt+6+5jEdX6OXWsHH0+n4+4gr3EJen+pkWGo+4MYaeFO2/71p5V0D7BB6tz6O9Kzel0r8n8BSS9u4SO1FYrFvA==; 31:JpFyUI3K8kFeQ6nnYSB11VLfYiWxO6ZRIcwBcYaMu1gZLjCgfbid0uyWM6UUqs9g3C3vEdX+LYY9nIOFMtJ/5a8P9bp4GX6XcdB118Ox2SkMlcZxlfWSh8loY7+SF92Q45/MdGRSJ+HulKP3nxFFtDfIVsXWFRxjCNji0SdY9T4PoE09p15t8KRYENRjgrY6y5UZgi/TUAwXrOrKFSSkc0SI/eAv0huq5vRgFHq1Z1U= X-MS-TrafficTypeDiagnostic: DB6PR0802MB2488: NoDisclaimer: True X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 20:uL2/lzHr3sIDAhiMFJuX51uqFs/aLSLwxe6dGYjSLl9+nlpzC2G6u2fEuZgogSg8RemhDIfGNvW0yfkmSDmPTvuz4PXUWWbOir+iNQEcQfE++MlJ/S3qpdeOEdmpPBROwQ1z2LfAR00lJQa3AN58RN0iIT11V7WiN/uu7RLRM90=; 4:W60WPCUOULLsB675wzxnYD7cPUOOiOYG+9QhtMJFLKzG24N7lKZsVLV8sZQRMjb0pvALIQ70smjAJ+jLPGhjLHWiwCe4GuSe7FK+T9LjU1QbXYAIuZfugHRhhahbMVGjsVQulJAKLgJhIRee3yCugcI4B7SGgPTo5kWmskXu60FG+ul4S+LY6hSUpmC2/fnxbkyN3/ZJO7BNRfDfgHw3IfKeqq9mNojBk3bUzuM2pESUkh3u4vEromkhbJ8b2HyVBTljt/VAqRdrx8MsCqlIIMXjkM0Qsf7xh1nTwwBp/iU= X-Exchange-Antispam-Report-Test: UriScan:(180628864354917); X-Microsoft-Antispam-PRVS: X-Exchange-Antispam-Report-CFA-Test: BCL:0; PCL:0; RULEID:(100000700101)(100105000095)(100000701101)(100105300095)(100000702101)(100105100095)(102415395)(6040450)(2401047)(8121501046)(5005006)(100000703101)(100105400095)(10201501046)(3002001)(93006095)(93001095)(6055026)(6041248)(20161123562025)(20161123564025)(20161123555025)(201703131423075)(201702281528075)(201703061421075)(201703061406153)(20161123560025)(20161123558100)(6072148)(201708071742011)(100000704101)(100105200095)(100000705101)(100105500095); SRVR:DB6PR0802MB2488; BCL:0; PCL:0; RULEID:(100000800101)(100110000095)(100000801101)(100110300095)(100000802101)(100110100095)(100000803101)(100110400095)(100000804101)(100110200095)(100000805101)(100110500095); SRVR:DB6PR0802MB2488; X-Forefront-PRVS: 04359FAD81 X-Forefront-Antispam-Report: SFV:NSPM; SFS:(10009020)(4630300001)(7370300001)(6009001)(6049001)(346002)(39860400002)(376002)(199003)(189002)(377424004)(83506001)(25786009)(4326008)(101416001)(59896002)(50986999)(478600001)(105586002)(106356001)(2906002)(87266999)(80316001)(189998001)(65816999)(54356999)(76176999)(84326002)(6116002)(3846002)(72206003)(2476003)(270700001)(568964002)(68736007)(4610100001)(5890100001)(6486002)(77096006)(7350300001)(2950100002)(36756003)(5660300001)(7736002)(305945005)(33656002)(66066001)(53936002)(16526017)(16576012)(65956001)(316002)(65806001)(16586007)(58126008)(64126003)(97736004)(8676002)(6916009)(21490400003)(81166006)(81156014)(86362001)(41533002); DIR:OUT; SFP:1101; SCL:1; SRVR:DB6PR0802MB2488; H:[10.2.206.69]; FPR:; SPF:None; PTR:InfoNoRecords; A:1; MX:1; LANG:en; Received-SPF: None (protection.outlook.com: arm.com does not designate permitted sender hosts) X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 23:mfVvUF8BFJ9FP1dAxNStDOAxB4DyUN1hV6/kasYZ9aR0jto4R+7u9qI/7ZgUW1jNkrEa37cbdvr9PGOQKaxJ0K0tBYuVFytHA5f0UOfVyd4t3EqybioLFYTISnY2s3R6xtiBFnZIDj7LR0j3B5uWefHoPOKNGo3qN/rh+FOdQlWApOZXyne2yhycwlBmdXhXp4NamZ1ICsiqycfxUe3NvtaFLEsrnPHXYZ/aU+e2zwitKE1SjfBqfTRmbfX3iJFQEdELea2ssvhm+GTG7uWWRavcDd9MsT4k/dqATd+C/k6nlzR2Pd53zohZF2IeU48etcpTg4UBWLJQgjfErbRc9pe/L/CCdQydK8aPQaiYQUvKbA+aWVU0zQAfu6wZiXORNGTYlDSx+qiSZwRl6jBACrphxPxnUqLfjWanHXIRXlge7iuhFNG38eCij0+X877Bmv097hG6zzwvlmEGE8eL6pYkU3clsBSrk/F6oe7oEbn0lBzBZ3Hu4LkrRr/st8AwxLPGhRGCOMYX81dZ7WgKTJWkTcMHLZW6lNh8ipMbpZxZpJyiyexawIAyNhYxlNP//Lt27DwIZC8DBRxSyzXzMlDMn7zgKUpPp+rXorPy0WgnHdAkhM73iz0FF1IFfdyLOcZJ29qd4E197Zz23A1CxR9wsnRz2GsnMjuT/8YiPlQCM/K7T4KAsmU8FNGkXkaH4WQLjOpJBDvwGjdC6y70T3roKCQrWR6U5rzKhy9WhjdmBN7Hm9KAFX8NTiLiRSKO4g8is4zAgacRorrN+iDpiZo377eCZczRVE5mSiNGg+KHQk2xPHWU2Gr0ZxAKJ3mNitBDgLb7hSyFEkilm/n1veFm5KqhMVWx6+kQD9EBp7TSXFXGT6Z3AkYcJxf2DKLQBsbASoA46yJ7sX9GhegnJ/0EwJ5cLqyXe2etbzmSJeWANR18iw8M8aaziSJOIjgnYaQVDJ07TzuvB6RUpiImEDmSegHt9xRZ3NXAMSsFAh2mkrZipZuNkd8YJoEftcxp1o8J/R5iK1MYuP1PRNgQk7scs+Ez8k90YupKOrbTtnHH8aFyZXvwtDSffOtf4TZWLa2qPG4ElFx8teBIdiGX9sxqzMQic5YI6nViKqZ+oVPgj+sCVL1vEtu2xeldmzqnRF2QG35kVBe8rGlMGUQCyXg1ClpW0hz/KmnTjC+5yFzm0iVQWqAoPyytwSKLgcDD/Oy4xkfTw6SxS1FPBUy/VfQjphFR+8RU2E3K7tX+jjFneLxdcWD0zMFuxTZ6oDsx7SbgS864VODUrfgZXGcm4Mv9GqkUscUnPkYd9sXUUZfQ4MUKbds/qDCF8sl5QNjeXvkftweqvfHoibj/Bd9rRP3hL3aYJOBgoaiPqq2SZsVSkPxygF67Y1uJuCmqKEnBYaRpH3N/5sWXAZ1lXjRn62oS7Yv+p2tWroYTtZ3IuHoLn7bviaiXYZPB3OxbEQOgBNbsOW0rJRDZBD0bAELCsA== X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 6:GWFE9kxT15hwbL1p0IAw3EjosC5ay9JrfxN3jTo94Y878uvmd8pw5LK4PW9IWHZChVAg0cF9sHp7crZncfmXKFQQJ7H7/nmxo4xI+MnSfIpsgiU6cy23OX8npadYHuikc8DO/WJcApllmq7sq+30uI30tT2CtOFs9OkwzcNcVWEbHPpG331XvTmdMNySxvSz69wRCXfhILAet2ErpJZyIGZwJp4rVaW08y52lIsW3gVG6lLbVPIS8oXpxLDAzohQBI2d9jS18gl9Igwm5siFQf6bf28cqc9w/cA7b6j0Tsz7xNigKn0CaFopfnsTZQB5rU2uJAH90k/8dAOzdF1UWg==; 5:kb1qYuirhMl5q+278n5g/qxchwiZm4KL+B7ehzUKL/KmnJDAV1rtJ7ZUK4xbsOSjn0WmX3ZqMna/kR7SQqJ6ep8187spM4crKSfkKHdIVsgxrqbzVniCRbeILS6+HTecQ3SYO8zhCKG9EYVY2NF5bg==; 24:BMXIBlHY0cOkSNeDgQPCsZA2SSRr0J29ICR1yd9q9HkviYV+vH9OUupkcZ3aqKo4OFEHH82Ju4pYtbR/hamMdWq0JHB5w2XDpp0+NsMR/Q8=; 7:yPodgbJ5BTvRysQz7eAQvbIZlPezJc258pozMedz3T7qNpGSSZZ82XM2OUI+4CNE0Dr6KBCAm6fQ1tXCrrIKLuDqOOaCEasQXOdKjWBQAli/SCtoN87mhH7oI9saZRkU9tC2XQ07jn1OxNDm7jECCI6RrK1ZofJY+GeoLP8s/cvJ4yJPMt62UVo8fzdzlAgOuMPVTyXX3IHmH8CVBm2pDAEDHHpem26a9fLg+kYhgBo= SpamDiagnosticOutput: 1:99 SpamDiagnosticMetadata: NSPM X-OriginatorOrg: arm.com X-MS-Exchange-CrossTenant-OriginalArrivalTime: 19 Sep 2017 12:58:20.8049 (UTC) X-MS-Exchange-CrossTenant-FromEntityHeader: Hosted X-MS-Exchange-Transport-CrossTenantHeadersStamped: DB6PR0802MB2488 with the logf benchmark (without wrapper, aarch64): reciprocal-throughput: 2.2x faster latency: 1.9x faster old worst case error: 0.89 ulp new worst case error: 0.82 ulp aarch64 .text size: -356 bytes aarch64 .rodata size: +248 bytes Uses double precision arithmetics and a lookup table to allow smaller polynomial and avoid the use of division. Data is in a separate translation unit with fixed layout to prevent the compiler generating suboptimal literal access. Errors are handled inline according to POSIX rules, but this patch keeps the wrapper with SVID compatible error handling. Needs libm-test-ulps adjustment for clogf in non-nearest rounding mode. 2017-09-19 Szabolcs Nagy * math/Makefile (type-float-routines): Add e_logf_data. * sysdeps/ieee754/flt-32/e_logf.c: New implementation. * sysdeps/ieee754/flt-32/e_logf_data.c: New file. * sysdeps/ieee754/flt-32/math_config.h (__logf_data): Define. (LOGF_TABLE_BITS, LOGF_POLY_ORDER): Define. * sysdeps/i386/fpu/e_logf_data.c: New file. * sysdeps/ia64/fpu/e_logf_data.c: New file. * sysdeps/m68k/m680x0/fpu/e_logf_data.c: New file. diff --git a/math/Makefile b/math/Makefile index 04586156f8..919fec13ef 100644 --- a/math/Makefile +++ b/math/Makefile @@ -115,7 +115,7 @@ type-double-routines := branred doasin dosincos halfulp mpa mpatan2 \ # float support type-float-suffix := f -type-float-routines := k_rem_pio2f math_errf e_exp2f_data +type-float-routines := k_rem_pio2f math_errf e_exp2f_data e_logf_data # _Float128 support type-float128-suffix := f128 diff --git a/sysdeps/i386/fpu/e_logf_data.c b/sysdeps/i386/fpu/e_logf_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/i386/fpu/e_logf_data.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ia64/fpu/e_logf_data.c b/sysdeps/ia64/fpu/e_logf_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/ia64/fpu/e_logf_data.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ieee754/flt-32/e_logf.c b/sysdeps/ieee754/flt-32/e_logf.c index cf75e11781..b8d262441f 100644 --- a/sysdeps/ieee754/flt-32/e_logf.c +++ b/sysdeps/ieee754/flt-32/e_logf.c @@ -1,85 +1,87 @@ -/* e_logf.c -- float version of e_log.c. - * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com. - */ +/* Single-precision log function. + Copyright (C) 2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. -/* - * ==================================================== - * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. - * - * Developed at SunPro, a Sun Microsystems, Inc. business. - * Permission to use, copy, modify, and distribute this - * software is freely granted, provided that this notice - * is preserved. - * ==================================================== - */ + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ #include -#include +#include +#include "math_config.h" -static const float -ln2_hi = 6.9313812256e-01, /* 0x3f317180 */ -ln2_lo = 9.0580006145e-06, /* 0x3717f7d1 */ -two25 = 3.355443200e+07, /* 0x4c000000 */ -Lg1 = 6.6666668653e-01, /* 3F2AAAAB */ -Lg2 = 4.0000000596e-01, /* 3ECCCCCD */ -Lg3 = 2.8571429849e-01, /* 3E924925 */ -Lg4 = 2.2222198546e-01, /* 3E638E29 */ -Lg5 = 1.8183572590e-01, /* 3E3A3325 */ -Lg6 = 1.5313838422e-01, /* 3E1CD04F */ -Lg7 = 1.4798198640e-01; /* 3E178897 */ +/* +LOGF_TABLE_BITS = 4 +LOGF_POLY_ORDER = 4 -static const float zero = 0.0; +ULP error: 0.818 (nearest rounding.) +Relative error: 1.957 * 2^-26 (before rounding.) +*/ + +#define T __logf_data.tab +#define A __logf_data.poly +#define Ln2 __logf_data.ln2 +#define N (1 << LOGF_TABLE_BITS) +#define OFF 0x3f330000 float -__ieee754_logf(float x) +__ieee754_logf (float x) { - float hfsq,f,s,z,R,w,t1,t2,dk; - int32_t k,ix,i,j; + /* double_t for better performance on targets with FLT_EVAL_METHOD==2. */ + double_t z, r, r2, y, y0, invc, logc; + uint32_t ix, iz, tmp; + int k, i; + + ix = asuint (x); +#if WANT_ROUNDING + /* Fix sign of zero with downward rounding when x==1. */ + if (__glibc_unlikely (ix == 0x3f800000)) + return 0; +#endif + if (__glibc_unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) + { + /* x < 0x1p-126 or inf or nan. */ + if (ix * 2 == 0) + return __math_divzerof (1); + if (ix == 0x7f800000) /* log(inf) == inf. */ + return x; + if ((ix & 0x80000000) || ix * 2 >= 0xff000000) + return __math_invalidf (x); + /* x is subnormal, normalize it. */ + ix = asuint (x * 0x1p23f); + ix -= 23 << 23; + } + + /* x = 2^k z; where z is in range [OFF,2*OFF] and exact. + The range is split into N subintervals. + The ith subinterval contains z and c is near its center. */ + tmp = ix - OFF; + i = (tmp >> (23 - LOGF_TABLE_BITS)) % N; + k = (int32_t) tmp >> 23; /* arithmetic shift */ + iz = ix - (tmp & 0x1ff << 23); + invc = T[i].invc; + logc = T[i].logc; + z = (double_t) asfloat (iz); - GET_FLOAT_WORD(ix,x); + /* log(x) = log1p(z/c-1) + log(c) + k*Ln2 */ + r = z * invc - 1; + y0 = logc + (double_t) k * Ln2; - k=0; - if (ix < 0x00800000) { /* x < 2**-126 */ - if (__builtin_expect((ix&0x7fffffff)==0, 0)) - return -two25/zero; /* log(+-0)=-inf */ - if (__builtin_expect(ix<0, 0)) - return (x-x)/(x-x); /* log(-#) = NaN */ - k -= 25; x *= two25; /* subnormal number, scale up x */ - GET_FLOAT_WORD(ix,x); - } - if (__builtin_expect(ix >= 0x7f800000, 0)) return x+x; - k += (ix>>23)-127; - ix &= 0x007fffff; - i = (ix+(0x95f64<<3))&0x800000; - SET_FLOAT_WORD(x,ix|(i^0x3f800000)); /* normalize x or x/2 */ - k += (i>>23); - f = x-(float)1.0; - if((0x007fffff&(15+ix))<16) { /* |f| < 2**-20 */ - if(f==zero) { - if(k==0) return zero; else {dk=(float)k; - return dk*ln2_hi+dk*ln2_lo;} - } - R = f*f*((float)0.5-(float)0.33333333333333333*f); - if(k==0) return f-R; else {dk=(float)k; - return dk*ln2_hi-((R-dk*ln2_lo)-f);} - } - s = f/((float)2.0+f); - dk = (float)k; - z = s*s; - i = ix-(0x6147a<<3); - w = z*z; - j = (0x6b851<<3)-ix; - t1= w*(Lg2+w*(Lg4+w*Lg6)); - t2= z*(Lg1+w*(Lg3+w*(Lg5+w*Lg7))); - i |= j; - R = t2+t1; - if(i>0) { - hfsq=(float)0.5*f*f; - if(k==0) return f-(hfsq-s*(hfsq+R)); else - return dk*ln2_hi-((hfsq-(s*(hfsq+R)+dk*ln2_lo))-f); - } else { - if(k==0) return f-s*(f-R); else - return dk*ln2_hi-((s*(f-R)-dk*ln2_lo)-f); - } + /* Pipelined polynomial evaluation to approximate log1p(r). */ + r2 = r * r; + y = A[1] * r + A[2]; + y = A[0] * r2 + y; + y = y * r2 + (y0 + r); + return (float) y; } strong_alias (__ieee754_logf, __logf_finite) diff --git a/sysdeps/ieee754/flt-32/e_logf_data.c b/sysdeps/ieee754/flt-32/e_logf_data.c new file mode 100644 index 0000000000..ce2a0f8a5c --- /dev/null +++ b/sysdeps/ieee754/flt-32/e_logf_data.c @@ -0,0 +1,44 @@ +/* Data definition for logf. + Copyright (C) 2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#include "math_config.h" + +const struct logf_data __logf_data = { + .tab = { + { 0x1.661ec79f8f3bep+0, -0x1.57bf7808caadep-2 }, + { 0x1.571ed4aaf883dp+0, -0x1.2bef0a7c06ddbp-2 }, + { 0x1.49539f0f010bp+0, -0x1.01eae7f513a67p-2 }, + { 0x1.3c995b0b80385p+0, -0x1.b31d8a68224e9p-3 }, + { 0x1.30d190c8864a5p+0, -0x1.6574f0ac07758p-3 }, + { 0x1.25e227b0b8eap+0, -0x1.1aa2bc79c81p-3 }, + { 0x1.1bb4a4a1a343fp+0, -0x1.a4e76ce8c0e5ep-4 }, + { 0x1.12358f08ae5bap+0, -0x1.1973c5a611cccp-4 }, + { 0x1.0953f419900a7p+0, -0x1.252f438e10c1ep-5 }, + { 0x1p+0, 0x0p+0 }, + { 0x1.e608cfd9a47acp-1, 0x1.aa5aa5df25984p-5 }, + { 0x1.ca4b31f026aap-1, 0x1.c5e53aa362eb4p-4 }, + { 0x1.b2036576afce6p-1, 0x1.526e57720db08p-3 }, + { 0x1.9c2d163a1aa2dp-1, 0x1.bc2860d22477p-3 }, + { 0x1.886e6037841edp-1, 0x1.1058bc8a07ee1p-2 }, + { 0x1.767dcf5534862p-1, 0x1.4043057b6ee09p-2 }, + }, + .ln2 = 0x1.62e42fefa39efp-1, + .poly = { + -0x1.00ea348b88334p-2, 0x1.5575b0be00b6ap-2, -0x1.ffffef20a4123p-2, + } +}; diff --git a/sysdeps/ieee754/flt-32/math_config.h b/sysdeps/ieee754/flt-32/math_config.h index 1759c447c8..e1feafccfe 100644 --- a/sysdeps/ieee754/flt-32/math_config.h +++ b/sysdeps/ieee754/flt-32/math_config.h @@ -111,4 +111,16 @@ extern const struct exp2f_data double poly_scaled[EXP2F_POLY_ORDER]; } __exp2f_data; +#define LOGF_TABLE_BITS 4 +#define LOGF_POLY_ORDER 4 +extern const struct logf_data +{ + struct + { + double invc, logc; + } tab[1 << LOGF_TABLE_BITS]; + double ln2; + double poly[LOGF_POLY_ORDER]; +} __logf_data; + #endif diff --git a/sysdeps/m68k/m680x0/fpu/e_logf_data.c b/sysdeps/m68k/m680x0/fpu/e_logf_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/e_logf_data.c @@ -0,0 +1 @@ +/* Not needed. */ From patchwork Tue Sep 19 13:00:05 2017 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit X-Patchwork-Submitter: Szabolcs Nagy X-Patchwork-Id: 815517 Return-Path: X-Original-To: incoming@patchwork.ozlabs.org Delivered-To: patchwork-incoming@bilbo.ozlabs.org Authentication-Results: ozlabs.org; spf=pass (mailfrom) smtp.mailfrom=sourceware.org (client-ip=209.132.180.131; helo=sourceware.org; envelope-from=libc-alpha-return-84726-incoming=patchwork.ozlabs.org@sourceware.org; receiver=) Authentication-Results: ozlabs.org; dkim=pass (1024-bit key; secure) header.d=sourceware.org header.i=@sourceware.org header.b="QBdkdMkM"; dkim-atps=neutral Received: from sourceware.org (server1.sourceware.org [209.132.180.131]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (No client certificate requested) by ozlabs.org (Postfix) with ESMTPS id 3xxNK828Vbz9s7g for ; Tue, 19 Sep 2017 23:00:44 +1000 (AEST) DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; q=dns; s=default; b= VI8dMO4poYjupmK8tSkFH0gipBQgBGx7J9/gWjTIyIAXGgOJLLXaF6Sr+bou1P7f kn7nRBbmG5X5q8x0uszU747h4U/4fDCsjqhtj+DqgunK9cBvj69jvf4PiwiDpOSa 4wFOSw72liTSOkiUwE3m342lQIOLSdQogWAjD2DgVA8= DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; s=default; bh=bG7o yzegXfdD11FtlJLh+5tIJGs=; b=QBdkdMkMpwxOzCPuQ/KTDMasYThVRPwCaC0x 6guwPkaj50HL4ayLwRn1CQznA1jTTG2mgbCxWXSuWPh0UHeqhX+aeIXOswFwRVTF iJ5C/biVih3zF1sNU9k8qYB8RaHTCDEls7NQsjKPrjvioXHKkuiBfBRMN6Iur/B5 QIBhpiE= Received: (qmail 32415 invoked by alias); 19 Sep 2017 13:00:33 -0000 Mailing-List: contact libc-alpha-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Unsubscribe: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: libc-alpha-owner@sourceware.org Delivered-To: mailing list libc-alpha@sourceware.org Received: (qmail 29663 invoked by uid 89); 19 Sep 2017 13:00:27 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-23.9 required=5.0 tests=AWL, BAYES_00, GIT_PATCH_0, GIT_PATCH_1, GIT_PATCH_2, GIT_PATCH_3, KAM_ASCII_DIVIDERS, RCVD_IN_DNSWL_NONE, SPF_HELO_PASS, SPF_PASS autolearn=ham version=3.3.2 spammy= X-HELO: EUR01-DB5-obe.outbound.protection.outlook.com Authentication-Results: spf=none (sender IP is ) smtp.mailfrom=Szabolcs.Nagy@arm.com; Message-ID: <59C114D5.5090009@arm.com> Date: Tue, 19 Sep 2017 14:00:05 +0100 From: Szabolcs Nagy User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:31.0) Gecko/20100101 Thunderbird/31.8.0 MIME-Version: 1.0 To: GNU C Library CC: nd@arm.com Subject: [PATCH 5/9] New generic log2f References: <59C1123F.9080003@arm.com> In-Reply-To: <59C1123F.9080003@arm.com> X-ClientProxiedBy: DB6P189CA0022.EURP189.PROD.OUTLOOK.COM (2603:10a6:6:2e::35) To VI1PR0802MB2495.eurprd08.prod.outlook.com (2603:10a6:800:b6::23) X-MS-PublicTrafficType: Email X-MS-Office365-Filtering-Correlation-Id: 0a04a308-ecd7-4858-24fd-08d4ff5e5b68 X-MS-Office365-Filtering-HT: Tenant X-Microsoft-Antispam: UriScan:; BCL:0; PCL:0; RULEID:(300000500095)(300135000095)(300000501095)(300135300095)(22001)(300000502095)(300135100095)(2017030254152)(48565401081)(300000503095)(300135400095)(49563074)(201703131423075)(201703031133081)(201702281549075)(300000504095)(300135200095)(300000505095)(300135600095)(300000506095)(300135500095); SRVR:VI1PR0802MB2495; X-Microsoft-Exchange-Diagnostics: 1; VI1PR0802MB2495; 3:2AMBoEqLgrK9W2H0hI4ieVtgv/iudyPUDb+KDLipOO4AQ3hgOdIv3zt7fH+YVUK/Od27tBGlOHfqNBXFvmAH1jLutjTXuItghzKJbrqMcllD7knK/2guA/wV83C53YISRM/su9cneteqIeRd6aja2TcxG68RXQUcsnDUhsHZZtjmDOYkLfO1OJ3njtyU3qL0/NwFoPFcX5V/Ur8r8rY5mMoZ6mUS1lQyftIV+v9EiTTwGfKt9kc77a1KKjtS8SA3; 25:4kna9W+fvAtgMxBdHYSBpqI/IIPWpvM7o22osZVzyHkOEDTmZ5NGbwOqyFzJGqXWLIOn8P/jD3volNq7R4s3vcgtl74aW6519UfsC7RwfSfnRRv46ihjbh3HVLHxDRCta0iB/FfoKEWHInwWe2HzNKr6ucpgpekW4SOHNiecT5vUEa45U4FLmhkAKkUma5jtqkudur7Fd5lcmNTsMymtU/D/zRWl/w9SvxKKpSepjhVeGuqrbhPFjNbSgk3VI4h2q4fjogBbbBROWHiW2sEqqzei2gYKnCg7DuvsOXdvtbq1Y5Om6HeCsbdEg+waW7Q5MzOAeDXyYqI/z4gTeLsmiQ==; 31:aupgmoc37eJ4jNga3HzB/iBQDtJVDGpP1Xt8NigX7uvwzKjsGvGt4vOWlUbTO8u4omSe7nIkP3GHGfYG+H0fOz+WN4B8Y8f4x4+zJxhafwBJ1DIB9eJkxmthOTXVU9uUXr2KXNzJ161m3cfEFrJxqaRq/V6IXmi5Qvw0Z9IbbL1iqkz2lpY6JXrcjMLt0rlDmkHzPSNPDdFdQO3WghvblIbBLr6uJ9EHdU7zS4vVX78= X-MS-TrafficTypeDiagnostic: VI1PR0802MB2495: NoDisclaimer: True X-Microsoft-Exchange-Diagnostics: 1; VI1PR0802MB2495; 20:tHLPcmK/dsJOVKTTQp2gHiAbQEXADToBF6gAhs322ait3qigf3tqZ5HcRft/r3sCWHX/22eMKDzhBDxr6ESeUEmAG1o/DvWM+XrITrdtEpvTX1arJ52cUA0Ou5UqnIDUy4DenzBp3zYX6rB3AVoxgKarhYs+gesDrG393fbGjoA=; 4:s7G2coaIBlVw1z5oGwHy0dn5jmdrapvLXFDsvXakvs0kNUDGUTtsD2ccAFl++s9n6n72qcNdCvNzAL7eLeEXBNXFci7di2iWYAHahn9uLuCD8o4VBEfpB95B30CbOorJMX3pFeRXxOzqroVCkIE1mm5IKw/aIciUkcbpciFtpxHcyihyevawmyzI0dLKvUz9OowVDSaQG/eKjUk9d79aXItoTYgfTdr7GpLINBEoD7xwPJ0h22X0iut0E9S4gnmp7/szOsnPUXvWM8iUFXnokNBj/fwWzRzybYtWmR3wnTY= X-Exchange-Antispam-Report-Test: UriScan:(180628864354917); X-Microsoft-Antispam-PRVS: X-Exchange-Antispam-Report-CFA-Test: BCL:0; PCL:0; RULEID:(100000700101)(100105000095)(100000701101)(100105300095)(100000702101)(100105100095)(102415395)(6040450)(2401047)(5005006)(8121501046)(93006095)(93001095)(10201501046)(3002001)(100000703101)(100105400095)(6055026)(6041248)(20161123555025)(20161123562025)(20161123564025)(20161123560025)(201703131423075)(201702281528075)(201703061421075)(201703061406153)(20161123558100)(6072148)(201708071742011)(100000704101)(100105200095)(100000705101)(100105500095); SRVR:VI1PR0802MB2495; BCL:0; PCL:0; RULEID:(100000800101)(100110000095)(100000801101)(100110300095)(100000802101)(100110100095)(100000803101)(100110400095)(100000804101)(100110200095)(100000805101)(100110500095); SRVR:VI1PR0802MB2495; X-Forefront-PRVS: 04359FAD81 X-Forefront-Antispam-Report: SFV:NSPM; SFS:(10009020)(7370300001)(4630300001)(6009001)(6049001)(39860400002)(346002)(376002)(199003)(377424004)(189002)(84326002)(2906002)(58126008)(6666003)(25786009)(68736007)(59896002)(4610100001)(5660300001)(77096006)(83506001)(4326008)(72206003)(106356001)(86362001)(80316001)(5890100001)(21490400003)(6486002)(2950100002)(76176999)(53936002)(7736002)(50986999)(54356999)(16586007)(65816999)(6916009)(270700001)(305945005)(105586002)(64126003)(478600001)(66066001)(189998001)(65806001)(65956001)(33656002)(8676002)(16526017)(101416001)(2476003)(36756003)(16576012)(81156014)(97736004)(81166006)(3846002)(316002)(568964002)(6116002)(7350300001)(87266999)(5000100001)(41533002); DIR:OUT; SFP:1101; SCL:1; SRVR:VI1PR0802MB2495; H:[10.2.206.69]; FPR:; SPF:None; PTR:InfoNoRecords; A:1; MX:1; LANG:en; Received-SPF: None (protection.outlook.com: arm.com does not designate permitted sender hosts) X-Microsoft-Exchange-Diagnostics: 1; VI1PR0802MB2495; 23:0o0ruk0W8jyoq0dlPufwEaDHEV8a5RJqyZwXmdxObgarfRoqT+s27MyIcdDNF+bm7BqcRSqyVViliZgQT5PiPgLgjOBkfjWBB3HsIObR8866pdcKoq5lVtni8fCu5KLrC2MPAedS12+6zqCQyD0LLqT9xRlG/Z21+eXGHPuE1zHaOikxylUGyycE7SYG5g7Nk7Cf0u5qZyA9+JmhhIasn+dRWXwUz2r32PyMx9XsqU5tVPubLtb4v3k++QODzY1XNxw3MT1tqdKdYTaWiykx0TU6oUUc6WOFkXfr6joCaq6/m1bN5ND0YKADy0KrMgJsP/DCpZJLpb0EXOiyE40Wj+fUdKIKggtpUcb2lf2ADBN7YT9YMERfoC76k+9ZWYvwF9jcJlZ3RVypzuYQ0TaOg70+apG7M1zrX6T2/3GuNzg1xWRp+elXAY+wL9svyV2aje7SzZ6PE+Dxb1qUXWbZl9J8skmLykqZk1hRFMtTtdJvKEUFeiJ9pp5HB22WTJ6TBxCWBLlxCrAyPemsy15dnjM7AAxfurGedQ7vagXFYSG9zebWFa2BawgIMYs4Zq6WnrborVBwb8B0E626fFlCoGm7X32kOf/QqryLKz+BK4qQEcP3JVtAd3ByMjOwzXIpR9Bs6tRBb/ucebP00phdPMhTOK6wLXpKpBNhICfP/qkIYGpqd2IWyZGjvPRmZCvy3KUXYB5sKnTMVzrbGLpv23gQKS8a2VYtCOtZIUO3s+aPZV7fVQJOv8Xe9bvzlOUJO21CFM2hh9RQMuvPrqTIGvejxU48GEDedeQrlREd790LuzvC0Kxu7BYVOsodeAWzAo6GLH6CXuDqLJWGERT1fdBjzrJoD0C1zhAstQma/iY5yzWAmQvYNMxS4eu1G7LlLo//pZD2AVoDC6Js5hS0ZNXxZimh9TkK9aT9LJgECfBr98ghWBZ0oOrW8ynYL+f8uwrfpRqefoOzY/rticGDN/+jZBh2HtG4I5xY389rS+QZeMVkUl3BS9BFktdmjsvg52QIlDlXoTHJcT+gXOj/enEX4c3kW1ifFJKj9UOE1fmDvV3ykOAP5viXGND8eLNgCXDHZUEByRFbx+bFP25sUiPuN2q7UB8/cJ3nnSQoMlT3TJ2PoDlZ7VjHujFlM9byL8xbyHZRLOS8561A1mOqAHi51FIHJ0i1L9jr/8F03Q8ZRCpVB96dYCNk3mKO06FC+VkeOgP3rkQmnDY75gaMVwwM6cEj77adiQ8teJgBlk6ZpFYkfDVbGBWt5+Kq585GdNpQFGEA24uGlO4iCpl+DDqp+dGiZXpkXmNpaV0VPiT3mhK9Ogd0U++82z6Jaql/gRB0vQ+z6uUmET5dh+JAf+KsmK7md/FtUnpClxKM9amoaBtnuLJm5vACec/3SPPzY2s9tNE1TrWwlsG6hqYGr6hCJGUT0xkE17O1dut3TU9UeAAIHqv2D0ZV+mRNTt+V0gLb0PKONeGCsWMcByvSUz8yaV92i7Q8hwIDBMmatB/OZQdCUlj4yD+6A8IRTgVU X-Microsoft-Exchange-Diagnostics: 1; VI1PR0802MB2495; 6:7F4in4Nd2WLqoDlGzxiZ/M4tR7ypxXXZNJPb0580+9JFJsOXI9AVnfqhkLHKRgAMvSc53+2i/UDJjta3ZRHvpkNZ1ifXusgej+VUfSk//clzamPyYG5g4hG1l+BwyYYEFv96O3Tg1h02t2a4CZ5ClhwnVF8IJ1q2QL67Mut1dfsgluC+MdKqAKYyjxzlv2tzbxKnT+tnTfnI2D3OS7AaY2ypqVeV7mVpQe1bXDeUa5BuPvzG/iisZhIaLpNgaR4FsH6MKyvoH5kWT0pyNmarmgkSSmA9E/kiGibUrBXuk9B647P2w40x1vAxvyzPX9nnG4khSt/9FZJliWKvRrIirQ==; 5:y7P3f+KK3ezK4k9sF3UNdJwdFY1TJQ+IXoyYrvGe4zw5WbxI0Yb25ABDRHp0iQQVi6Mb8RvAa0IcGRuGFAfd51CmHlsN6AtNxOQWOc86zxCWZOb6bknWgKjXUcGzWEDUTd9NmT8EIL31Qx//8D8P/w==; 24:HjP//1xm3kqdCLzz/+xFsjRAGWEGupcpGm32Zu61c37uuAtRjQ6HtdRkZ97/2nAr8bK0DKxYd4ebXrRjR02cmoy/nIJ3KP0MI/Y7FjJYndA=; 7:vFGQStFZZJ/co4bxedANbhd+7Fr/Vqq64Hso7/DmKkxV7sAalBilyZH8LD2H+eSOyFjdqDl9cm390ERQCZikoojCSBc9LZBoLxKZAjV03w5EsF23Okb9DXM/S7TK8eD61QCDcyg82NT5o5XRBtC+zUqbYKw3BJP4BUlK9+KMk6pdWuMnZp0Z8ccs885LCMfQq6a6zk/mQ6RoAWxpAcmBK9ShDu4erKrXo3DbF+VJVhM= SpamDiagnosticOutput: 1:99 SpamDiagnosticMetadata: NSPM X-OriginatorOrg: arm.com X-MS-Exchange-CrossTenant-OriginalArrivalTime: 19 Sep 2017 13:00:08.8361 (UTC) X-MS-Exchange-CrossTenant-FromEntityHeader: Hosted X-MS-Exchange-Transport-CrossTenantHeadersStamped: VI1PR0802MB2495 Similar to the new logf: double precision arithmetics and a small lookup table is used. The argument reduction step is the same as for the new logf. with naive ubenchmark (on aarch64): reciprocal-throughput: 2.7x faster latency: 2x faster old worst case error: 1.72 ulp new worst case error: 0.75 ulp aarch64 .text size: -252 bytes aarch64 .rodata size: +244 bytes 2017-09-19 Szabolcs Nagy * math/Makefile (type-float-routines): Add e_log2f_data. * sysdeps/ieee754/flt-32/e_log2f.c: New implementation. * sysdeps/ieee754/flt-32/e_log2f_data.c: New file. * sysdeps/ieee754/flt-32/math_config.h (__log2f_data): Define. (LOG2F_TABLE_BITS, LOG2F_POLY_ORDER): Define. * sysdeps/i386/fpu/e_log2f_data.c: New file. * sysdeps/ia64/fpu/e_log2f_data.c: New file. * sysdeps/m68k/m680x0/fpu/e_log2f_data.c: New file. diff --git a/math/Makefile b/math/Makefile index 919fec13ef..b4b3101592 100644 --- a/math/Makefile +++ b/math/Makefile @@ -115,7 +115,8 @@ type-double-routines := branred doasin dosincos halfulp mpa mpatan2 \ # float support type-float-suffix := f -type-float-routines := k_rem_pio2f math_errf e_exp2f_data e_logf_data +type-float-routines := k_rem_pio2f math_errf e_exp2f_data e_logf_data \ + e_log2f_data # _Float128 support type-float128-suffix := f128 diff --git a/sysdeps/i386/fpu/e_log2f_data.c b/sysdeps/i386/fpu/e_log2f_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/i386/fpu/e_log2f_data.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ia64/fpu/e_log2f_data.c b/sysdeps/ia64/fpu/e_log2f_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/ia64/fpu/e_log2f_data.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ieee754/flt-32/e_log2f.c b/sysdeps/ieee754/flt-32/e_log2f.c index 782d901094..6c42f27843 100644 --- a/sysdeps/ieee754/flt-32/e_log2f.c +++ b/sysdeps/ieee754/flt-32/e_log2f.c @@ -1,86 +1,88 @@ -/* e_logf.c -- float version of e_log.c. - * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com. - * adapted for log2 by Ulrich Drepper - */ +/* Single-precision log2 function. + Copyright (C) 2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. -/* - * ==================================================== - * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. - * - * Developed at SunPro, a Sun Microsystems, Inc. business. - * Permission to use, copy, modify, and distribute this - * software is freely granted, provided that this notice - * is preserved. - * ==================================================== - */ + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ #include -#include -#include +#include +#include "math_config.h" + +/* +LOG2F_TABLE_BITS = 4 +LOG2F_POLY_ORDER = 4 -static const float -ln2 = 0.69314718055994530942, -two25 = 3.355443200e+07, /* 0x4c000000 */ -Lg1 = 6.6666668653e-01, /* 3F2AAAAB */ -Lg2 = 4.0000000596e-01, /* 3ECCCCCD */ -Lg3 = 2.8571429849e-01, /* 3E924925 */ -Lg4 = 2.2222198546e-01, /* 3E638E29 */ -Lg5 = 1.8183572590e-01, /* 3E3A3325 */ -Lg6 = 1.5313838422e-01, /* 3E1CD04F */ -Lg7 = 1.4798198640e-01; /* 3E178897 */ +ULP error: 0.752 (nearest rounding.) +Relative error: 1.9 * 2^-26 (before rounding.) +*/ -static const float zero = 0.0; +#define N (1 << LOG2F_TABLE_BITS) +#define T __log2f_data.tab +#define A __log2f_data.poly +#define OFF 0x3f330000 float -__ieee754_log2f(float x) +__ieee754_log2f (float x) { - float hfsq,f,s,z,R,w,t1,t2,dk; - int32_t k,ix,i,j; + /* double_t for better performance on targets with FLT_EVAL_METHOD==2. */ + double_t z, r, r2, p, y, y0, invc, logc; + uint32_t ix, iz, top, tmp; + int k, i; + + ix = asuint (x); +#if WANT_ROUNDING + /* Fix sign of zero with downward rounding when x==1. */ + if (__glibc_unlikely (ix == 0x3f800000)) + return 0; +#endif + if (__glibc_unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) + { + /* x < 0x1p-126 or inf or nan. */ + if (ix * 2 == 0) + return __math_divzerof (1); + if (ix == 0x7f800000) /* log2(inf) == inf. */ + return x; + if ((ix & 0x80000000) || ix * 2 >= 0xff000000) + return __math_invalidf (x); + /* x is subnormal, normalize it. */ + ix = asuint (x * 0x1p23f); + ix -= 23 << 23; + } + + /* x = 2^k z; where z is in range [OFF,2*OFF] and exact. + The range is split into N subintervals. + The ith subinterval contains z and c is near its center. */ + tmp = ix - OFF; + i = (tmp >> (23 - LOG2F_TABLE_BITS)) % N; + top = tmp & 0xff800000; + iz = ix - top; + k = (int32_t) tmp >> 23; /* arithmetic shift */ + invc = T[i].invc; + logc = T[i].logc; + z = (double_t) asfloat (iz); - GET_FLOAT_WORD(ix,x); + /* log2(x) = log1p(z/c-1)/ln2 + log2(c) + k */ + r = z * invc - 1; + y0 = logc + (double_t) k; - k=0; - if (ix < 0x00800000) { /* x < 2**-126 */ - if (__builtin_expect((ix&0x7fffffff)==0, 0)) - return -two25/__fabsf (x); /* log(+-0)=-inf */ - if (__builtin_expect(ix<0, 0)) - return (x-x)/(x-x); /* log(-#) = NaN */ - k -= 25; x *= two25; /* subnormal number, scale up x */ - GET_FLOAT_WORD(ix,x); - } - if (__builtin_expect(ix >= 0x7f800000, 0)) return x+x; - k += (ix>>23)-127; - ix &= 0x007fffff; - i = (ix+(0x95f64<<3))&0x800000; - SET_FLOAT_WORD(x,ix|(i^0x3f800000)); /* normalize x or x/2 */ - k += (i>>23); - dk = (float)k; - f = x-(float)1.0; - if((0x007fffff&(15+ix))<16) { /* |f| < 2**-20 */ - if(f==zero) - { - if (FIX_INT_FP_CONVERT_ZERO && dk == 0.0f) - dk = 0.0f; - return dk; - } - R = f*f*((float)0.5-(float)0.33333333333333333*f); - return dk-(R-f)/ln2; - } - s = f/((float)2.0+f); - z = s*s; - i = ix-(0x6147a<<3); - w = z*z; - j = (0x6b851<<3)-ix; - t1= w*(Lg2+w*(Lg4+w*Lg6)); - t2= z*(Lg1+w*(Lg3+w*(Lg5+w*Lg7))); - i |= j; - R = t2+t1; - if(i>0) { - hfsq=(float)0.5*f*f; - return dk-((hfsq-(s*(hfsq+R)))-f)/ln2; - } else { - return dk-((s*(f-R))-f)/ln2; - } + /* Pipelined polynomial evaluation to approximate log1p(r)/ln2. */ + r2 = r * r; + y = A[1] * r + A[2]; + y = A[0] * r2 + y; + p = A[3] * r + y0; + y = y * r2 + p; + return (float) y; } strong_alias (__ieee754_log2f, __log2f_finite) diff --git a/sysdeps/ieee754/flt-32/e_log2f_data.c b/sysdeps/ieee754/flt-32/e_log2f_data.c new file mode 100644 index 0000000000..e39de3ba56 --- /dev/null +++ b/sysdeps/ieee754/flt-32/e_log2f_data.c @@ -0,0 +1,44 @@ +/* Data definition for log2f. + Copyright (C) 2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#include "math_config.h" + +const struct log2f_data __log2f_data = { + .tab = { + { 0x1.661ec79f8f3bep+0, -0x1.efec65b963019p-2 }, + { 0x1.571ed4aaf883dp+0, -0x1.b0b6832d4fca4p-2 }, + { 0x1.49539f0f010bp+0, -0x1.7418b0a1fb77bp-2 }, + { 0x1.3c995b0b80385p+0, -0x1.39de91a6dcf7bp-2 }, + { 0x1.30d190c8864a5p+0, -0x1.01d9bf3f2b631p-2 }, + { 0x1.25e227b0b8eap+0, -0x1.97c1d1b3b7afp-3 }, + { 0x1.1bb4a4a1a343fp+0, -0x1.2f9e393af3c9fp-3 }, + { 0x1.12358f08ae5bap+0, -0x1.960cbbf788d5cp-4 }, + { 0x1.0953f419900a7p+0, -0x1.a6f9db6475fcep-5 }, + { 0x1p+0, 0x0p+0 }, + { 0x1.e608cfd9a47acp-1, 0x1.338ca9f24f53dp-4 }, + { 0x1.ca4b31f026aap-1, 0x1.476a9543891bap-3 }, + { 0x1.b2036576afce6p-1, 0x1.e840b4ac4e4d2p-3 }, + { 0x1.9c2d163a1aa2dp-1, 0x1.40645f0c6651cp-2 }, + { 0x1.886e6037841edp-1, 0x1.88e9c2c1b9ff8p-2 }, + { 0x1.767dcf5534862p-1, 0x1.ce0a44eb17bccp-2 }, + }, + .poly = { + -0x1.712b6f70a7e4dp-2, 0x1.ecabf496832ep-2, -0x1.715479ffae3dep-1, + 0x1.715475f35c8b8p0, + } +}; diff --git a/sysdeps/ieee754/flt-32/math_config.h b/sysdeps/ieee754/flt-32/math_config.h index e1feafccfe..4e27b99bf4 100644 --- a/sysdeps/ieee754/flt-32/math_config.h +++ b/sysdeps/ieee754/flt-32/math_config.h @@ -123,4 +123,15 @@ extern const struct logf_data double poly[LOGF_POLY_ORDER]; } __logf_data; +#define LOG2F_TABLE_BITS 4 +#define LOG2F_POLY_ORDER 4 +extern const struct log2f_data +{ + struct + { + double invc, logc; + } tab[1 << LOG2F_TABLE_BITS]; + double poly[LOG2F_POLY_ORDER]; +} __log2f_data; + #endif diff --git a/sysdeps/m68k/m680x0/fpu/e_log2f_data.c b/sysdeps/m68k/m680x0/fpu/e_log2f_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/e_log2f_data.c @@ -0,0 +1 @@ +/* Not needed. */ From patchwork Tue Sep 19 13:01:45 2017 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit X-Patchwork-Submitter: Szabolcs Nagy X-Patchwork-Id: 815521 Return-Path: X-Original-To: incoming@patchwork.ozlabs.org Delivered-To: patchwork-incoming@bilbo.ozlabs.org Authentication-Results: ozlabs.org; spf=pass (mailfrom) smtp.mailfrom=sourceware.org (client-ip=209.132.180.131; helo=sourceware.org; envelope-from=libc-alpha-return-84727-incoming=patchwork.ozlabs.org@sourceware.org; receiver=) Authentication-Results: ozlabs.org; dkim=pass (1024-bit key; secure) header.d=sourceware.org header.i=@sourceware.org header.b="gBLHuZsG"; dkim-atps=neutral Received: from sourceware.org (server1.sourceware.org [209.132.180.131]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (No client certificate requested) by ozlabs.org (Postfix) with ESMTPS id 3xxNLv73Fqz9s7g for ; Tue, 19 Sep 2017 23:02:15 +1000 (AEST) DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; q=dns; s=default; b= Xc3Q6DmGRMA4FEpp4Lo6lhut0iobumvRs3KsOAdSCvDhnVMCtkP4/cJ7DnyyHAXA dMkf4Q69vm0Qkb608EbflzQ9/Vvekw0Qqn3MB5AgevRCF0cIawkrCrJQCBl5kEBk JtXSWXZA7FNH7o/05z15jMls6TvtqSmbzwfF1KdI4k0= DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; s=default; bh=/z74 qOuGBfRvZjgKrjHEmbTC6gQ=; b=gBLHuZsGl0tXc073W5kXPbqaTraY9l7RENfn gAVplIRfithSgEczaTwN8NWjKmabARzdQYVM57nS7rXoHbqJMHI5j/x2KXehCPI0 TNENUCPOX/60Z+83zNrYlk/tzWyr2ASbWu07t4JwudGCUN2W93bWQoKlaxsQ9wYq EaW+gYg= Received: (qmail 63594 invoked by alias); 19 Sep 2017 13:01:59 -0000 Mailing-List: contact libc-alpha-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Unsubscribe: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: libc-alpha-owner@sourceware.org Delivered-To: mailing list libc-alpha@sourceware.org Received: (qmail 63471 invoked by uid 89); 19 Sep 2017 13:01:58 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-23.9 required=5.0 tests=AWL, BAYES_00, GIT_PATCH_0, GIT_PATCH_1, GIT_PATCH_2, GIT_PATCH_3, KAM_ASCII_DIVIDERS, RCVD_IN_DNSWL_NONE, SPF_HELO_PASS, SPF_PASS autolearn=ham version=3.3.2 spammy=ski, zz X-HELO: EUR01-VE1-obe.outbound.protection.outlook.com Authentication-Results: spf=none (sender IP is ) smtp.mailfrom=Szabolcs.Nagy@arm.com; Message-ID: <59C11539.4020503@arm.com> Date: Tue, 19 Sep 2017 14:01:45 +0100 From: Szabolcs Nagy User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:31.0) Gecko/20100101 Thunderbird/31.8.0 MIME-Version: 1.0 To: GNU C Library CC: nd@arm.com Subject: [PATCH 6/9] New generic powf References: <59C1123F.9080003@arm.com> In-Reply-To: <59C1123F.9080003@arm.com> X-ClientProxiedBy: DB6PR0902CA0034.eurprd09.prod.outlook.com (2603:10a6:6:2::47) To AM5PR0802MB2481.eurprd08.prod.outlook.com (2603:10a6:203:97::23) X-MS-PublicTrafficType: Email X-MS-Office365-Filtering-Correlation-Id: 693679d2-e724-4fe6-3acb-08d4ff5e9698 X-MS-Office365-Filtering-HT: Tenant X-Microsoft-Antispam: UriScan:; BCL:0; PCL:0; RULEID:(300000500095)(300135000095)(300000501095)(300135300095)(300000502095)(300135100095)(22001)(2017030254152)(48565401081)(300000503095)(300135400095)(49563074)(201703131423075)(201703031133081)(201702281549075)(300000504095)(300135200095)(300000505095)(300135600095)(300000506095)(300135500095); SRVR:AM5PR0802MB2481; X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2481; 3:ViPZVlhvSZJpryFUT22hMEBoDpDyvmIScYVXA98xyOH1Ei2b8sKHVHi5m8un6jhm0XwAI27R1wr3Pve9te6lRmO9CJRJoq+7zrgQXmHinEfHIz3NpPitIBJIBbeUoEnRIxNhO7p5RWNDWImfL5kEMm9CaxDDAtWwEPMocXfTGhCBNbu0Zs4gOxQjx7z5FoBJitwCGHhb7JqaWMZquJe7C77LGN62/nTQkcy/H/LXw39LBIoxcGKE48hI38awoeZn; 25:gVxoDtzB6Gaxo6nX93do6dzOp3j8PgYKGkyxR/qJrRiwxQvtrpRnAgeLWEktO+NnOHWob8cgp53A1oMEJTynmq2/H+pkKFwqeMlEqZCjP8jCagSaXmYYdhEMpu8pVZ2LPBm7dIRN6YKcYHJiYXKxCJQvIyLNWmH4qon2lRN7Dt2fxvS3VznmMjJgLLbolOWUgNvI/2QFUyeeJVnak42L1mygeTjPeoMIn1kBfz7XbB9V3IpNuCTZBbxsSiMDI06BjqmzrN2n5Jc486or5gMvSpJ3vsm/Lq7GP9fkBjitGrRaUu68t+LoQA/CO3YfDtLo3Ghscl9edjbmA4fv3Hx14g==; 31:jIfuY7kzf1SZ4+U81rTiEzdzqR+IIHs0i9M9lJ3uY7jIunn/RhzGJymC9j6mG4fA65CadLNHjmvex8dbFiBsAFSPd8GJdS7BpFPNIoNOn4JFX4e8er0MDe2hsWKmrRLlUCG/L3VCR8yA4oSWBdOVRILD94C7/BMsmlnyZaMT4g36l1ZBZkx58mxDJsYiklDr3L5ksWmtmUJTRpVhIzaAAEaAT+IMneI5XHCXbv7HF0w= X-MS-TrafficTypeDiagnostic: AM5PR0802MB2481: NoDisclaimer: True X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2481; 20:y3WN4lEoU8TqzgJVIkbfbUbdUzMJr5eZhP3rzvTE4HLuqLzc9nCJQbBaPDRlv0jbM5S6jRROhH8u28roGcxZ5D+1lzNeAGXaN8pcxrhQEwW0CFkJpnXcfyST31mC6OptM17lnH1dp4dK0G0r01spnu9CV7YJaw3surow3cCZt2M=; 4:ujP0t71Yl/YhmPRhs3uykiE356QRN02uMQ3/jNirxQUIY2H+s7Ou/JI68L19QcOh9TvxhWcKxyqYZvI+qnFQp8Dzck+FRrJ3EXHKo6UwBu7HfjkplMdUON7ujaPRYqmnOs0Vep5iSO33J9VF3EnPlMVr1QeJzqrrygIYp1QdeVmrHedMQYxhqnn5yI7DucTaa2fDJjRsV1eogZyoAqWTs/EnTCBX30FpGktMIB4jxrF+ckKkpA4opy5kDJEtRCxZT3PMcMSGLchaqHZQVjUwh6nLBEtKJ2gS1y5m99Fex/4= X-Exchange-Antispam-Report-Test: UriScan:(180628864354917); X-Microsoft-Antispam-PRVS: X-Exchange-Antispam-Report-CFA-Test: BCL:0; PCL:0; RULEID:(100000700101)(100105000095)(100000701101)(100105300095)(100000702101)(100105100095)(102415395)(6040450)(2401047)(8121501046)(5005006)(93006095)(93001095)(10201501046)(100000703101)(100105400095)(3002001)(6055026)(6041248)(20161123560025)(20161123555025)(20161123558100)(20161123562025)(20161123564025)(201703131423075)(201702281528075)(201703061421075)(201703061406153)(6072148)(201708071742011)(100000704101)(100105200095)(100000705101)(100105500095); SRVR:AM5PR0802MB2481; BCL:0; PCL:0; RULEID:(100000800101)(100110000095)(100000801101)(100110300095)(100000802101)(100110100095)(100000803101)(100110400095)(100000804101)(100110200095)(100000805101)(100110500095); SRVR:AM5PR0802MB2481; X-Forefront-PRVS: 04359FAD81 X-Forefront-Antispam-Report: SFV:NSPM; SFS:(10009020)(4630300001)(7370300001)(6049001)(6009001)(346002)(39860400002)(376002)(377424004)(199003)(189002)(65956001)(65806001)(66066001)(53936002)(4326008)(25786009)(80316001)(105586002)(106356001)(5890100001)(68736007)(33656002)(305945005)(7736002)(81156014)(8676002)(81166006)(72206003)(478600001)(77096006)(6486002)(2476003)(6116002)(3846002)(5660300001)(6666003)(6916009)(270700001)(4610100001)(568964002)(2950100002)(21490400003)(36756003)(59896002)(189998001)(2906002)(58126008)(101416001)(16576012)(16526017)(50986999)(83506001)(16586007)(76176999)(84326002)(64126003)(86362001)(87266999)(54356999)(7350300001)(316002)(65816999)(97736004)(41533002); DIR:OUT; SFP:1101; SCL:1; SRVR:AM5PR0802MB2481; H:[10.2.206.69]; FPR:; SPF:None; PTR:InfoNoRecords; A:1; MX:1; LANG:en; Received-SPF: None (protection.outlook.com: arm.com does not designate permitted sender hosts) X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2481; 23:nM9Ndscqk5TF0OMGTvTfg1N/6NpzkrMvAP/D13vc698OiQcuUpAql35oxdNB3VulgIitNBc3z0GMEf6LVCzzb+DUbvV3BVIwjCg3XsGaExvOnhNEMsHgfw4PO/7QDl5w4z2OG0Sfgo7iTzCZ/lD8jNAa4Bps8dd8llYDhFjbA0300k5xwiXWwQKhsMpZl6Fwe3tsXT7ate5TKJNX2H+uVD8CGGT46laWFw4EZleOKoVBKHmXMidQxILFFLAZgDsbkRue/UpDj+LRDbVwdciq4oM4AGW7r2B9Cv2cQ+3zA6a6tJdv5wJZRdL7GAFbaGxCUlmxqj31CSlV4IE4BVyWyQ6iicJPoJa/ifhKFJvkL5A+Tx7b/wGiu/IyMQtHt0Ql3xI3x8znR+6XYuPRzwo2EhSVMIcPanaRxLFL3UPMURr9tSIwn6/wwmDDY5Vyp3ZPZEcUGzzj7yEz3EKudQpUf0bxTHJ0ks6GKgDyvZpbqlh8xTwkpDkD73gB1KCQi+qRiTpudOjQ8x0931BHMtkCSLoppoPUp69qiquSfwqLzsMwDijJVmyr97z2KczUybYp7Q0vokebLqdH6qwB2QIOyhFMyQjb3CmwQ+Yg+sFj2AgYyE4F80aIjkU6wntKFVJ02d/jV9MjVCfRBimpG1FbzlDIfLHUTXdQOJxMSvRFhepqhjpZuszIwrqfxtJTLw290TWPgLd/TDMkUCo4vZ4p/lZtDEEPOZjodoS6hBvRD3Sq98mE1HdsVv9X4E2ZT5vMuesKUv84it0GpZtidY1THL600PuC+HzgMWA4ll3Hnr4hgi1Si1pnS7Onx3L/vHn65xWr325/bU4W9iMmghqQDRRvYt6wP9WfmWNP7xWlvyEBOsSP9C/wtFixGhAGGjrU9tOw0eOew+/GxwjHyED5+o7ly/o3kSAU5upJhyjADGBM4q22VKkiSEVoWN8W6JYND+2JNHx6KA1u0F5qAVPPECjtCiHjWSOa3KxFVQiFYNqMWza+/h9OmKtoHsTkeKjjg2kidnz9OQRaz8PVT20G2cYK0OgWLq6octOFEU8gWeGxyHyAyzxM0Si/g35IP2APUr8zL2Huo560SznBiH5y0hwgUAyPJykNIoFM5Q+y3wxtf7jTzKDqvM2PnkT8ag2RT8XViYeT4yncvo1MXdGf+6g53bNGo65iFEKlmsA7AWcV8C/FvqueDmxrk5PVW9gGAdbJEKlqNdqmFOVjXRyi1oFTX69DBmbRtVwnHAEiRskdhCP+kLbYe4i2NLLTpSHd4f21tZUSfoF+rRkhRCU7K+K9XDIyFJexGMxcO5B0hv/qUmgJ2mYC58mMY1db9++CYxMZnYu3eYdvXOqANQjJldEupBvM7qXFbPbRKIr5J2Eokdf6RGMgHmunYl2SMunsCOjOa0qCnWj0bXegOQ3vZ4e78spXRQYCA4WePuBPnlCy6WKXLF+fZg8663AO7zy6nug+JSVCpq1qFx33P4XFDMMOnXGIjdJC1WZgPtwNTF4= X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2481; 6:90v5TYYTjKzm5YXKEsBM11Pl+7zm0q6RZIKg/+2d9a32yNGhCHtjXxaeU7L0d5uz9x+0Ga9lb0W09nE+Uf07yyqfBxogYDwJysUefRrP7ymsfbbOJN37AFfFbTneb8SEVBFZ2OF3CGTBdm/9pTrxO64yLUwWHN9GeupqoYqjpkG6CrnqWQt5LXdfrNoVWrzQG6qSSEGJgoxzdslanA+GDJ92a5widOBUNSPuyE9pk7JQrPjKRQ3ktj7oZCJ7+tRs/cX8ytECi+0tohqzQh2LBuealt7O2QVoNyHoSrSYYt5b1dX2QL/OZwJuv3D8vBDoqDKtmjwCyCcTuLTyQaOZKw==; 5:NqF2CxQrYuE/TyU22c+PHTHAa/amcnLG72VZVki2c/dOYE3bJOD/HHJCu524aaI5nSQyXeq6Hl0RN7JYBNOroz5IqKA7WsX9E4XZDqdvcVa51NkdNyVDocJ5kwC5jgWWmqNmAcakhG4fDPGBp2lQHw==; 24:MysX6pM4RNVTi1SMSIkaZnfPtP4DqPD6ILZI18NikWVhZQ0s0oohBM90iE5faRHJUK86IzMawFNR3uaWO+JE9x62BJp/uwbGPdqWbUVZG0w=; 7:R7dDmhwQvxsj/mhJGj5cIx7WeE+A15drbbuPyuGV/wJq2tsfAGfGjrigy1sUnrOTM6i8FkTPDMBWaugVgVstLWvXVzmCbIZ45oaKqJQt7AQpSq7ITeNQLOYYY27qSc2HtEKBDwwb4HHN/emPS/zmTaJ2xDB65Lc5FydT1mc8htwuMnT//r2sMhc3e3sxF0afxg3IGWxhJZj6mUMKJAO5D180rkra8LAH8J99arCrm4g= SpamDiagnosticOutput: 1:99 SpamDiagnosticMetadata: NSPM X-OriginatorOrg: arm.com X-MS-Exchange-CrossTenant-OriginalArrivalTime: 19 Sep 2017 13:01:48.1837 (UTC) X-MS-Exchange-CrossTenant-FromEntityHeader: Hosted X-MS-Exchange-Transport-CrossTenantHeadersStamped: AM5PR0802MB2481 with powf benchmark (on aarch64): reciprocal-throughput: 4.2x faster latency: 2.6x faster old worst-case error: 1.11 ulp new worst-case error: 0.82 ulp aarch64 .text size: -780 bytes aarch64 .rodata size: +144 bytes powf(x,y) is implemented as exp2(y*log2(x)) with the same algorithms that are used in exp2f and log2f, except that the log2f polynomial is larger for extra precision and its output (and exp2f input) may be scaled by a power of 2 (POWF_SCALE) to simplify the argument reduction step of exp2 (possible when efficient round and convert toint operation is available). The special case handling tries to minimize the checks in the hot path. When the input of exp2_inline is checked, int arithmetics is used as that was faster on the tested aarch64 cores. 2017-09-19 Szabolcs Nagy * math/Makefile (type-float-routines): Add e_powf_log2_data. * sysdeps/ieee754/flt-32/e_powf.c: New implementation. * sysdeps/ieee754/flt-32/e_powf_log2_data.c: New file. * sysdeps/ieee754/flt-32/math_config.h (__powf_data): Define. (issignalingf_inline): Likewise. (POWF_LOG2_TABLE_BITS): Likewise. (POWF_LOG2_POLY_ORDER): Likewise. (POWF_SCALE_BITS): Likewise. (POWF_SCALE): Likewise. * sysdeps/i386/fpu/e_powf_log2_data.c: New file. * sysdeps/ia64/fpu/e_powf_log2_data.c: New file. * sysdeps/m68k/m680x0/fpu/e_powf_log2_data.c: New file. diff --git a/math/Makefile b/math/Makefile index b4b3101592..6c8aa3e413 100644 --- a/math/Makefile +++ b/math/Makefile @@ -116,7 +116,7 @@ type-double-routines := branred doasin dosincos halfulp mpa mpatan2 \ # float support type-float-suffix := f type-float-routines := k_rem_pio2f math_errf e_exp2f_data e_logf_data \ - e_log2f_data + e_log2f_data e_powf_log2_data # _Float128 support type-float128-suffix := f128 diff --git a/sysdeps/i386/fpu/e_powf_log2_data.c b/sysdeps/i386/fpu/e_powf_log2_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/i386/fpu/e_powf_log2_data.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ia64/fpu/e_powf_log2_data.c b/sysdeps/ia64/fpu/e_powf_log2_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/ia64/fpu/e_powf_log2_data.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ieee754/flt-32/e_powf.c b/sysdeps/ieee754/flt-32/e_powf.c index ce8e11f1ea..644a18d05e 100644 --- a/sysdeps/ieee754/flt-32/e_powf.c +++ b/sysdeps/ieee754/flt-32/e_powf.c @@ -1,7 +1,5 @@ -/* e_powf.c -- float version of e_pow.c. - * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com. - */ -/* Copyright (C) 2017 Free Software Foundation, Inc. +/* Single-precision pow function. + Copyright (C) 2017 Free Software Foundation, Inc. This file is part of the GNU C Library. The GNU C Library is free software; you can redistribute it and/or @@ -18,210 +16,202 @@ License along with the GNU C Library; if not, see . */ -/* - * ==================================================== - * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. - * - * Developed at SunPro, a Sun Microsystems, Inc. business. - * Permission to use, copy, modify, and distribute this - * software is freely granted, provided that this notice - * is preserved. - * ==================================================== - */ - #include -#include - -static const float huge = 1.0e+30, tiny = 1.0e-30; - -static const float -bp[] = {1.0, 1.5,}, -zero = 0.0, -one = 1.0, -two = 2.0, -two24 = 16777216.0, /* 0x4b800000 */ - /* poly coefs for (3/2)*(log(x)-2s-2/3*s**3 */ -L1 = 6.0000002384e-01, /* 0x3f19999a */ -L2 = 4.2857143283e-01, /* 0x3edb6db7 */ -L3 = 3.3333334327e-01, /* 0x3eaaaaab */ -L4 = 2.7272811532e-01, /* 0x3e8ba305 */ -L5 = 2.3066075146e-01, /* 0x3e6c3255 */ -L6 = 2.0697501302e-01, /* 0x3e53f142 */ -P1 = 1.6666667163e-01, /* 0x3e2aaaab */ -P2 = -2.7777778450e-03, /* 0xbb360b61 */ -P3 = 6.6137559770e-05, /* 0x388ab355 */ -P4 = -1.6533901999e-06, /* 0xb5ddea0e */ -P5 = 4.1381369442e-08, /* 0x3331bb4c */ -ovt = 4.2995665694e-08; /* -(128-log2(ovfl+.5ulp)) */ - -static const double - dp[] = { 0.0, 0x1.2b803473f7ad1p-1, }, /* log2(1.5) */ - lg2 = M_LN2, - cp = 2.0/3.0/M_LN2, - invln2 = 1.0/M_LN2; +#include +#include "math_config.h" -float -__ieee754_powf(float x, float y) +/* +POWF_LOG2_POLY_ORDER = 5 +EXP2F_TABLE_BITS = 5 + +ULP error: 0.82 (~ 0.5 + relerr*2^24) +relerr: 1.27 * 2^-26 (Relative error ~= 128*Ln2*relerr_log2 + relerr_exp2) +relerr_log2: 1.83 * 2^-33 (Relative error of logx.) +relerr_exp2: 1.69 * 2^-34 (Relative error of exp2(ylogx).) +*/ + +#define N (1 << POWF_LOG2_TABLE_BITS) +#define T __powf_log2_data.tab +#define A __powf_log2_data.poly +#define OFF 0x3f330000 + +/* Subnormal input is normalized so ix has negative biased exponent. + Output is multiplied by N (POWF_SCALE) if TOINT_INTRINICS is set. */ +static inline double_t +log2_inline (uint32_t ix) { - float z, ax, s; - double d1, d2; - int32_t i,j,k,yisint,n; - int32_t hx,hy,ix,iy; - - GET_FLOAT_WORD(hy,y); - iy = hy&0x7fffffff; - - /* y==zero: x**0 = 1 */ - if(iy==0 && !issignaling (x)) return one; - - /* x==+-1 */ - if(x == 1.0 && !issignaling (y)) return one; - if(x == -1.0 && isinf(y)) return one; - - GET_FLOAT_WORD(hx,x); - ix = hx&0x7fffffff; - - /* +-NaN return x+y */ - if(__builtin_expect(ix > 0x7f800000 || - iy > 0x7f800000, 0)) - return x+y; - - /* special value of y */ - if (__builtin_expect(iy==0x7f800000, 0)) { /* y is +-inf */ - if (ix==0x3f800000) - return y - y; /* inf**+-1 is NaN */ - else if (ix > 0x3f800000)/* (|x|>1)**+-inf = inf,0 */ - return (hy>=0)? y: zero; - else /* (|x|<1)**-,+inf = inf,0 */ - return (hy<0)?-y: zero; - } - if(iy==0x3f800000) { /* y is +-1 */ - if(hy<0) return one/x; else return x; - } - if(hy==0x40000000) return x*x; /* y is 2 */ - if(hy==0x3f000000) { /* y is 0.5 */ - if(__builtin_expect(hx>=0, 1)) /* x >= +0 */ - return __ieee754_sqrtf(x); - } + /* double_t for better performance on targets with FLT_EVAL_METHOD==2. */ + double_t z, r, r2, r4, p, q, y, y0, invc, logc; + uint32_t iz, top, tmp; + int k, i; + + /* x = 2^k z; where z is in range [OFF,2*OFF] and exact. + The range is split into N subintervals. + The ith subinterval contains z and c is near its center. */ + tmp = ix - OFF; + i = (tmp >> (23 - POWF_LOG2_TABLE_BITS)) % N; + top = tmp & 0xff800000; + iz = ix - top; + k = (int32_t) top >> (23 - POWF_SCALE_BITS); /* arithmetic shift */ + invc = T[i].invc; + logc = T[i].logc; + z = (double_t) asfloat (iz); + + /* log2(x) = log1p(z/c-1)/ln2 + log2(c) + k */ + r = z * invc - 1; + y0 = logc + (double_t) k; + + /* Pipelined polynomial evaluation to approximate log1p(r)/ln2. */ + r2 = r * r; + y = A[0] * r + A[1]; + p = A[2] * r + A[3]; + r4 = r2 * r2; + q = A[4] * r + y0; + q = p * r2 + q; + y = y * r4 + q; + return y; +} - /* determine if y is an odd int when x < 0 - * yisint = 0 ... y is not an integer - * yisint = 1 ... y is an odd int - * yisint = 2 ... y is an even int - */ - yisint = 0; - if(hx<0) { - if(iy>=0x4b800000) yisint = 2; /* even integer y */ - else if(iy>=0x3f800000) { - k = (iy>>23)-0x7f; /* exponent */ - j = iy>>(23-k); - if((j<<(23-k))==iy) yisint = 2-(j&1); - } - } +#undef N +#undef T +#define N (1 << EXP2F_TABLE_BITS) +#define T __exp2f_data.tab +#define SIGN_BIAS (1 << (EXP2F_TABLE_BITS + 11)) + +/* The output of log2 and thus the input of exp2 is either scaled by N + (in case of fast toint intrinsics) or not. The unscaled xd must be + in [-1021,1023], sign_bias sets the sign of the result. */ +static inline double_t +exp2_inline (double_t xd, unsigned long sign_bias) +{ + uint64_t ki, ski, t; + /* double_t for better performance on targets with FLT_EVAL_METHOD==2. */ + double_t kd, z, r, r2, y, s; + +#if TOINT_INTRINSICS +# define C __exp2f_data.poly_scaled + /* N*x = k + r with r in [-1/2, 1/2] */ + kd = roundtoint (xd); /* k */ + ki = converttoint (xd); +#else +# define C __exp2f_data.poly +# define SHIFT __exp2f_data.shift_scaled + /* x = k/N + r with r in [-1/(2N), 1/(2N)] */ + kd = (double) (xd + SHIFT); /* Rounding to double precision is required. */ + ki = asuint64 (kd); + kd -= SHIFT; /* k/N */ +#endif + r = xd - kd; + + /* exp2(x) = 2^(k/N) * 2^r ~= s * (C0*r^3 + C1*r^2 + C2*r + 1) */ + t = T[ki % N]; + ski = ki + sign_bias; + t += ski << (52 - EXP2F_TABLE_BITS); + s = asdouble (t); + z = C[0] * r + C[1]; + r2 = r * r; + y = C[2] * r + 1; + y = z * r2 + y; + y = y * s; + return y; +} - ax = fabsf(x); - /* special value of x */ - if(__builtin_expect(ix==0x7f800000||ix==0||ix==0x3f800000, 0)){ - z = ax; /*x is +-0,+-inf,+-1*/ - if(hy<0) z = one/z; /* z = (1/|x|) */ - if(hx<0) { - if(((ix-0x3f800000)|yisint)==0) { - z = (z-z)/(z-z); /* (-1)**non-int is NaN */ - } else if(yisint==1) - z = -z; /* (x<0)**odd = -(|x|**odd) */ - } - return z; - } +/* Returns 0 if not int, 1 if odd int, 2 if even int. */ +static inline int +checkint (uint32_t iy) +{ + int e = iy >> 23 & 0xff; + if (e < 0x7f) + return 0; + if (e > 0x7f + 23) + return 2; + if (iy & ((1 << (0x7f + 23 - e)) - 1)) + return 0; + if (iy & (1 << (0x7f + 23 - e))) + return 1; + return 2; +} - /* (x<0)**(non-int) is NaN */ - if(__builtin_expect(((((uint32_t)hx>>31)-1)|yisint)==0, 0)) - return (x-x)/(x-x); - - /* |y| is huge */ - if(__builtin_expect(iy>0x4d000000, 0)) { /* if |y| > 2**27 */ - /* over/underflow if x is not close to one */ - if(ix<0x3f7ffff8) return (hy<0)? huge*huge:tiny*tiny; - if(ix>0x3f800007) return (hy>0)? huge*huge:tiny*tiny; - /* now |1-x| is tiny <= 2**-20, suffice to compute - log(x) by x-x^2/2+x^3/3-x^4/4 */ - d2 = ax-1; /* d2 has 20 trailing zeros. */ - d2 = d2 * invln2 - - (d2 * d2) * (0.5 - d2 * (0.333333333333 - d2 * 0.25)) * invln2; - } else { - /* Avoid internal underflow for tiny y. The exact value - of y does not matter if |y| <= 2**-32. */ - if (iy < 0x2f800000) - SET_FLOAT_WORD (y, (hy & 0x80000000) | 0x2f800000); - n = 0; - /* take care subnormal number */ - if(ix<0x00800000) - {ax *= two24; n -= 24; GET_FLOAT_WORD(ix,ax); } - n += ((ix)>>23)-0x7f; - j = ix&0x007fffff; - /* determine interval */ - ix = j|0x3f800000; /* normalize ix */ - if(j<=0x1cc471) k=0; /* |x|= 2u * 0x7f800000 - 1; +} - s = one; /* s (sign of result -ve**odd) = -1 else = 1 */ - if(((((uint32_t)hx>>31)-1)|(yisint-1))==0) - s = -one; /* (-ve)**(odd int) */ - - /* compute y * d2 */ - d1 = y * d2; - z = d1; - GET_FLOAT_WORD(j,z); - if (__builtin_expect(j>0x43000000, 0)) /* if z > 128 */ - return s*huge*huge; /* overflow */ - else if (__builtin_expect(j==0x43000000, 0)) { /* if z == 128 */ - if(ovt>(z-d1)) return s*huge*huge; /* overflow */ +float +__ieee754_powf (float x, float y) +{ + unsigned long sign_bias = 0; + uint32_t ix, iy; + + ix = asuint (x); + iy = asuint (y); + if (__glibc_unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000 + || zeroinfnan (iy))) + { + /* Either (x < 0x1p-126 or inf or nan) or (y is 0 or inf or nan). */ + if (__glibc_unlikely (zeroinfnan (iy))) + { + if (2 * iy == 0) + return issignalingf_inline (x) ? x + y : 1.0f; + if (ix == 0x3f800000) + return issignalingf_inline (y) ? x + y : 1.0f; + if (2 * ix > 2u * 0x7f800000 || 2 * iy > 2u * 0x7f800000) + return x + y; + if (2 * ix == 2 * 0x3f800000) + return 1.0f; + if ((2 * ix < 2 * 0x3f800000) == !(iy & 0x80000000)) + return 0.0f; /* |x|<1 && y==inf or |x|>1 && y==-inf. */ + return y * y; + } + if (__glibc_unlikely (zeroinfnan (ix))) + { + float_t x2 = x * x; + if (ix & 0x80000000 && checkint (iy) == 1) + { + x2 = -x2; + sign_bias = 1; + } +#if WANT_ERRNO + if (2 * ix == 0 && iy & 0x80000000) + return __math_divzerof (sign_bias); +#endif + return iy & 0x80000000 ? 1 / x2 : x2; } - else if (__builtin_expect((j&0x7fffffff)>0x43160000, 0))/* z <= -150 */ - return s*tiny*tiny; /* underflow */ - else if (__builtin_expect((uint32_t) j==0xc3160000, 0)){/* z == -150*/ - if(0.0<=(z-d1)) return s*tiny*tiny; /* underflow */ + /* x and y are non-zero finite. */ + if (ix & 0x80000000) + { + /* Finite x < 0. */ + int yint = checkint (iy); + if (yint == 0) + return __math_invalidf (x); + if (yint == 1) + sign_bias = SIGN_BIAS; + ix &= 0x7fffffff; } - /* - * compute 2**d1 - */ - i = j&0x7fffffff; - k = (i>>23)-0x7f; - n = 0; - if(i>0x3f000000) { /* if |z| > 0.5, set n = [z+0.5] */ - n = j+(0x00800000>>(k+1)); - k = ((n&0x7fffffff)>>23)-0x7f; /* new k for n */ - SET_FLOAT_WORD(z,n&~(0x007fffff>>k)); - n = ((n&0x007fffff)|0x00800000)>>(23-k); - if(j<0) n = -n; - d1 -= z; + if (ix < 0x00800000) + { + /* Normalize subnormal x so exponent becomes negative. */ + ix = asuint (x * 0x1p23f); + ix &= 0x7fffffff; + ix -= 23 << 23; } - d1 = d1 * lg2; - d2 = d1*d1; - d2 = d1 - d2*(P1+d2*(P2+d2*(P3+d2*(P4+d2*P5)))); - d2 = (d1*d2)/(d2-two); - z = one - (d2-d1); - GET_FLOAT_WORD(j,z); - j += (n<<23); - if((j>>23)<=0) /* subnormal output */ - { - z = __scalbnf (z, n); - float force_underflow = z * z; - math_force_eval (force_underflow); - } - else SET_FLOAT_WORD(z,j); - return s*z; + } + double_t logx = log2_inline (ix); + double_t ylogx = y * logx; /* Note: cannot overflow, y is single prec. */ + if (__glibc_unlikely ((asuint64 (ylogx) >> 47 & 0xffff) + >= asuint64 (126.0 * POWF_SCALE) >> 47)) + { + /* |y*log(x)| >= 126. */ + if (ylogx > 0x1.fffffffd1d571p+6 * POWF_SCALE) + return __math_oflowf (sign_bias); + if (ylogx <= -150.0 * POWF_SCALE) + return __math_uflowf (sign_bias); +#if WANT_ERRNO_UFLOW + if (ylogx < -149.0 * POWF_SCALE) + return __math_may_uflowf (sign_bias); +#endif + } + return (float) exp2_inline (ylogx, sign_bias); } strong_alias (__ieee754_powf, __powf_finite) diff --git a/sysdeps/ieee754/flt-32/e_powf_log2_data.c b/sysdeps/ieee754/flt-32/e_powf_log2_data.c new file mode 100644 index 0000000000..7cff06f59b --- /dev/null +++ b/sysdeps/ieee754/flt-32/e_powf_log2_data.c @@ -0,0 +1,45 @@ +/* Data definition for powf. + Copyright (C) 2017 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#include "math_config.h" + +const struct powf_log2_data __powf_log2_data = { + .tab = { + { 0x1.661ec79f8f3bep+0, -0x1.efec65b963019p-2 * POWF_SCALE }, + { 0x1.571ed4aaf883dp+0, -0x1.b0b6832d4fca4p-2 * POWF_SCALE }, + { 0x1.49539f0f010bp+0, -0x1.7418b0a1fb77bp-2 * POWF_SCALE }, + { 0x1.3c995b0b80385p+0, -0x1.39de91a6dcf7bp-2 * POWF_SCALE }, + { 0x1.30d190c8864a5p+0, -0x1.01d9bf3f2b631p-2 * POWF_SCALE }, + { 0x1.25e227b0b8eap+0, -0x1.97c1d1b3b7afp-3 * POWF_SCALE }, + { 0x1.1bb4a4a1a343fp+0, -0x1.2f9e393af3c9fp-3 * POWF_SCALE }, + { 0x1.12358f08ae5bap+0, -0x1.960cbbf788d5cp-4 * POWF_SCALE }, + { 0x1.0953f419900a7p+0, -0x1.a6f9db6475fcep-5 * POWF_SCALE }, + { 0x1p+0, 0x0p+0 * POWF_SCALE }, + { 0x1.e608cfd9a47acp-1, 0x1.338ca9f24f53dp-4 * POWF_SCALE }, + { 0x1.ca4b31f026aap-1, 0x1.476a9543891bap-3 * POWF_SCALE }, + { 0x1.b2036576afce6p-1, 0x1.e840b4ac4e4d2p-3 * POWF_SCALE }, + { 0x1.9c2d163a1aa2dp-1, 0x1.40645f0c6651cp-2 * POWF_SCALE }, + { 0x1.886e6037841edp-1, 0x1.88e9c2c1b9ff8p-2 * POWF_SCALE }, + { 0x1.767dcf5534862p-1, 0x1.ce0a44eb17bccp-2 * POWF_SCALE }, + }, + .poly = { + 0x1.27616c9496e0bp-2 * POWF_SCALE, -0x1.71969a075c67ap-2 * POWF_SCALE, + 0x1.ec70a6ca7baddp-2 * POWF_SCALE, -0x1.7154748bef6c8p-1 * POWF_SCALE, + 0x1.71547652ab82bp0 * POWF_SCALE, + } +}; diff --git a/sysdeps/ieee754/flt-32/math_config.h b/sysdeps/ieee754/flt-32/math_config.h index 4e27b99bf4..2dcc5a57d7 100644 --- a/sysdeps/ieee754/flt-32/math_config.h +++ b/sysdeps/ieee754/flt-32/math_config.h @@ -21,6 +21,7 @@ #include #include +#include #include #ifndef WANT_ROUNDING @@ -90,6 +91,15 @@ asdouble (uint64_t i) return u.f; } +static inline int +issignalingf_inline (float x) +{ + uint32_t ix = asuint (x); + if (HIGH_ORDER_BIT_IS_SET_FOR_SNAN) + return (ix & 0x7fc00000) == 0x7fc00000; + return 2 * (ix ^ 0x00400000) > 2u * 0x7fc00000; +} + #define NOINLINE __attribute__ ((noinline)) attribute_hidden float __math_oflowf (unsigned long); @@ -134,4 +144,21 @@ extern const struct log2f_data double poly[LOG2F_POLY_ORDER]; } __log2f_data; +#define POWF_LOG2_TABLE_BITS 4 +#define POWF_LOG2_POLY_ORDER 5 +#if TOINT_INTRINSICS +#define POWF_SCALE_BITS EXP2F_TABLE_BITS +#else +#define POWF_SCALE_BITS 0 +#endif +#define POWF_SCALE ((double) (1 << POWF_SCALE_BITS)) +extern const struct powf_log2_data +{ + struct + { + double invc, logc; + } tab[1 << POWF_LOG2_TABLE_BITS]; + double poly[POWF_LOG2_POLY_ORDER]; +} __powf_log2_data; + #endif diff --git a/sysdeps/m68k/m680x0/fpu/e_powf_log2_data.c b/sysdeps/m68k/m680x0/fpu/e_powf_log2_data.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/e_powf_log2_data.c @@ -0,0 +1 @@ +/* Not needed. */ From patchwork Tue Sep 19 13:03:33 2017 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit X-Patchwork-Submitter: Szabolcs Nagy X-Patchwork-Id: 815523 Return-Path: X-Original-To: incoming@patchwork.ozlabs.org Delivered-To: patchwork-incoming@bilbo.ozlabs.org Authentication-Results: ozlabs.org; spf=pass (mailfrom) smtp.mailfrom=sourceware.org (client-ip=209.132.180.131; helo=sourceware.org; envelope-from=libc-alpha-return-84728-incoming=patchwork.ozlabs.org@sourceware.org; receiver=) Authentication-Results: ozlabs.org; dkim=pass (1024-bit key; secure) header.d=sourceware.org header.i=@sourceware.org header.b="DaGnu3zU"; dkim-atps=neutral Received: from sourceware.org (server1.sourceware.org [209.132.180.131]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (No client certificate requested) by ozlabs.org (Postfix) with ESMTPS id 3xxNNp6GKfz9s7g for ; Tue, 19 Sep 2017 23:03:54 +1000 (AEST) DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; q=dns; s=default; b= aFnd8AiTmoNap32JKZHnvnHr9s0qpjf2Va7WbH6/YT7WEqXBrEa/Sd1Q/BYzG/Ih 8AzRxnOwvgqYbS394LFTrL28c9PnPPYwtl/lpf5VYPA28wbdVRH0vimqVK3e10QY MgSJF0AT/JLBUlCkYISlrbfLW80mGtSDaP2odrZCREs= DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; s=default; bh=FTdN fDhiArqjg9GmuP3UbaHzAJQ=; b=DaGnu3zU21uUD+6N06Sw8jps36ftjujx82Fs GJwtG/R6djXQkuUp1E1GRDPCHZqKiZcdnNbLbNF+FRR4QmOp+Salp8xwrfLyNY0/ vzzJiOKY1WOHda1l0pp1WNd1cDLeE84VhItB2LqDe0vH5VZiOc5f7ZYyvRXRALbO DWONcyA= Received: (qmail 66851 invoked by alias); 19 Sep 2017 13:03:46 -0000 Mailing-List: contact libc-alpha-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Unsubscribe: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: libc-alpha-owner@sourceware.org Delivered-To: mailing list libc-alpha@sourceware.org Received: (qmail 66837 invoked by uid 89); 19 Sep 2017 13:03:45 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-24.6 required=5.0 tests=AWL, BAYES_00, GIT_PATCH_0, GIT_PATCH_1, GIT_PATCH_2, GIT_PATCH_3, RCVD_IN_DNSWL_NONE, SPF_HELO_PASS, SPF_PASS autolearn=ham version=3.3.2 spammy= X-HELO: EUR03-DB5-obe.outbound.protection.outlook.com Authentication-Results: spf=none (sender IP is ) smtp.mailfrom=Szabolcs.Nagy@arm.com; Message-ID: <59C115A5.5050800@arm.com> Date: Tue, 19 Sep 2017 14:03:33 +0100 From: Szabolcs Nagy User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:31.0) Gecko/20100101 Thunderbird/31.8.0 MIME-Version: 1.0 To: GNU C Library CC: nd@arm.com Subject: [PATCH 7/9] New symbol version for logf, log2f and powf without SVID compat References: <59C1123F.9080003@arm.com> In-Reply-To: <59C1123F.9080003@arm.com> X-ClientProxiedBy: AM5PR06CA0034.eurprd06.prod.outlook.com (2603:10a6:206:2::47) To AM5PR0802MB2484.eurprd08.prod.outlook.com (2603:10a6:203:9b::23) X-MS-PublicTrafficType: Email X-MS-Office365-Filtering-Correlation-Id: 860c6d7e-2a96-461d-c20c-08d4ff5ed672 X-MS-Office365-Filtering-HT: Tenant X-Microsoft-Antispam: UriScan:; BCL:0; PCL:0; RULEID:(300000500095)(300135000095)(300000501095)(300135300095)(300000502095)(300135100095)(22001)(2017030254152)(48565401081)(300000503095)(300135400095)(49563074)(201703131423075)(201703031133081)(201702281549075)(300000504095)(300135200095)(300000505095)(300135600095)(300000506095)(300135500095); SRVR:AM5PR0802MB2484; X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2484; 3:oFazd08aEu7Tp4UcdcJkp/osxIocKyGaRApQmW7LCTgjvK3tMKchO5JIhu0riACr3J8Xwcb51lJvkbJRisTYBgur/+ZR6zPdJ2zCNo5LLb9lSaBX4irY1yX9g6mJSgxP1lZlDGS9FboMFZc8MlaLl7mwkuQEXvXwcvaMC7oBLzG1IP2AbhmBe91NfGkOhQpO0ogmo7skSdfXTdS8Unqe4hbTSpiwDZJyfktbYJsCzWL5M1p4Lt+KZN19u5ds9XBl; 25:/VsFEsR/E0oEaxG0Xu7o0fFxkqLQKjXZeMdu+S/j71w2yb2kKKZGGiDw67/e2JFwv5wCYBVJ0NuOKbmpANcuCtmWVH6z/pboBl4f8T70nV+2JA3rbHWfu2bS3AM+Nh+xuhMwSuOMmoK7hlQPCvgHd3Uw9KG7swOS0Fhw7p1Ijgmcb656IYwJUp+t3LoajY0YR/mgwJa6mccPJAoiBCiYs5vxKPSHuHDx6TIL1sESyCR6/RqWCr/r7bD3xnyIvAly+zlTMPHNtteg50sYg6x6fkunbj+HndFa5Q4bW1gDIIer3YLPgjvBHikue5OKW6tviEeAOKMHdMZhGScfBYyZzQ==; 31:2YlkndIw5K6OJWC/L++WssK7cdoMUMTQjK6zqW7UDLLFlYY4mtmdkIm6aEVul6eZpilJj+He2MA+nv/Nrke7oCYtTvNq5qvRPccLvRBA3mLcrq/sMGml3DTQjrs1V4PyfkxRHlzD+uKyBWbkXlea4kEaNzLjz34yiAFzOXtebzWxAJxythCfz0x2frKoLYohkwD5k79ExUZMmiv9D1dAWsx09zOruZWX1yvDrD4n2Iw= X-MS-TrafficTypeDiagnostic: AM5PR0802MB2484: NoDisclaimer: True X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2484; 20:olaUCtlq1BEwo2sKvq4fuGh0kNYvh3zq/A9No3+ktKDoUKSBkFNUlN66abLI1GUNLwfHHYZ0eRXXkRfbrGCxx/Srvtv+q7k69DaofGYaHXGX3HkYo32vbjodQHPqj/oceJlJtic90BcUa9SdkxgckdXwn4kdfd954GmOgwI+X7g=; 4:devN8gGz8P2Jz+fYwtAY8s1NpGZfVYqhOFplZQGCpikhweJ105eJTTZZQ2hOXnFFDVVOrAWgwfx/OFNBSRvjjVksR4KAryw5PTn50EmSi1v27gIUnHHxJDJwGMLAsXgg/LOWKPE7TRRNAu5FPsCj0ylVYUfW094ddSfFYDZHLV34WtCBN34Bqoak654EQ5oS29dKA9k0zp0k8W55T4gYhDo9V1+OPyI18rhGNqPwhuPiHtMH4tbs2Ug0szmKqazHskOScj1bZciFrDVNe+GO0Bnl1D1hm8imH8CUzZPk/Bk= X-Exchange-Antispam-Report-Test: UriScan:(180628864354917); X-Microsoft-Antispam-PRVS: X-Exchange-Antispam-Report-CFA-Test: BCL:0; PCL:0; RULEID:(100000700101)(100105000095)(100000701101)(100105300095)(100000702101)(100105100095)(102415395)(6040450)(2401047)(8121501046)(5005006)(93006095)(93001095)(100000703101)(100105400095)(3002001)(10201501046)(6055026)(6041248)(20161123560025)(20161123558100)(20161123562025)(20161123564025)(20161123555025)(201703131423075)(201702281528075)(201703061421075)(201703061406153)(6072148)(201708071742011)(100000704101)(100105200095)(100000705101)(100105500095); SRVR:AM5PR0802MB2484; BCL:0; PCL:0; RULEID:(100000800101)(100110000095)(100000801101)(100110300095)(100000802101)(100110100095)(100000803101)(100110400095)(100000804101)(100110200095)(100000805101)(100110500095); SRVR:AM5PR0802MB2484; X-Forefront-PRVS: 04359FAD81 X-Forefront-Antispam-Report: SFV:NSPM; SFS:(10009020)(7370300001)(4630300001)(6009001)(6049001)(39860400002)(376002)(346002)(199003)(377424004)(189002)(86362001)(4610100001)(2950100002)(77096006)(6486002)(6916009)(76176999)(65816999)(84326002)(87266999)(478600001)(54356999)(36756003)(83506001)(189998001)(50986999)(101416001)(5660300001)(72206003)(65806001)(25786009)(66066001)(65956001)(5890100001)(68736007)(105586002)(59896002)(106356001)(81166006)(81156014)(2906002)(4326008)(64126003)(3846002)(6116002)(305945005)(2476003)(7350300001)(8676002)(270700001)(80316001)(16526017)(316002)(16576012)(97736004)(16586007)(7736002)(33656002)(58126008)(568964002)(21490400003)(53936002); DIR:OUT; SFP:1101; SCL:1; SRVR:AM5PR0802MB2484; H:[10.2.206.69]; FPR:; SPF:None; PTR:InfoNoRecords; A:1; MX:1; LANG:en; Received-SPF: None (protection.outlook.com: arm.com does not designate permitted sender hosts) X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2484; 23:VEbVAlOBJnvux6gYGcX89vkwb4FMk0Qypf3B5p73kTJlzLIQfFx1hMO/ecvlxUD6L/YmQXKaKtfayooVc8q9Zx9MeKZugIMrAeF6euIPt72cMNk6yz+tGtbl04aPpEkMH+FCJHGIk8SQ83ng86QJZaBagcPDuGxYAabxBiyGH/bfYXqS8G6Qpptm6x2qic2xdklidy8pNABQtkAWxzvyNQJwi4wxAz7TILqeXKGpXB5uWiYJ0illkBUB3KRYMH/F6Uqxxw6Um/9mHoF4gAUFaP1tHw6xXOyGxoWNEKgbngb5jGBqJpD+4cqS50FNt8qz5VdGhpXLI4HjIoYqquHwjCCKxQ0206RJJNleQzrYLUr4pfqXq3/oqclauv344nM9A93Zrz5z/2/DmodR8r9d9a8WSECymKv7LbiMP0308M9Rq8AUxiwlUsHpRB4QRDpN0qkVDeuZfGsvvxCKpKiAeVd+H8SAwCTtmdIZCmdAwWUuKuwoLRwVVjGNTJ5HIW9lpqdCpy6aR8GAXz5Y00DUf/S6RDo2mSeoiYV6sNeWCUOLXgy4lej2L3RpXhTdWDXw6COkFfMHWhIGqN6ptKLeowud6Ab644QWbH+HPt0e8adriEQkmMmPRPbQhZ3Ey8qdcEnGLjRKDWyxL4uTtMeyIVOdL+bAhlGSUJUDwwDJd4W5geBYtfEe6C0/nOvb1YN0lYFCVBqwbzi23pkAOjZ2lyZWQPEshJTemgnL5X/qqN19IZqa82rth5F6TKnjrKIvfjc05cMc+saQaTrDCkbBsap2rByQIZlWiXfCMpdNncVg/vVPGvY77WvaUBhHB/s28wm/wieRvVO+Pl85HZ0FZqS8k48ONoR98SMhv8vPpuTLbw1G33N9s3azSRv/EMQc9lFtOIfSc4HjKlJV7etbMz7bI6HqBvwAph2/J6VQ8uBSWlHOsNPkwJhn8EA8wJEDWbPximrQ+akzeM3ahBvMN9aH/S4p5uGkwURVEiCOGsfechiSOsYYLR+gnHSHaxDk4okgzIDIOp2BoQyw1oTcB8Tap49Tt0z1P84wn+MUG7kisFqSnVjvlGrvMGgU6vQUY4ftD0yvkYhJ/FOzQuZ9MgfFIDMCWWxFYEM5xLXwlTXM/3xssVsJkNjHwryygUxqjydBNnQ7oB3G2+v6dYpMocNzFJXDHM6zSh6hHsi3hZIYf4hPB5i3iRQqeJVLtRpYRcJAfMklpOhoGJg4dcXrzCGW6ju//dKRWsWTOF/V8rhkwv3mtTrhSGyHzZlcF7v/ChQTVNMl4HfaO8MxbBSQHSJuR/04I0SdIp+NEcQ5GOTTy6pHMxI97+RUrTSxEJSdOaEHSFd/zeGlAVWSVm8Kfrna2ktOctLqoR0oyKj0X5hgfF5/xu4DTZjamiDP//NP425cQ8lfgMGx2qFzVn+wJjtaXDTwqvd8dHigdj+sGSvqxrkYk/c5jdSFsRPAU962 X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2484; 6:2xT9BqifsyxRdH0fjgSAdJjU248hJnlIJKnMq2lmIrQZp3bACvtMOJ0M1MNlml2CW5pADWQA5HuUsa2gIAel/F8oe48x2g3CT0G28zN7PnkFvu3JIcq2aey/fZHyVkn011ak+dR0sSGNRfaFBaT1D+GXwLboGOFUjE2SC3MmOEwQ8uYXF9AooO1xsNaMp2RAQYXW1Cx8WzQCIvWKmq4Ae8IU6G1+dye2/kURDaqF+MshJiJA5UdIFhwWYHSu6ixwdrVfkZfmJTdC/0cC0El7NXIAJTgSWzwBjpeYrvTyXu5KKK77JANi7ZMa8qYEjGe4d5m0MWNU9k13k39ptvxs0w==; 5:/96/DJiAvjjMIOr4TqfMShATSJWZSFfPHQJF/dl/4pRj+xUwre5+tfID/tinnhoBGKso0KujRoEL3F1GYC8pU/+s34fCxnP2eW9q8i34GT8Sg5BX+ofofETcjGmgg+FGyTPN7bfA1oid7E7inkJM2g==; 24:dPmkGLycyeNUVceSArZ3Ge23ZMNepI/6WvhYM//Eovg57NqFlKKOXTtciOE62xFEwpD5ytMHU2sHamrbW+au+ClQUNc/+NiB7oZcUj/Ld6o=; 7:tCp5/rOUU6FdZN1UBGt9QFLP0PDBjEdEMhF8d5gpOsAK1nBIz1vkBmrgJgAKw/3T/4rERxbQYnxBdKo+JuHvpVaN5sYG8qk3Q5J2xO+bF3oYEUSHNMXO+6kP8iU86Q9f0W9xXM0njPvAb27hRxfykJp0qQdJrEBWa4kb9g1p4b+H6Bm85KgE2k6+m/mCSr01YxfRqndjSGWuo2MnRm1zbXJQg28er+YoJZ2c2b+ltM8= SpamDiagnosticOutput: 1:99 SpamDiagnosticMetadata: NSPM X-OriginatorOrg: arm.com X-MS-Exchange-CrossTenant-OriginalArrivalTime: 19 Sep 2017 13:03:35.4686 (UTC) X-MS-Exchange-CrossTenant-FromEntityHeader: Hosted X-MS-Exchange-Transport-CrossTenantHeadersStamped: AM5PR0802MB2484 This patch changes the logf, log2f and powf error handling semantics to only set errno accoring to POSIX rules. New symbol version is introduced at GLIBC_2.27. The old wrappers are kept for compat symbols. ia64 needed assembly change to have the new and compat versioned symbol map to the same function. All linux libm abilists are updated. 2017-09-19 Szabolcs Nagy * math/Versions (logf): New libm symbol at GLIBC_2.27. (log2f): Likewise. (powf): Likewise. * math/w_log2f.c: New file. * math/w_logf.c: New file. * math/w_powf.c: New file. * math/w_log2f_compat.c (__log2f_compat): For compat symbol only. * math/w_logf_compat.c (__logf_compat): Likewise. * math/w_powf_compat.c (__powf_compat): Likewise. * sysdeps/ia64/fpu/e_log2f.S: Add versioned symbols. * sysdeps/ia64/fpu/e_logf.S: Likewise. * sysdeps/ia64/fpu/e_powf.S: Likewise. * sysdeps/unix/sysv/linux/aarch64/libm.abilist: Update. * sysdeps/unix/sysv/linux/alpha/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/arm/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/hppa/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/i386/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/ia64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/m68k/coldfire/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/m68k/m680x0/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/microblaze/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/mips/mips32/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/mips/mips64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/nios2/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/powerpc/powerpc32/fpu/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/powerpc/powerpc32/nofpu/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/powerpc/powerpc64/libm-le.abilist: Likewise. * sysdeps/unix/sysv/linux/powerpc/powerpc64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/s390/s390-32/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/s390/s390-64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/sh/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/sparc/sparc32/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/sparc/sparc64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/tile/tilegx/tilegx32/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/tile/tilegx/tilegx64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/tile/tilepro/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/x86_64/64/libm.abilist: Likewise. * sysdeps/unix/sysv/linux/x86_64/x32/libm.abilist: Likewise. diff --git a/math/Versions b/math/Versions index 380f6a2a1a..2fbdb2f266 100644 --- a/math/Versions +++ b/math/Versions @@ -230,6 +230,6 @@ libm { fromfpx; fromfpxf; fromfpxl; ufromfpx; ufromfpxf; ufromfpxl; } GLIBC_2.27 { - expf; exp2f; + expf; exp2f; logf; log2f; powf; } } diff --git a/math/w_log2f.c b/math/w_log2f.c new file mode 100644 index 0000000000..cda0c3a644 --- /dev/null +++ b/math/w_log2f.c @@ -0,0 +1,7 @@ +#include +#undef __USE_WRAPPER_TEMPLATE +#define __USE_WRAPPER_TEMPLATE 1 +#undef declare_mgen_alias +#define declare_mgen_alias(a, b) +#include +versioned_symbol (libm, __log2f, log2f, GLIBC_2_27); diff --git a/math/w_log2f_compat.c b/math/w_log2f_compat.c index 295c1620f7..3caa310c51 100644 --- a/math/w_log2f_compat.c +++ b/math/w_log2f_compat.c @@ -23,10 +23,10 @@ #include -#if LIBM_SVID_COMPAT +#if LIBM_SVID_COMPAT && SHLIB_COMPAT (libm, GLIBC_2_1, GLIBC_2_27) /* wrapper log2f(x) */ float -__log2f (float x) +__log2f_compat (float x) { if (__builtin_expect (islessequal (x, 0.0f), 0) && _LIB_VERSION != _IEEE_) { @@ -44,5 +44,5 @@ __log2f (float x) return __ieee754_log2f (x); } -libm_alias_float (__log2, log2) +compat_symbol (libm, __log2f_compat, log2f, GLIBC_2_1); #endif diff --git a/math/w_logf.c b/math/w_logf.c new file mode 100644 index 0000000000..d960e016d7 --- /dev/null +++ b/math/w_logf.c @@ -0,0 +1,7 @@ +#include +#undef __USE_WRAPPER_TEMPLATE +#define __USE_WRAPPER_TEMPLATE 1 +#undef declare_mgen_alias +#define declare_mgen_alias(a, b) +#include +versioned_symbol (libm, __logf, logf, GLIBC_2_27); diff --git a/math/w_logf_compat.c b/math/w_logf_compat.c index 7cdacdf921..936b3a6e67 100644 --- a/math/w_logf_compat.c +++ b/math/w_logf_compat.c @@ -23,10 +23,10 @@ #include -#if LIBM_SVID_COMPAT +#if LIBM_SVID_COMPAT && SHLIB_COMPAT (libm, GLIBC_2_0, GLIBC_2_27) /* wrapper logf(x) */ float -__logf (float x) +__logf_compat (float x) { if (__builtin_expect (islessequal (x, 0.0f), 0) && _LIB_VERSION != _IEEE_) { @@ -44,5 +44,5 @@ __logf (float x) return __ieee754_logf (x); } -libm_alias_float (__log, log) +compat_symbol (libm, __logf_compat, logf, GLIBC_2_0); #endif diff --git a/math/w_powf.c b/math/w_powf.c new file mode 100644 index 0000000000..a18348329e --- /dev/null +++ b/math/w_powf.c @@ -0,0 +1,7 @@ +#include +#undef __USE_WRAPPER_TEMPLATE +#define __USE_WRAPPER_TEMPLATE 1 +#undef declare_mgen_alias +#define declare_mgen_alias(a, b) +#include +versioned_symbol (libm, __powf, powf, GLIBC_2_27); diff --git a/math/w_powf_compat.c b/math/w_powf_compat.c index 39e818af7e..7745639efe 100644 --- a/math/w_powf_compat.c +++ b/math/w_powf_compat.c @@ -22,10 +22,10 @@ #include -#if LIBM_SVID_COMPAT +#if LIBM_SVID_COMPAT && SHLIB_COMPAT (libm, GLIBC_2_0, GLIBC_2_27) /* wrapper powf */ float -__powf (float x, float y) +__powf_compat (float x, float y) { float z = __ieee754_powf (x, y); if (__glibc_unlikely (!isfinite (z))) @@ -60,5 +60,5 @@ __powf (float x, float y) return z; } -libm_alias_float (__pow, pow) +compat_symbol (libm, __powf_compat, powf, GLIBC_2_0); #endif diff --git a/sysdeps/ia64/fpu/e_log2f.S b/sysdeps/ia64/fpu/e_log2f.S index 2c3f18f360..9b754d1043 100644 --- a/sysdeps/ia64/fpu/e_log2f.S +++ b/sysdeps/ia64/fpu/e_log2f.S @@ -252,7 +252,7 @@ LOCAL_OBJECT_END(T_table) .section .text -GLOBAL_LIBM_ENTRY(log2f) +GLOBAL_LIBM_ENTRY(__log2f) { .mfi alloc r32=ar.pfs,1,4,4,0 @@ -491,7 +491,13 @@ SPECIAL_log2f: br.ret.sptk b0;; } -GLOBAL_LIBM_END(log2f) +GLOBAL_LIBM_END(__log2f) +#ifdef SHARED +.symver __log2f,log2f@@GLIBC_2.27 +.weak __log2f_compat +.set __log2f_compat,__log2f +.symver __log2f_compat,log2f@GLIBC_2.2 +#endif LOCAL_LIBM_ENTRY(__libm_error_region) diff --git a/sysdeps/ia64/fpu/e_logf.S b/sysdeps/ia64/fpu/e_logf.S index 2dda2186d0..d5f5437793 100644 --- a/sysdeps/ia64/fpu/e_logf.S +++ b/sysdeps/ia64/fpu/e_logf.S @@ -1088,6 +1088,12 @@ logf_libm_err: nop.i 0 };; GLOBAL_IEEE754_END(logf) +#ifdef SHARED +.symver logf,logf@@GLIBC_2.27 +.weak __logf_compat +.set __logf_compat,__logf +.symver __logf_compat,logf@GLIBC_2.2 +#endif // Stack operations when calling error support. diff --git a/sysdeps/ia64/fpu/e_powf.S b/sysdeps/ia64/fpu/e_powf.S index d61bc79e5e..388391624f 100644 --- a/sysdeps/ia64/fpu/e_powf.S +++ b/sysdeps/ia64/fpu/e_powf.S @@ -868,7 +868,7 @@ data8 0xEAC0C6E7DD24392F , 0x00003FFF LOCAL_OBJECT_END(pow_tbl2) .section .text -GLOBAL_LIBM_ENTRY(powf) +GLOBAL_LIBM_ENTRY(__powf) // Get exponent of x. Will be used to calculate K. { .mfi @@ -2002,7 +2002,13 @@ POW_OVER_UNDER_ERROR: } ;; -GLOBAL_LIBM_END(powf) +GLOBAL_LIBM_END(__powf) +#ifdef SHARED +.symver __powf,powf@@GLIBC_2.27 +.weak __powf_compat +.set __powf_compat,__powf +.symver __powf_compat,powf@GLIBC_2.2 +#endif LOCAL_LIBM_ENTRY(__libm_error_region) diff --git a/sysdeps/unix/sysv/linux/aarch64/libm.abilist b/sysdeps/unix/sysv/linux/aarch64/libm.abilist index 10102eeaff..3f0190ae03 100644 --- a/sysdeps/unix/sysv/linux/aarch64/libm.abilist +++ b/sysdeps/unix/sysv/linux/aarch64/libm.abilist @@ -463,3 +463,6 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F diff --git a/sysdeps/unix/sysv/linux/alpha/libm.abilist b/sysdeps/unix/sysv/linux/alpha/libm.abilist index e09a115aa9..78edc5e3d9 100644 --- a/sysdeps/unix/sysv/linux/alpha/libm.abilist +++ b/sysdeps/unix/sysv/linux/alpha/libm.abilist @@ -473,6 +473,9 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.3.4 GLIBC_2.3.4 A GLIBC_2.3.4 __c1_cabsf F GLIBC_2.3.4 __c1_cacosf F diff --git a/sysdeps/unix/sysv/linux/arm/libm.abilist b/sysdeps/unix/sysv/linux/arm/libm.abilist index 8095876449..b3fd4a27b2 100644 --- a/sysdeps/unix/sysv/linux/arm/libm.abilist +++ b/sysdeps/unix/sysv/linux/arm/libm.abilist @@ -120,6 +120,9 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A GLIBC_2.4 _LIB_VERSION D 0x4 GLIBC_2.4 __clog10 F diff --git a/sysdeps/unix/sysv/linux/hppa/libm.abilist b/sysdeps/unix/sysv/linux/hppa/libm.abilist index 19d40ef50d..ffa61bf1b1 100644 --- a/sysdeps/unix/sysv/linux/hppa/libm.abilist +++ b/sysdeps/unix/sysv/linux/hppa/libm.abilist @@ -432,5 +432,8 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A GLIBC_2.4 exp2l F diff --git a/sysdeps/unix/sysv/linux/i386/libm.abilist b/sysdeps/unix/sysv/linux/i386/libm.abilist index 791fba28e5..1a7e6bf449 100644 --- a/sysdeps/unix/sysv/linux/i386/libm.abilist +++ b/sysdeps/unix/sysv/linux/i386/libm.abilist @@ -614,4 +614,7 @@ GLIBC_2.26 ynf128 F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A diff --git a/sysdeps/unix/sysv/linux/ia64/libm.abilist b/sysdeps/unix/sysv/linux/ia64/libm.abilist index 65a0fbe56a..7e15735eae 100644 --- a/sysdeps/unix/sysv/linux/ia64/libm.abilist +++ b/sysdeps/unix/sysv/linux/ia64/libm.abilist @@ -543,4 +543,7 @@ GLIBC_2.26 ynf128 F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A diff --git a/sysdeps/unix/sysv/linux/m68k/coldfire/libm.abilist b/sysdeps/unix/sysv/linux/m68k/coldfire/libm.abilist index 8095876449..b3fd4a27b2 100644 --- a/sysdeps/unix/sysv/linux/m68k/coldfire/libm.abilist +++ b/sysdeps/unix/sysv/linux/m68k/coldfire/libm.abilist @@ -120,6 +120,9 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A GLIBC_2.4 _LIB_VERSION D 0x4 GLIBC_2.4 __clog10 F diff --git a/sysdeps/unix/sysv/linux/m68k/m680x0/libm.abilist b/sysdeps/unix/sysv/linux/m68k/m680x0/libm.abilist index 5e692dda7b..aae61169f9 100644 --- a/sysdeps/unix/sysv/linux/m68k/m680x0/libm.abilist +++ b/sysdeps/unix/sysv/linux/m68k/m680x0/libm.abilist @@ -474,4 +474,7 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A diff --git a/sysdeps/unix/sysv/linux/microblaze/libm.abilist b/sysdeps/unix/sysv/linux/microblaze/libm.abilist index 65f1d5b451..0d3b4b1e90 100644 --- a/sysdeps/unix/sysv/linux/microblaze/libm.abilist +++ b/sysdeps/unix/sysv/linux/microblaze/libm.abilist @@ -431,3 +431,6 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F diff --git a/sysdeps/unix/sysv/linux/mips/mips32/libm.abilist b/sysdeps/unix/sysv/linux/mips/mips32/libm.abilist index c32ea5b96a..d32d58d4e4 100644 --- a/sysdeps/unix/sysv/linux/mips/mips32/libm.abilist +++ b/sysdeps/unix/sysv/linux/mips/mips32/libm.abilist @@ -433,6 +433,9 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A GLIBC_2.4 exp2l F _gp_disp _gp_disp A diff --git a/sysdeps/unix/sysv/linux/mips/mips64/libm.abilist b/sysdeps/unix/sysv/linux/mips/mips64/libm.abilist index 18b2aa2404..f33ba0576b 100644 --- a/sysdeps/unix/sysv/linux/mips/mips64/libm.abilist +++ b/sysdeps/unix/sysv/linux/mips/mips64/libm.abilist @@ -465,4 +465,7 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A diff --git a/sysdeps/unix/sysv/linux/nios2/libm.abilist b/sysdeps/unix/sysv/linux/nios2/libm.abilist index e492a68e9d..0fe34e98fa 100644 --- a/sysdeps/unix/sysv/linux/nios2/libm.abilist +++ b/sysdeps/unix/sysv/linux/nios2/libm.abilist @@ -431,3 +431,6 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F diff --git a/sysdeps/unix/sysv/linux/powerpc/powerpc32/fpu/libm.abilist b/sysdeps/unix/sysv/linux/powerpc/powerpc32/fpu/libm.abilist index ad8f0372f4..ed013deefd 100644 --- a/sysdeps/unix/sysv/linux/powerpc/powerpc32/fpu/libm.abilist +++ b/sysdeps/unix/sysv/linux/powerpc/powerpc32/fpu/libm.abilist @@ -476,6 +476,9 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A GLIBC_2.4 __clog10l F GLIBC_2.4 __finitel F diff --git a/sysdeps/unix/sysv/linux/powerpc/powerpc32/nofpu/libm.abilist b/sysdeps/unix/sysv/linux/powerpc/powerpc32/nofpu/libm.abilist index 9c26b5b809..6f2873dc80 100644 --- a/sysdeps/unix/sysv/linux/powerpc/powerpc32/nofpu/libm.abilist +++ b/sysdeps/unix/sysv/linux/powerpc/powerpc32/nofpu/libm.abilist @@ -475,6 +475,9 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A GLIBC_2.4 __clog10l F GLIBC_2.4 __finitel F diff --git a/sysdeps/unix/sysv/linux/powerpc/powerpc64/libm-le.abilist b/sysdeps/unix/sysv/linux/powerpc/powerpc64/libm-le.abilist index 8e36699f28..723be46c20 100644 --- a/sysdeps/unix/sysv/linux/powerpc/powerpc64/libm-le.abilist +++ b/sysdeps/unix/sysv/linux/powerpc/powerpc64/libm-le.abilist @@ -608,3 +608,6 @@ GLIBC_2.26 ynf128 F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F diff --git a/sysdeps/unix/sysv/linux/powerpc/powerpc64/libm.abilist b/sysdeps/unix/sysv/linux/powerpc/powerpc64/libm.abilist index 9ca0c3ccfc..f3aeac2e1e 100644 --- a/sysdeps/unix/sysv/linux/powerpc/powerpc64/libm.abilist +++ b/sysdeps/unix/sysv/linux/powerpc/powerpc64/libm.abilist @@ -151,6 +151,9 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.3 GLIBC_2.3 A GLIBC_2.3 _LIB_VERSION D 0x4 GLIBC_2.3 __clog10 F diff --git a/sysdeps/unix/sysv/linux/s390/s390-32/libm.abilist b/sysdeps/unix/sysv/linux/s390/s390-32/libm.abilist index 8a79f0137f..2b758e80fd 100644 --- a/sysdeps/unix/sysv/linux/s390/s390-32/libm.abilist +++ b/sysdeps/unix/sysv/linux/s390/s390-32/libm.abilist @@ -463,6 +463,9 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A GLIBC_2.4 __clog10l F GLIBC_2.4 __finitel F diff --git a/sysdeps/unix/sysv/linux/s390/s390-64/libm.abilist b/sysdeps/unix/sysv/linux/s390/s390-64/libm.abilist index df81853618..62c9bb57a8 100644 --- a/sysdeps/unix/sysv/linux/s390/s390-64/libm.abilist +++ b/sysdeps/unix/sysv/linux/s390/s390-64/libm.abilist @@ -461,6 +461,9 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A GLIBC_2.4 __clog10l F GLIBC_2.4 __finitel F diff --git a/sysdeps/unix/sysv/linux/sh/libm.abilist b/sysdeps/unix/sysv/linux/sh/libm.abilist index 6b6a42dc9c..a57fbc0eac 100644 --- a/sysdeps/unix/sysv/linux/sh/libm.abilist +++ b/sysdeps/unix/sysv/linux/sh/libm.abilist @@ -432,5 +432,8 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A GLIBC_2.4 exp2l F diff --git a/sysdeps/unix/sysv/linux/sparc/sparc32/libm.abilist b/sysdeps/unix/sysv/linux/sparc/sparc32/libm.abilist index 24d67d22e1..f8f10e5952 100644 --- a/sysdeps/unix/sysv/linux/sparc/sparc32/libm.abilist +++ b/sysdeps/unix/sysv/linux/sparc/sparc32/libm.abilist @@ -467,6 +467,9 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A GLIBC_2.4 __clog10l F GLIBC_2.4 __finitel F diff --git a/sysdeps/unix/sysv/linux/sparc/sparc64/libm.abilist b/sysdeps/unix/sysv/linux/sparc/sparc64/libm.abilist index 2fdccc0de3..b5412c9b48 100644 --- a/sysdeps/unix/sysv/linux/sparc/sparc64/libm.abilist +++ b/sysdeps/unix/sysv/linux/sparc/sparc64/libm.abilist @@ -464,4 +464,7 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A diff --git a/sysdeps/unix/sysv/linux/tile/tilegx/tilegx32/libm.abilist b/sysdeps/unix/sysv/linux/tile/tilegx/tilegx32/libm.abilist index 98bc348f91..b711e87026 100644 --- a/sysdeps/unix/sysv/linux/tile/tilegx/tilegx32/libm.abilist +++ b/sysdeps/unix/sysv/linux/tile/tilegx/tilegx32/libm.abilist @@ -432,3 +432,6 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F diff --git a/sysdeps/unix/sysv/linux/tile/tilegx/tilegx64/libm.abilist b/sysdeps/unix/sysv/linux/tile/tilegx/tilegx64/libm.abilist index 98bc348f91..b711e87026 100644 --- a/sysdeps/unix/sysv/linux/tile/tilegx/tilegx64/libm.abilist +++ b/sysdeps/unix/sysv/linux/tile/tilegx/tilegx64/libm.abilist @@ -432,3 +432,6 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F diff --git a/sysdeps/unix/sysv/linux/tile/tilepro/libm.abilist b/sysdeps/unix/sysv/linux/tile/tilepro/libm.abilist index 98bc348f91..b711e87026 100644 --- a/sysdeps/unix/sysv/linux/tile/tilepro/libm.abilist +++ b/sysdeps/unix/sysv/linux/tile/tilepro/libm.abilist @@ -432,3 +432,6 @@ GLIBC_2.25 ufromfpxl F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F diff --git a/sysdeps/unix/sysv/linux/x86_64/64/libm.abilist b/sysdeps/unix/sysv/linux/x86_64/64/libm.abilist index e6fd3fe3df..201c2ab1b1 100644 --- a/sysdeps/unix/sysv/linux/x86_64/64/libm.abilist +++ b/sysdeps/unix/sysv/linux/x86_64/64/libm.abilist @@ -603,4 +603,7 @@ GLIBC_2.26 ynf128 F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F GLIBC_2.4 GLIBC_2.4 A diff --git a/sysdeps/unix/sysv/linux/x86_64/x32/libm.abilist b/sysdeps/unix/sysv/linux/x86_64/x32/libm.abilist index afa7b98697..10e389a96a 100644 --- a/sysdeps/unix/sysv/linux/x86_64/x32/libm.abilist +++ b/sysdeps/unix/sysv/linux/x86_64/x32/libm.abilist @@ -602,3 +602,6 @@ GLIBC_2.26 ynf128 F GLIBC_2.27 GLIBC_2.27 A GLIBC_2.27 exp2f F GLIBC_2.27 expf F +GLIBC_2.27 log2f F +GLIBC_2.27 logf F +GLIBC_2.27 powf F From patchwork Tue Sep 19 13:04:57 2017 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit X-Patchwork-Submitter: Szabolcs Nagy X-Patchwork-Id: 815526 Return-Path: X-Original-To: incoming@patchwork.ozlabs.org Delivered-To: patchwork-incoming@bilbo.ozlabs.org Authentication-Results: ozlabs.org; spf=pass (mailfrom) smtp.mailfrom=sourceware.org (client-ip=209.132.180.131; helo=sourceware.org; envelope-from=libc-alpha-return-84729-incoming=patchwork.ozlabs.org@sourceware.org; receiver=) Authentication-Results: ozlabs.org; dkim=pass (1024-bit key; secure) header.d=sourceware.org header.i=@sourceware.org header.b="Q2+FnJ1R"; dkim-atps=neutral Received: from sourceware.org (server1.sourceware.org [209.132.180.131]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (No client certificate requested) by ozlabs.org (Postfix) with ESMTPS id 3xxNQK0gPqz9s7g for ; Tue, 19 Sep 2017 23:05:12 +1000 (AEST) DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; q=dns; s=default; b= o1ap5al743lEcocvp0aaQIK6MbKzIm7h3dcM086g5QEeTDocP28GkLrYU0hJdbnW FzgnoR5tiaSmOzHa+JAI5rmlySpjZ+y/3rVYm7VgXGrK96qOW1984UsCYDDEy53i hG80Kn2eqkb0Ibsq2BpmLKeCLeOC8Or79qdsfGZPjdw= DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; s=default; bh=DflM KDt8hUialwwkoeOuJ9eaK5E=; b=Q2+FnJ1RePi5t9GOvnBM2uLZhmTogIYTvkHN tdCAwLVY5XlzXnEgGkZ5wLFDCzcaorJh+9btQXsCk7mDIyB3MZEyRm/lZKMhehh9 H0ii8WyfYn2CWA5KvZ6Y6VhEzhepPIsQ3eiBEImZX5g1xa9qWK/o4/puhnfeIbH2 rMEkRmA= Received: (qmail 70169 invoked by alias); 19 Sep 2017 13:05:05 -0000 Mailing-List: contact libc-alpha-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Unsubscribe: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: libc-alpha-owner@sourceware.org Delivered-To: mailing list libc-alpha@sourceware.org Received: (qmail 70139 invoked by uid 89); 19 Sep 2017 13:05:04 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-24.7 required=5.0 tests=AWL, BAYES_00, GIT_PATCH_0, GIT_PATCH_1, GIT_PATCH_2, GIT_PATCH_3, RCVD_IN_DNSWL_NONE, SPF_HELO_PASS, SPF_PASS autolearn=ham version=3.3.2 spammy= X-HELO: EUR01-VE1-obe.outbound.protection.outlook.com Authentication-Results: spf=none (sender IP is ) smtp.mailfrom=Szabolcs.Nagy@arm.com; Message-ID: <59C115F9.8010800@arm.com> Date: Tue, 19 Sep 2017 14:04:57 +0100 From: Szabolcs Nagy User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:31.0) Gecko/20100101 Thunderbird/31.8.0 MIME-Version: 1.0 To: GNU C Library CC: nd@arm.com Subject: [PATCH 8/9] Do not wrap logf, log2f and powf References: <59C1123F.9080003@arm.com> In-Reply-To: <59C1123F.9080003@arm.com> X-ClientProxiedBy: AM5PR06CA0001.eurprd06.prod.outlook.com (2603:10a6:206:2::14) To AM5PR0802MB2481.eurprd08.prod.outlook.com (2603:10a6:203:97::23) X-MS-PublicTrafficType: Email X-MS-Office365-Filtering-Correlation-Id: dd822979-9611-4e93-5ca8-08d4ff5f084f X-MS-Office365-Filtering-HT: Tenant X-Microsoft-Antispam: UriScan:; BCL:0; PCL:0; RULEID:(300000500095)(300135000095)(300000501095)(300135300095)(300000502095)(300135100095)(22001)(2017030254152)(48565401081)(300000503095)(300135400095)(49563074)(201703131423075)(201703031133081)(201702281549075)(300000504095)(300135200095)(300000505095)(300135600095)(300000506095)(300135500095); SRVR:AM5PR0802MB2481; X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2481; 3:ieeUt3fDsdL8xkf/EYCYLYffEbOJKVybXOThEw51ithGm68jYXtRyQWkJ2iNL2YENwyxBSS3jwE4U07xHzfuE+IqCc7RWmuFhLdoPeBvi0XnkKGgO++2tbYQyoRj4uPrQDDsNNjM2DWJnR3e/1qfl5xoeO6rWlltIsRsAquB252QTaCx3899YPA27yWC+batUGXCaYm/CAkB+xFWT8zBX/weqPXdRTQjFbkWEzzGfv57tt4OGn5cgUYt6nBV9lhT; 25:DUeD+0iM0oxZzkPfCO+UvHzWKYfQdfCul5lDm36TO3TczKBkMuB20rAduMzflkKE2eVSUtu7o8TcUV229m0M1jXontkN/cxhpjRZolbqbgCw/pSszHwXpqehhCpMhTBpKOF/X83rKpfO4Snrav4pPi0BrpvRHjvFalbVZwk+j+A4uvtn7+H7YyOFIlF1yz5HYAEHm+SI5TzFCOc2KDBDcyQjD48UiHoZdrSzWvx3eIRXp9r9PrB1c/SyCroIKJ97cgBLPNRoFh+aYuQn8zO8RgYqKUkCvB5c5/3heGkH6C2ii6C5b0ZVCYNBrObBzqehn7tpEu4DfVouF5q/xQQigg==; 31:oepuq40ltmfpfSPhoLgX89cq1zg0VU13sMRXTlm1T+m5LkTiGny3gkWfsGXR3iUywITlTUX2WQjEyZucMlPk2KMCj6h+aVwYMi8fAhRyakgZHgLXVQNyUnl+ieUtpyPadg7b1F0RtGY4fSLmcQ2H60WxlKJxe3WE0Ta6U4FqEP/Q2IGUh3SSN4i2DpnRIeAJ6wLRT+S4HFysK/1hD8tTVfLbjj3ScNGpX8smZlJdEk4= X-MS-TrafficTypeDiagnostic: AM5PR0802MB2481: NoDisclaimer: True X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2481; 20:PoIfepQmh9o9ZR8ORF85X7YQsoyNXg3GUjYc0ntdQKZYIAkBCHtF8XEnoIxBwGaOiFNJX/04Cxkm/ngns+26jcvHKmV1OIIFVAsV1jhywtYJneGFInTTOfo58/gR9pJoG5iDp/yQHb7eukpoifvnKCfOKQVezT8jXocw1+0zbdo=; 4:Es7J0vQEYcV2pZ7wbkUxAw0bkMJHOZHv/uFgSqOS82Dc41Vo4QruiHehPtEokVeNrTsuXaE1gOgTu2rHAA38bVO7hXH5OlEW/HyeKhFDYmLBuJ0JuyDbkad6Jy/6vIs0wPsPBdbfH8Fw8cAdSC6sS8c8v47G4500l8ZI1bKgkon0L+GMk9pjSnHc1/kHE/imUgRcIA/BXr5vkWi11UDQ7n0t01noNKZc3b3Tb4VQXEoJq8TS99Ji5jE/KTmgYEoqk92GFDBNuHRbvFYgIUG5lDRM9BK7Rw/T2pkImrsETKI= X-Exchange-Antispam-Report-Test: UriScan:(180628864354917); X-Microsoft-Antispam-PRVS: X-Exchange-Antispam-Report-CFA-Test: BCL:0; PCL:0; RULEID:(100000700101)(100105000095)(100000701101)(100105300095)(100000702101)(100105100095)(102415395)(6040450)(2401047)(8121501046)(5005006)(93006095)(93001095)(10201501046)(100000703101)(100105400095)(3002001)(6055026)(6041248)(20161123560025)(20161123555025)(20161123558100)(20161123562025)(20161123564025)(201703131423075)(201702281528075)(201703061421075)(201703061406153)(6072148)(201708071742011)(100000704101)(100105200095)(100000705101)(100105500095); SRVR:AM5PR0802MB2481; BCL:0; PCL:0; RULEID:(100000800101)(100110000095)(100000801101)(100110300095)(100000802101)(100110100095)(100000803101)(100110400095)(100000804101)(100110200095)(100000805101)(100110500095); SRVR:AM5PR0802MB2481; X-Forefront-PRVS: 04359FAD81 X-Forefront-Antispam-Report: SFV:NSPM; SFS:(10009020)(4630300001)(7370300001)(6049001)(6009001)(346002)(39860400002)(376002)(377424004)(199003)(189002)(65956001)(65806001)(66066001)(53936002)(4326008)(25786009)(80316001)(105586002)(106356001)(5890100001)(68736007)(33656002)(305945005)(7736002)(81156014)(8676002)(81166006)(72206003)(478600001)(77096006)(6486002)(2476003)(6116002)(3846002)(5660300001)(6916009)(270700001)(4610100001)(568964002)(2950100002)(21490400003)(5000100001)(36756003)(59896002)(189998001)(2906002)(58126008)(101416001)(16576012)(16526017)(50986999)(83506001)(16586007)(76176999)(84326002)(64126003)(86362001)(87266999)(54356999)(7350300001)(316002)(65816999)(97736004); DIR:OUT; SFP:1101; SCL:1; SRVR:AM5PR0802MB2481; H:[10.2.206.69]; FPR:; SPF:None; PTR:InfoNoRecords; A:1; MX:1; LANG:en; Received-SPF: None (protection.outlook.com: arm.com does not designate permitted sender hosts) X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2481; 23:OF1qhrKxkMQjjcFGClkb3r/UEHPOjD938fEn34OC1K5vDwHmABn7bvkwMU5EdSXlmHbur4+ksVv6gqOGHqBSSIIphsWy+e1Nrm/EPTCEnihuSO/W/q9wkmvV9dXv7S5h/UcIILKqQzbNqR+zjiqn33/MAAJxATx6ZMnKEZfeCP86gLR46CPDqEhENWadTt8d5W8SpOxhYN4wtMISJL6+yqB2Y+8trqW1DspvyunqJd2pS8KjtUXMu0gGDu75usHZVztHpfCsun3tFhgwU/2Hw1CG6fgru1QL/ImKh0dVjv8pJv1mRxqsb3TCGdLupDKNUB/mtI0/JTJMGTzs6vnxYc8nk37RG7T/ofVe5tARa25au+KEIg/fg0tE0z4m1I/CzHUm/iIzOdjwtjvR32yeE1DZcEko2XL9TzsJzuWkvL0BPPpIIr5AmwF9LL/Eeh99iTaIdYRQ/mrRnY25YqZ5cNfk6j7PRTTr9N0zqLGSxkEAGSJmoNKXoZCSsoSogUSBpPyVpo/E6e5zpUQ53zjxP6nSOGEcukVaadjDPsPCuuhJKxGNmHAWEDPTp7IhpLweWqHNhzWHAUIk5copmlaD2xaVas4sBcHE1vYmyH3qta+ubpCawihdu6NPmJoWjlgrZz0r97Q2E027zckNHfJKotIa5YA9hC1G+zV233V6kHNw+aeHre3Ezd8NsSd1Ur3RQSmY/7D0TNYvoFJrSR4ImLDIEVig0yfbC9sWH3tZApoRFR8BCqh/Dm9ZgUz5QpPyB662bfy2njsep2HIFg6dapxR0Rn6ikegOq9UHQHhnQr7N9RFpcd2RCGUoE+VzgIlE35Bsm3xzLfDFFHwbhlpiwaoD9FnB2YV21jQYpydtVRoWz/p0CVGvX4z8rEitNOeFcYXc3ewqkv8us+z+kF8ITbLLavtet4t/ZJJFcp+X0z8/oA23NgPgLEONv2ObcSrKS3B50axoZxj1sAoSVh0XS5s1uLaBbhTCtGW4wljzkLPrD457BvmmjTvErbdLzNMHGDrF95DmXgZ5DycDjx3AwesiBOg8OCMHzW1fQJkGBeTBA7brGdXsJ+NRxQfCHY2Ui2d3F+dBtSoUc2gIfRN8jVuSa8Yy1nOJFl2ixSQ3P/1+wdg+gvH3FtNNCzkTb4QMxX6LH4LScb90ozJRn1p2UVyaSaBACZ+VBiwNRddPdNjvstOGpiyNFGMFDMQiPqg9MKM6ZMyQhZE6wITtKiJMLh+AmyC3xIE1BaJuUnCB1RrumtUiVC5qrMsBGag2s2M3jM5CbiY6WtyO0OqY3e/xOiQf7Y+T8EE5Bpw+A6B8Z5/PNO9VX5Ve3ip7V9QvlBvbzv6ThVzUmNHXSRm24EXjGngo2WHSneHFHYIJDJj7VVJnqDcH4fh/Km5/rGqwlAW+89VvtIGJcA4LxAw4kd+itRzMhGLRBBP63Ww+NhZWmiA6GoKLR3glna3bHFoZxe2j0jZeirCv8+HEm+mCygVIA== X-Microsoft-Exchange-Diagnostics: 1; AM5PR0802MB2481; 6:uTqpdjrpyHZk6dlG0cijcG1ZrigZUSnyg0/tmF3mXQowBb5dvzrGrh5NKM29xztzTwAMuJAECY4qFtjbZWZl4ShMe/gWTWJM14HCkABWzfCxtJZaYL5Swmay08OTLgmjUS8FvvYcYDm9e7lqXRlgQ7amMzP515tRIMtMcdpryhUi5C9evz/CUVzfWzndzZyxBW8Tl4XoY0J5lZDgF7nYpdcol91K2vADX9/d1rpuWTLDv+cTcmGE9bTQqLED6iJAWAiauwxv7kbEmqyJGskr4yufGZ2Uh+NOHUJB2lzBlqZsbqvNsLfLG/Dl7nQhC3eFokcjOZOiRitjmwqJyKjqxw==; 5:4l/Hqwor/EdUUBcIn+N7wZJx8epyCgS1Th/hvl2DIFZD7jC82fGTGQ2F4uoYYqnz/JJqESSsHJ/XXaTSCs3PukLFRqiTnil2biaHN4XpaJI5ENuKSFbECTKv0MIbhB2pSRgS1uxECbZzHbVNymmuhw==; 24:Qu3QhmYFuyZ8fPKJnrg+LQpj5hquxvjO7n+C4yddAjwVnUuwVjS+HSHhVOvdBTv+nKJObNXjqqaglx/pTczxL7zGgSW2TZcwt9h+GMAxuic=; 7:mv0JqWBEVV1tn77NT7T05WOWTbPhjVoOtavT5yIPj6PadOG6vvBeDsUYGgnyBuyIytsL9VOgU0Orj78kIisCqFltEfW4VG0/4yepEkzHT012SqeR7BYufc/6yomuux4C8o5iVwNRcYVjm4dITcpp/MKUsZNZYtpBMwWu7IyGLmOQCZf7p9yQQM+AQyrSD6w77O2vIDprgH+dPZQHU9zd/1+3KN1LyjqyMRTyY4IFCyk= SpamDiagnosticOutput: 1:99 SpamDiagnosticMetadata: NSPM X-OriginatorOrg: arm.com X-MS-Exchange-CrossTenant-OriginalArrivalTime: 19 Sep 2017 13:04:58.9664 (UTC) X-MS-Exchange-CrossTenant-FromEntityHeader: Hosted X-MS-Exchange-Transport-CrossTenantHeadersStamped: AM5PR0802MB2481 The new generic logf, log2f and powf code don't need wrappers any more, they set errno inline, so only use the wrappers on targets that need it. 2017-09-19 Szabolcs Nagy * sysdeps/ieee754/flt-32/e_log2f.c (__log2f): Define without wrapper. * sysdeps/ieee754/flt-32/e_logf.c (__logf): Likewise * sysdeps/ieee754/flt-32/e_powf.c (__powf): Likewise * sysdeps/ieee754/flt-32/w_log2f.c: New file. * sysdeps/ieee754/flt-32/w_logf.c: New file. * sysdeps/ieee754/flt-32/w_powf.c: New file. * sysdeps/i386/fpu/w_log2f.c: New file. * sysdeps/i386/fpu/w_logf.c: New file. * sysdeps/i386/fpu/w_powf.c: New file. * sysdeps/m68k/m680x0/fpu/w_log2f.c: New file. * sysdeps/m68k/m680x0/fpu/w_logf.c: New file. * sysdeps/m68k/m680x0/fpu/w_powf.c: New file. diff --git a/sysdeps/i386/fpu/w_log2f.c b/sysdeps/i386/fpu/w_log2f.c new file mode 100644 index 0000000000..3f5c71cec2 --- /dev/null +++ b/sysdeps/i386/fpu/w_log2f.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/i386/fpu/w_logf.c b/sysdeps/i386/fpu/w_logf.c new file mode 100644 index 0000000000..ea48d1356e --- /dev/null +++ b/sysdeps/i386/fpu/w_logf.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/i386/fpu/w_powf.c b/sysdeps/i386/fpu/w_powf.c new file mode 100644 index 0000000000..d133216f5b --- /dev/null +++ b/sysdeps/i386/fpu/w_powf.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/ieee754/flt-32/e_log2f.c b/sysdeps/ieee754/flt-32/e_log2f.c index 6c42f27843..a3506dcfd1 100644 --- a/sysdeps/ieee754/flt-32/e_log2f.c +++ b/sysdeps/ieee754/flt-32/e_log2f.c @@ -18,6 +18,7 @@ #include #include +#include #include "math_config.h" /* @@ -34,7 +35,7 @@ Relative error: 1.9 * 2^-26 (before rounding.) #define OFF 0x3f330000 float -__ieee754_log2f (float x) +__log2f (float x) { /* double_t for better performance on targets with FLT_EVAL_METHOD==2. */ double_t z, r, r2, p, y, y0, invc, logc; @@ -85,4 +86,6 @@ __ieee754_log2f (float x) y = y * r2 + p; return (float) y; } -strong_alias (__ieee754_log2f, __log2f_finite) +strong_alias (__log2f, __ieee754_log2f) +strong_alias (__log2f, __log2f_finite) +versioned_symbol (libm, __log2f, log2f, GLIBC_2_27); diff --git a/sysdeps/ieee754/flt-32/e_logf.c b/sysdeps/ieee754/flt-32/e_logf.c index b8d262441f..ba537d5236 100644 --- a/sysdeps/ieee754/flt-32/e_logf.c +++ b/sysdeps/ieee754/flt-32/e_logf.c @@ -18,6 +18,7 @@ #include #include +#include #include "math_config.h" /* @@ -35,7 +36,7 @@ Relative error: 1.957 * 2^-26 (before rounding.) #define OFF 0x3f330000 float -__ieee754_logf (float x) +__logf (float x) { /* double_t for better performance on targets with FLT_EVAL_METHOD==2. */ double_t z, r, r2, y, y0, invc, logc; @@ -84,4 +85,6 @@ __ieee754_logf (float x) y = y * r2 + (y0 + r); return (float) y; } -strong_alias (__ieee754_logf, __logf_finite) +strong_alias (__logf, __ieee754_logf) +strong_alias (__logf, __logf_finite) +versioned_symbol (libm, __logf, logf, GLIBC_2_27); diff --git a/sysdeps/ieee754/flt-32/e_powf.c b/sysdeps/ieee754/flt-32/e_powf.c index 644a18d05e..4368f49be2 100644 --- a/sysdeps/ieee754/flt-32/e_powf.c +++ b/sysdeps/ieee754/flt-32/e_powf.c @@ -18,6 +18,7 @@ #include #include +#include #include "math_config.h" /* @@ -139,7 +140,7 @@ zeroinfnan (uint32_t ix) } float -__ieee754_powf (float x, float y) +__powf (float x, float y) { unsigned long sign_bias = 0; uint32_t ix, iy; @@ -214,4 +215,6 @@ __ieee754_powf (float x, float y) } return (float) exp2_inline (ylogx, sign_bias); } -strong_alias (__ieee754_powf, __powf_finite) +strong_alias (__powf, __ieee754_powf) +strong_alias (__powf, __powf_finite) +versioned_symbol (libm, __powf, powf, GLIBC_2_27); diff --git a/sysdeps/ieee754/flt-32/w_log2f.c b/sysdeps/ieee754/flt-32/w_log2f.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/ieee754/flt-32/w_log2f.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ieee754/flt-32/w_logf.c b/sysdeps/ieee754/flt-32/w_logf.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/ieee754/flt-32/w_logf.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/ieee754/flt-32/w_powf.c b/sysdeps/ieee754/flt-32/w_powf.c new file mode 100644 index 0000000000..1cc8931700 --- /dev/null +++ b/sysdeps/ieee754/flt-32/w_powf.c @@ -0,0 +1 @@ +/* Not needed. */ diff --git a/sysdeps/m68k/m680x0/fpu/w_log2f.c b/sysdeps/m68k/m680x0/fpu/w_log2f.c new file mode 100644 index 0000000000..3f5c71cec2 --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/w_log2f.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/m68k/m680x0/fpu/w_logf.c b/sysdeps/m68k/m680x0/fpu/w_logf.c new file mode 100644 index 0000000000..ea48d1356e --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/w_logf.c @@ -0,0 +1 @@ +#include diff --git a/sysdeps/m68k/m680x0/fpu/w_powf.c b/sysdeps/m68k/m680x0/fpu/w_powf.c new file mode 100644 index 0000000000..d133216f5b --- /dev/null +++ b/sysdeps/m68k/m680x0/fpu/w_powf.c @@ -0,0 +1 @@ +#include From patchwork Tue Sep 19 13:06:31 2017 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit X-Patchwork-Submitter: Szabolcs Nagy X-Patchwork-Id: 815528 Return-Path: X-Original-To: incoming@patchwork.ozlabs.org Delivered-To: patchwork-incoming@bilbo.ozlabs.org Authentication-Results: ozlabs.org; spf=pass (mailfrom) smtp.mailfrom=sourceware.org (client-ip=209.132.180.131; helo=sourceware.org; envelope-from=libc-alpha-return-84730-incoming=patchwork.ozlabs.org@sourceware.org; receiver=) Authentication-Results: ozlabs.org; dkim=pass (1024-bit key; secure) header.d=sourceware.org header.i=@sourceware.org header.b="WWkYBB0O"; dkim-atps=neutral Received: from sourceware.org (server1.sourceware.org [209.132.180.131]) (using TLSv1.2 with cipher ECDHE-RSA-AES256-GCM-SHA384 (256/256 bits)) (No client certificate requested) by ozlabs.org (Postfix) with ESMTPS id 3xxNSK0TlGz9s4s for ; Tue, 19 Sep 2017 23:06:56 +1000 (AEST) DomainKey-Signature: a=rsa-sha1; c=nofws; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; q=dns; s=default; b= jY3gUELtIbw8rxUBssIC0m2KHCbJi5nbStYl7+SeRXr/+W8YiEWYKNi22fokR02b H9A5lSw6cUoJGpyoIbEo9iGmB4wY2jhyRgQJmFbD5hwNRR71KLdP3k7OBThT3iUm KGwrq2qrSCnE6mo9bSg44Q0TnQZyJnzFop4grrmqV1U= DKIM-Signature: v=1; a=rsa-sha1; c=relaxed; d=sourceware.org; h=list-id :list-unsubscribe:list-subscribe:list-archive:list-post :list-help:sender:message-id:date:from:mime-version:to:cc :subject:references:in-reply-to:content-type; s=default; bh=Y/uW qaFaxoVLt0tj+k5Qpeq/7nc=; b=WWkYBB0ONl1DmJroekZrCk00e54MUykj1Kme /T0DEZikS+IK+1YnnhhNUcMDa2dMWgrIYpabUG9R3wBsfKiGgpTIHN8AHwF+YYvq KkxT2JhOToG4KAqLw2Ql2I75az7UBg6s5xBW/1blBu2vShmOymV2VtD9w/mgrQc5 aRbdRp0= Received: (qmail 72153 invoked by alias); 19 Sep 2017 13:06:49 -0000 Mailing-List: contact libc-alpha-help@sourceware.org; run by ezmlm Precedence: bulk List-Id: List-Unsubscribe: List-Subscribe: List-Archive: List-Post: List-Help: , Sender: libc-alpha-owner@sourceware.org Delivered-To: mailing list libc-alpha@sourceware.org Received: (qmail 72128 invoked by uid 89); 19 Sep 2017 13:06:48 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-24.7 required=5.0 tests=AWL, BAYES_00, GIT_PATCH_0, GIT_PATCH_1, GIT_PATCH_2, GIT_PATCH_3, RCVD_IN_DNSWL_NONE, SPF_HELO_PASS, SPF_PASS autolearn=ham version=3.3.2 spammy= X-HELO: EUR01-HE1-obe.outbound.protection.outlook.com Authentication-Results: spf=none (sender IP is ) smtp.mailfrom=Szabolcs.Nagy@arm.com; Message-ID: <59C11657.30400@arm.com> Date: Tue, 19 Sep 2017 14:06:31 +0100 From: Szabolcs Nagy User-Agent: Mozilla/5.0 (X11; Linux x86_64; rv:31.0) Gecko/20100101 Thunderbird/31.8.0 MIME-Version: 1.0 To: GNU C Library CC: nd@arm.com Subject: [PATCH 9/9] AArch64: update libm-test-ulps References: <59C1123F.9080003@arm.com> In-Reply-To: <59C1123F.9080003@arm.com> X-ClientProxiedBy: DB6PR0202CA0048.eurprd02.prod.outlook.com (2603:10a6:4:a5::34) To DB6PR0802MB2488.eurprd08.prod.outlook.com (2603:10a6:4:a0::23) X-MS-PublicTrafficType: Email X-MS-Office365-Filtering-Correlation-Id: be0bd749-26fe-41f6-39f0-08d4ff5f411c X-MS-Office365-Filtering-HT: Tenant X-Microsoft-Antispam: UriScan:; BCL:0; PCL:0; RULEID:(300000500095)(300135000095)(300000501095)(300135300095)(22001)(300000502095)(300135100095)(2017030254152)(48565401081)(300000503095)(300135400095)(49563074)(201703131423075)(201703031133081)(201702281549075)(300000504095)(300135200095)(300000505095)(300135600095)(300000506095)(300135500095); SRVR:DB6PR0802MB2488; X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 3:AHWgM+tgkV09B2o3T/4GNAvmYb0pxnEiqD96PCthDesNYK+2zK/bwfcc9hdOHDHS5k7Z9p/344c1fbEue/+Vf50Z7cIJSGIAubDVzrgby4eNP/UQur+PMdmdG0E2qbchFmqKGm5zD2ZcJfttQYJZt92NsSargusXH4UBah6fkLTQr/b0k2+Tte+y2XhDCbdKijMhdBSkOFHMktRJlsQ6MGtJyntrScFJgYFFKKun8eVijF/UXJDZonHagE+nPZZv; 25:+6JgNtglDA1T5MS2Bm1kPIzVtrws4eE8itgpmEUHBsupPwOb4tY1y4QP3EL/daMI2KVFpLDGReVkRIdLoNRaDO3ECk+AFIbJdXfH4OhdEhLy9Ln7y6PFGBCTT6xVtPMdy1IFNAiBovOp5sOO1nHm/TlRVFqsGYELsOt9SgOVLdz4w40sWpEPM+Q8Ckk7Jfh492XrIcoa2ahrALuw6JzCp9n3MwX31YYvWEo/yRYzfDmy0xSgKvIu38StT6Z8fwicusByB0jDCnwBLXPorV92hkkN/NK/i4+7A+VXdT2lcUtG7cpJLbX3UVddrsKueq7vlhM/SeDisl2vi3MpnNwR0w==; 31:G2zFcVhk7gM3TBcnYE2UBRDHDsMb0xHq4sp6oY18glEmq0TIMu7xp1gKe2d645fqR+v+gIQv0vQKhWDQ/EeAv/W4/PNnk1JZjMEi0kkNqffmn9gO6QiivJ6IIEuzPLb86zFVHBf3cEKy4ZJ09IRj9EXxl8YJwtaIBmS8+0mJ6HDEtfeCrdmxDIVcKB5biCoXDsVO4n8SehBT5aRjTOe6QdnuoF99gnLJAVLN3D8YNaE= X-MS-TrafficTypeDiagnostic: DB6PR0802MB2488: NoDisclaimer: True X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 20:BJrd9yw8hiKdjqWR0gNalwVQ6XCNLJTxLamDa08FgDvkBAqU8ZXN1yvVq+/RIKY2hwCAu2QLHUL/exrzG0Pcc3NMHaEJ3WRhxF3QTKcb7zH7ohezxc2KjeAS9uvXyVO2NvFTJnVG2ho2HXcSkkxVVNUI7cqJRvZde+frFB8Vg0Q=; 4:khQEEJqfS9bxtFL7MNQ8ngeV+lB+uAKEbZHZM1F+SVpoC1RxNupNQWcAp62OHWy0XfEVBSRZESIYkU+X77oM6lnwjZUe59qza/9Ld05RU0TKgO85YHzgmg0v8y7YDFm/9cxj3iiYUAuS1/W1+RGYnrOH9FWxfmYNjwxHi/z7Cw6cM5BLp61nZx2Nm3lAHUYnlNa+3YqQM6HNngX2271uAbnlSSbLgIYG9VdGgsy4x2PmypYH5E2G1qS1yfCrEYnqsi0GP1+yizUbvZCriHMpW5aa3udKH9UldI+m6HlUmco= X-Exchange-Antispam-Report-Test: UriScan:(180628864354917); X-Microsoft-Antispam-PRVS: X-Exchange-Antispam-Report-CFA-Test: BCL:0; PCL:0; RULEID:(100000700101)(100105000095)(100000701101)(100105300095)(100000702101)(100105100095)(102415395)(6040450)(2401047)(8121501046)(5005006)(100000703101)(100105400095)(10201501046)(3002001)(93006095)(93001095)(6055026)(6041248)(20161123562025)(20161123564025)(20161123555025)(201703131423075)(201702281528075)(201703061421075)(201703061406153)(20161123560025)(20161123558100)(6072148)(201708071742011)(100000704101)(100105200095)(100000705101)(100105500095); SRVR:DB6PR0802MB2488; BCL:0; PCL:0; RULEID:(100000800101)(100110000095)(100000801101)(100110300095)(100000802101)(100110100095)(100000803101)(100110400095)(100000804101)(100110200095)(100000805101)(100110500095); SRVR:DB6PR0802MB2488; X-Forefront-PRVS: 04359FAD81 X-Forefront-Antispam-Report: SFV:NSPM; SFS:(10009020)(4630300001)(7370300001)(979002)(6009001)(6049001)(376002)(346002)(189002)(199003)(377424004)(36756003)(5660300001)(2950100002)(33656002)(66066001)(305945005)(7736002)(77096006)(6486002)(7350300001)(6666003)(6916009)(8676002)(81166006)(86362001)(81156014)(21490400003)(53936002)(16526017)(16576012)(58126008)(64126003)(97736004)(16586007)(65956001)(65806001)(316002)(478600001)(105586002)(106356001)(25786009)(4326008)(5000100001)(101416001)(83506001)(15650500001)(50986999)(564344004)(568964002)(270700001)(2476003)(4610100001)(230783001)(5890100001)(68736007)(54356999)(76176999)(2906002)(65816999)(189998001)(84326002)(6116002)(72206003)(3846002)(969003)(989001)(999001)(1009001)(1019001); DIR:OUT; SFP:1101; SCL:1; SRVR:DB6PR0802MB2488; H:[10.2.206.69]; FPR:; SPF:None; PTR:InfoNoRecords; MX:3; A:1; LANG:en; Received-SPF: None (protection.outlook.com: arm.com does not designate permitted sender hosts) X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 23:v+mvpuAkzKm5A9wSJkgZ6ye5lmW/C4TuYgPGCi9Qia/l+YqoksXXyyu1r1t1pMy4IYZr3BpoFaQDLxDgalJ1G9vvcAsAp7enSvDIibLRfOFDle8IKhKYKqYFMKw40HnUdhESwM39OdvGH6ZaR6OGJwV3vGES/EoFBl2xao672Wcsx70qOCIhNpXTG/+VKBDm0uq5CahMgz+Gs8DjWO+yLmdjv9GPAQoyh85lb3uuhcouSxRQaeiHydlHERyfdtTAl129JxkOdZhN1nBgEbtM77q2wN8Mkxq6jpYzaa8xxdydfvKkbQSciOram5PGKpjmjcgiDWboQ20E2N5hK5FFhw3Ko/q5kMbIQIF4jfOXi37MNkD+yn1/5ymYEsK/1qAY0+qTbK9xkkGVvDdKUNY++qahxjOKUCDi9iWox2gw1ispZgqDhUhT1DhyAiNQnvm3ye48+hQXWC8LY2VLeL55nkwrc7OTFMgELN29QoVIDgzUHEVt5xHnHiKAWLFDi8uUqfbPqOoKoQCVS+EDSMq5xhuP0JU4/TvrH4+o1sEGIEnWbP7cNxtZv14YtXOvJbkPbMrQbNM9MX7WdD1AnWGbfxivmUpLzMIwTHVi0MfUtqRyTr96Pv3X7J3Lv/oNcuBSsYV2X0KfkCydFOCkMKxq1GYZiDBs75qOkvmG+xISY6qo96pYMCnIg7KTak0VohQWvFMf+3YeBB7QXYjocHqpWruS9JUqRl2GS05aufu79JP+ySR6Jrr38/gH1Z+IyK2t2dvspW+3jvavyZ8hsAjze9wCnCe6Fia1E0PRa5ySru1LcYOxbRbwtl3C242phQlG09pgRFld7lO5utq2S17LqzFdDO6C+G6mxBpeDXDhi8urbwK8ZNIVj3DELAR/zLuc0SSlU1kc3tUomWtveK6R0bwytWMw6X5JRkiPiO4N6MwGywkvYQbA2p6J+7KybQ1JmCX4+yiEDfPfMhlzt/x/hGEpmwqEtF30wP77wVa37PlcMXTGQ7iNLxmQW2PR07MW9nSp7H05KXPWFjFG2Wg8HseH9k7QNWNs+9/coAmJMAqmMtlhshkvT34Yt3K3lCeJUMw6oQP9k2xCF+aEs26Yav7VpNqcLYu8y9tKZnxpQhI9s9PQFk+C8FhgyrwVVM0UA8QPodTwOiLQNLY+RfH74vE7XtfRHSWHA3mlaunsjIKExT5HReBbrwtcLFVbZFUdGftCWRx4yk6vi72Jx0Pa7m5JsY6mavw7r0uIRsp+KnaDS0T29RkFRtiSM8eJzRVSu/kTRaQf1iqX6DbPVGzFOu4MWKMsERN/UTY0tM3xANCP8y+LJSY76mHuPnrUzZeTRM6ZWr8ZcJSyjIv4iEweBkn6q4v0DLnlnqjyqOM9BQvpkR8yDfY8sAdYuCRTwAwCUMNoHAYjXjDz/YWLe9Qs9RpFL/MykrJ5I0kS2dLLqSFXEtqvBOOPhoggPiwqPRzBFTiBkh2FfgG5wzb+PdCUHesR5kSLs+KKG8MuyKjK+5PyZZxaCyws8JaCf/fsuWe1DMtkOEgdCEJep+WbBOcoLIxIk/JN/H1Qg7SFXmYO4Xd4PkWsWaBahH9yyPfF5e9a X-Microsoft-Exchange-Diagnostics: 1; DB6PR0802MB2488; 6:rCyGqnvEINKx9aShtpjvCg+3siw+0VsS7L20J59lj7lH+MTjbXBNyK9wDUjDk+6xT3L3sl5nnrMqCg4WPd/R7unWRgQi1+cMagRy+HnMIvv9Y51D9xzPF5w3J4GRQwmhTDXbZ4o92XRBx7fWCsCdWHgLQxyL/238w4EzmBMVSEQUbr+7lS9DYyjQpY4O3MxIFxOM6qcZvPF0xs0arEhSgE2AeG3PnMk2b3a/2zxIhtWe4HZD5n25fz6KtZd0HU6pwa0Ou9luGZ9JsCtu0e66DP0YUjlw8QR77/XV4rm4KVXcgFgw83qtxxjG4fE/uyKqA0T2e1lulp1UyWnIjXsU2A==; 5:eYhgwHdAJ76ExrcQEUi8m8/Iqe2Fzw8Lg84bMw4ca61PqxHsEfnNEsR5oMFUUClbzBEO623NYhjNuqcg0TS4XGgTzZSQAXdPTo4rQ5g2kCmi4S8jF1/9wjCyKPizUKM1nyu5FdoAAcUqry1OObSc6g==; 24:59PXVRibduFPBYMxtZjAs04uwzPnoVp4BoW7rRv2AqTTO6MLTypSOFTFb+vgJDdopGCR8LxP9p/6RkSQ/S3/4qSuwu0vyy/qVAXlffXcGAM=; 7:Nb6yM4wccux4/HkGu8u6l17uCfHMlmX8eddbcmdr/EJbo1rWntCjQaXyg7HLQiKTz3RSTRJH08wGaDCkmsI4VabVnVEMdZ3hSZtAB+q16kPNIbChIuqZlm0rwz+3Mv0yevOwDMQw2cWzIANG1BRQdeYqlPGHPFLpEPbTFrrakEnXCaOw8MSk8/w8L8/mdZJIFlxhRhPEtwPH6f0jLh8r/+8nMf0sgAjjTnh+lKiwghg= SpamDiagnosticOutput: 1:99 SpamDiagnosticMetadata: NSPM X-OriginatorOrg: arm.com X-MS-Exchange-CrossTenant-OriginalArrivalTime: 19 Sep 2017 13:06:34.1677 (UTC) X-MS-Exchange-CrossTenant-FromEntityHeader: Hosted X-MS-Exchange-Transport-CrossTenantHeadersStamped: DB6PR0802MB2488 Update for new expf and logf. 2017-09-19 Szabolcs Nagy * sysdeps/aarch64/libm-test-ulps: Update. diff --git a/sysdeps/aarch64/libm-test-ulps b/sysdeps/aarch64/libm-test-ulps index 797a8ef8cb..234fc6c489 100644 --- a/sysdeps/aarch64/libm-test-ulps +++ b/sysdeps/aarch64/libm-test-ulps @@ -917,9 +917,9 @@ ldouble: 2 Function: Real part of "clog10_downward": double: 5 -float: 4 +float: 5 idouble: 5 -ifloat: 4 +ifloat: 5 ildouble: 3 ldouble: 3 @@ -1583,15 +1583,21 @@ ldouble: 2 Function: "exp_downward": double: 1 +float: 1 idouble: 1 +ifloat: 1 Function: "exp_towardzero": double: 1 +float: 1 idouble: 1 +ifloat: 1 Function: "exp_upward": double: 1 +float: 1 idouble: 1 +ifloat: 1 Function: "expm1": double: 1