From patchwork Wed Apr 15 12:16:49 2015 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit X-Patchwork-Submitter: Wilco Dijkstra X-Patchwork-Id: 6227 Received: (qmail 45990 invoked by alias); 15 Apr 2015 12:16:57 -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 45978 invoked by uid 89); 15 Apr 2015 12:16:57 -0000 Authentication-Results: sourceware.org; auth=none X-Virus-Found: No X-Spam-SWARE-Status: No, score=-1.9 required=5.0 tests=AWL, BAYES_00, SPF_PASS autolearn=ham version=3.3.2 X-HELO: eu-smtp-delivery-143.mimecast.com From: "Wilco Dijkstra" To: Subject: [PATCH][AArch64] Inline __ieee754_sqrt(f) Date: Wed, 15 Apr 2015 13:16:49 +0100 Message-ID: <000101d07776$0c9c4ac0$25d4e040$@com> MIME-Version: 1.0 X-MC-Unique: PuxB3H-AT0SJSy9P2XNuEA-1 Inline __ieee754_sqrt and __ieee754_sqrtf. Also add external definitions. OK for commit? 2015-04-15 Wilco Dijkstra * sysdeps/aarch64/fpu/math_private.h (__ieee754_sqrt): New function. (__ieee754_sqrtf): New function. * sysdeps/aarch64/fpu/e_sqrt.c (__ieee754_sqrt): New function. * sysdeps/aarch64/fpu/e_sqrtf.c (__ieee754_sqrtf): New function. --- sysdeps/aarch64/fpu/e_sqrt.c | 28 ++++++++++++++++++++++++++++ sysdeps/aarch64/fpu/e_sqrtf.c | 28 ++++++++++++++++++++++++++++ sysdeps/aarch64/fpu/math_private.h | 16 ++++++++++++++++ 3 files changed, 72 insertions(+) create mode 100644 sysdeps/aarch64/fpu/e_sqrt.c create mode 100644 sysdeps/aarch64/fpu/e_sqrtf.c diff --git a/sysdeps/aarch64/fpu/e_sqrt.c b/sysdeps/aarch64/fpu/e_sqrt.c new file mode 100644 index 0000000..4f11ca2 --- /dev/null +++ b/sysdeps/aarch64/fpu/e_sqrt.c @@ -0,0 +1,28 @@ +/* Square root of floating point number. + Copyright (C) 2015 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 + +double +__ieee754_sqrt (double d) +{ + double res; + asm ("fsqrt %d0, %d1" : "=w" (res) : "w" (d)); + return res; +} +strong_alias (__ieee754_sqrt, __sqrt_finite) diff --git a/sysdeps/aarch64/fpu/e_sqrtf.c b/sysdeps/aarch64/fpu/e_sqrtf.c new file mode 100644 index 0000000..a2e99e1 --- /dev/null +++ b/sysdeps/aarch64/fpu/e_sqrtf.c @@ -0,0 +1,28 @@ +/* Single-precision floating point square root. + Copyright (C) 2015 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 + +float +__ieee754_sqrtf (float s) +{ + float res; + asm ("fsqrt %s0, %s1" : "=w" (res) : "w" (s)); + return res; +} +strong_alias (__ieee754_sqrtf, __sqrtf_finite) diff --git a/sysdeps/aarch64/fpu/math_private.h b/sysdeps/aarch64/fpu/math_private.h index 52a6ad9..b3c2509 100644 --- a/sysdeps/aarch64/fpu/math_private.h +++ b/sysdeps/aarch64/fpu/math_private.h @@ -22,6 +22,22 @@ #include #include +extern __always_inline double +__ieee754_sqrt (double d) +{ + double res; + asm __volatile__ ("fsqrt %d0, %d1" : "=w" (res) : "w" (d)); + return res; +} + +extern __always_inline float +__ieee754_sqrtf (float s) +{ + float res; + asm __volatile__ ("fsqrt %s0, %s1" : "=w" (res) : "w" (s)); + return res; +} + static __always_inline void libc_feholdexcept_aarch64 (fenv_t *envp) {