From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: by sourceware.org (Postfix, from userid 1791) id C4BDA3858C78; Wed, 27 Dec 2023 13:23:58 +0000 (GMT) DKIM-Filter: OpenDKIM Filter v2.11.0 sourceware.org C4BDA3858C78 DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=sourceware.org; s=default; t=1703683438; bh=r2HM01hlhUjhCIyJ+CO6PHaqyay6kaQhEFYwBre7j64=; h=From:To:Subject:Date:From; b=Siuxrj5VJ39sDRXPr/tTvlw1WiGWzHF/GM2UnY+FnOqqgfEcKiDCzuYRhsUPaGBPd 7FxdoueSv43CWsvkaxfWiXe2EiiTJ57fOOqWMvPNWJJE+oMlFjWKRVhqaqcmfCe5zF zYke1ZAt3cmU1gar+ryDFqC0inRidseFls+ZOpV0= Content-Type: text/plain; charset="us-ascii" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit From: Adhemerval Zanella To: glibc-cvs@sourceware.org Subject: [glibc/azanella/mips-hw-fp-round] mips: Implement ceilf with hardware floating-point rounding instruction X-Act-Checkin: glibc X-Git-Author: Adhemerval Zanella X-Git-Refname: refs/heads/azanella/mips-hw-fp-round X-Git-Oldrev: 78768e9002314aa427986fe6fb7aad14df877c5e X-Git-Newrev: ac1ef8bd961d5225e890741f15c10966172cbdf4 Message-Id: <20231227132358.C4BDA3858C78@sourceware.org> Date: Wed, 27 Dec 2023 13:23:58 +0000 (GMT) List-Id: https://sourceware.org/git/gitweb.cgi?p=glibc.git;h=ac1ef8bd961d5225e890741f15c10966172cbdf4 commit ac1ef8bd961d5225e890741f15c10966172cbdf4 Author: Adhemerval Zanella Date: Tue Dec 26 19:25:41 2023 -0300 mips: Implement ceilf with hardware floating-point rounding instruction Diff: --- sysdeps/mips/fpu/roundf_to_integer.h | 68 ++++++++++++++++++++++++++++++++++++ sysdeps/mips/fpu/s_ceilf.c | 36 +++++++++++++++++++ 2 files changed, 104 insertions(+) diff --git a/sysdeps/mips/fpu/roundf_to_integer.h b/sysdeps/mips/fpu/roundf_to_integer.h new file mode 100644 index 0000000000..49f5960f24 --- /dev/null +++ b/sysdeps/mips/fpu/roundf_to_integer.h @@ -0,0 +1,68 @@ +/* Round to integer generic implementation. + Copyright (C) 2023 Free Software . + 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 _ROUND_TO_INTEGER_H +#define _ROUND_TO_INTEGER_H + +#include +#include +#include + +enum round_mode +{ + CEIL, +}; + +static inline float +round_to_integer_float (enum round_mode mode, float x) +{ + uint32_t hx = asuint (x); + int ex = (hx & ~SIGN_MASK) >> MANTISSA_WIDTH; + + float r = x; + + fenv_t fe; + libc_feholdexceptf (&fe); + switch (mode) + { + case CEIL: + asm ("ceil.l.s %0, %0" : "+f" (r)); + break; + } + libc_fesetenvf (&fe); + + /* Inf or NaN. */ + if (__glibc_unlikely (ex == 0xff)) + return x + x; + + /* Number input and no fraction, return the number itself. */ + if (ex >= 127 + 25) + return x; + + asm ("cvt.s.l %0, %1" : "+f" (r)); + /* The copysign seems to generate better code. */ +#if 1 + return copysignf (r, x); +#else + if (ex >= 127) + return r; + return hx & SIGN_MASK ? -r : r; +#endif +} + +#endif diff --git a/sysdeps/mips/fpu/s_ceilf.c b/sysdeps/mips/fpu/s_ceilf.c new file mode 100644 index 0000000000..84335efffd --- /dev/null +++ b/sysdeps/mips/fpu/s_ceilf.c @@ -0,0 +1,36 @@ +/* Round to nearest integer, away from zero. MIPS version. + Copyright (C) 2023 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 + . */ + +#if ((__mips_fpr == 64) \ + && (__mips_hard_float == 1) \ + && ((__mips == 32 && __mips_isa_rev > 1) || __mips == 64)) + +#define NO_MATH_REDIRECT +#include +#include +#include + +float +__ceilf (float x) +{ + return round_to_integer_float (CEIL, x); +} +libm_alias_float (__ceil, ceil) +#else +# include +#endif