From: "H.J. Lu" <hjl.tools@gmail.com>
To: Adhemerval Zanella <adhemerval.zanella@linaro.org>
Cc: libc-alpha@sourceware.org, Joseph Myers <josmyers@redhat.com>,
Florian Weimer <fweimer@redhat.com>
Subject: Re: [PATCH v2 06/10] i386: Use generic exp10
Date: Wed, 27 Mar 2024 13:14:32 -0700 [thread overview]
Message-ID: <CAMe9rOo4FQeyjPphhti45rd3wmB5+wD4EaMaTVN+bYHGBAobhg@mail.gmail.com> (raw)
In-Reply-To: <20240327194024.1409677-7-adhemerval.zanella@linaro.org>
On Wed, Mar 27, 2024 at 12:40 PM Adhemerval Zanella
<adhemerval.zanella@linaro.org> wrote:
>
> The resulting performance is slight better (Ryzen 5900, gcc 13.2.1):
>
> * master
> "exp10": {
> "": {
> "duration": 3.70091e+09,
> "iterations": 5.8534e+07,
> "max": 91.279,
> "min": 62.6225,
> "mean": 63.2267
> }
> }
>
> * patch
> "exp10": {
> "": {
> "duration": 3.70793e+09,
> "iterations": 6.328e+07,
> "max": 259.592,
> "min": 52.1145,
> "mean": 58.5957
> }
> }
>
> Checked on i686-linux-gnu.
> ---
> sysdeps/i386/fpu/Versions | 1 +
> sysdeps/i386/fpu/e_exp10.S | 51 -----------------------
> sysdeps/i386/fpu/e_exp10.c | 2 +
> sysdeps/i386/fpu/e_exp_data.c | 1 -
> sysdeps/i386/fpu/w_exp10_compat.c | 8 ----
> sysdeps/ieee754/dbl-64/e_exp10.c | 7 +++-
> sysdeps/mach/hurd/i386/libm.abilist | 1 +
> sysdeps/unix/sysv/linux/i386/libm.abilist | 1 +
> 8 files changed, 10 insertions(+), 62 deletions(-)
> delete mode 100644 sysdeps/i386/fpu/e_exp10.S
> create mode 100644 sysdeps/i386/fpu/e_exp10.c
> delete mode 100644 sysdeps/i386/fpu/e_exp_data.c
> delete mode 100644 sysdeps/i386/fpu/w_exp10_compat.c
>
> diff --git a/sysdeps/i386/fpu/Versions b/sysdeps/i386/fpu/Versions
> index 9509f9b7c7..7326f25583 100644
> --- a/sysdeps/i386/fpu/Versions
> +++ b/sysdeps/i386/fpu/Versions
> @@ -5,6 +5,7 @@ libm {
> }
> GLIBC_2.40 {
> # No SVID compatible error handling.
> + exp10;
> fmod; fmodf;
> }
> }
> diff --git a/sysdeps/i386/fpu/e_exp10.S b/sysdeps/i386/fpu/e_exp10.S
> deleted file mode 100644
> index 902f70b77f..0000000000
> --- a/sysdeps/i386/fpu/e_exp10.S
> +++ /dev/null
> @@ -1,51 +0,0 @@
> -
> -#include <machine/asm.h>
> -#include <i386-math-asm.h>
> -#include <libm-alias-finite.h>
> -
> -DEFINE_DBL_MIN
> -
> -#ifdef PIC
> -# define MO(op) op##@GOTOFF(%ecx)
> -#else
> -# define MO(op) op
> -#endif
> -
> - .text
> -/* 10^x = 2^(x * log2(10)) */
> -ENTRY(__ieee754_exp10)
> -#ifdef PIC
> - LOAD_PIC_REG (cx)
> -#endif
> - fldl 4(%esp)
> -/* I added the following ugly construct because exp(+-Inf) resulted
> - in NaN. The ugliness results from the bright minds at Intel.
> - For the i686 the code can be written better.
> - -- drepper@cygnus.com. */
> - fxam /* Is NaN or +-Inf? */
> - fstsw %ax
> - movb $0x45, %dh
> - andb %ah, %dh
> - cmpb $0x05, %dh
> - je 1f /* Is +-Inf, jump. */
> - fldl2t
> - fmulp /* x * log2(10) */
> - fld %st
> - frndint /* int(x * log2(10)) */
> - fsubr %st,%st(1) /* fract(x * log2(10)) */
> - fxch
> - f2xm1 /* 2^(fract(x * log2(10))) - 1 */
> - fld1
> - faddp /* 2^(fract(x * log2(10))) */
> - fscale /* e^x */
> - fstp %st(1)
> - DBL_NARROW_EVAL_UFLOW_NONNEG_NAN
> - ret
> -
> -1: testl $0x200, %eax /* Test sign. */
> - jz 2f /* If positive, jump. */
> - fstp %st
> - fldz /* Set result to 0. */
> -2: ret
> -END (__ieee754_exp10)
> -libm_alias_finite (__ieee754_exp10, __exp10)
> diff --git a/sysdeps/i386/fpu/e_exp10.c b/sysdeps/i386/fpu/e_exp10.c
> new file mode 100644
> index 0000000000..340254fc6e
> --- /dev/null
> +++ b/sysdeps/i386/fpu/e_exp10.c
> @@ -0,0 +1,2 @@
> +#define EXP10_VERSION GLIBC_2_40
> +#include <sysdeps/ieee754/dbl-64/e_exp10.c>
> diff --git a/sysdeps/i386/fpu/e_exp_data.c b/sysdeps/i386/fpu/e_exp_data.c
> deleted file mode 100644
> index 1cc8931700..0000000000
> --- a/sysdeps/i386/fpu/e_exp_data.c
> +++ /dev/null
> @@ -1 +0,0 @@
> -/* Not needed. */
> diff --git a/sysdeps/i386/fpu/w_exp10_compat.c b/sysdeps/i386/fpu/w_exp10_compat.c
> deleted file mode 100644
> index 49a0e03385..0000000000
> --- a/sysdeps/i386/fpu/w_exp10_compat.c
> +++ /dev/null
> @@ -1,8 +0,0 @@
> -/* i386 provides an optimized __ieee754_exp10. */
> -#ifdef SHARED
> -# define NO_COMPAT_NEEDED 1
> -# include <math/w_exp10_compat.c>
> -#else
> -# include <math-type-macros-double.h>
> -# include <w_exp10_template.c>
> -#endif
> diff --git a/sysdeps/ieee754/dbl-64/e_exp10.c b/sysdeps/ieee754/dbl-64/e_exp10.c
> index 225fc74c4c..c63b852f72 100644
> --- a/sysdeps/ieee754/dbl-64/e_exp10.c
> +++ b/sysdeps/ieee754/dbl-64/e_exp10.c
> @@ -99,7 +99,7 @@ __exp10 (double x)
>
> /* Reduce x: z = x * N / log10(2), k = round(z). */
> double_t z = __exp_data.invlog10_2N * x;
> - double_t kd;
> + double kd;
> int64_t ki;
> #if TOINT_INTRINSICS
> kd = roundtoint (z);
> @@ -147,7 +147,10 @@ __exp10 (double x)
> strong_alias (__exp10, __ieee754_exp10)
> libm_alias_finite (__ieee754_exp10, __exp10)
> #if LIBM_SVID_COMPAT
> -versioned_symbol (libm, __exp10, exp10, GLIBC_2_39);
> +# ifndef EXP10_VERSION
> +# define EXP10_VERSION GLIBC_2_39
> +# endif
> +versioned_symbol (libm, __exp10, exp10, EXP10_VERSION);
> libm_alias_double_other (__exp10, exp10)
> #else
> libm_alias_double (__exp10, exp10)
> diff --git a/sysdeps/mach/hurd/i386/libm.abilist b/sysdeps/mach/hurd/i386/libm.abilist
> index 88e7538e51..01c5633663 100644
> --- a/sysdeps/mach/hurd/i386/libm.abilist
> +++ b/sysdeps/mach/hurd/i386/libm.abilist
> @@ -1181,5 +1181,6 @@ GLIBC_2.35 fsqrt F
> GLIBC_2.35 fsqrtl F
> GLIBC_2.35 hypot F
> GLIBC_2.35 hypotf F
> +GLIBC_2.40 exp10 F
> GLIBC_2.40 fmod F
> GLIBC_2.40 fmodf F
> diff --git a/sysdeps/unix/sysv/linux/i386/libm.abilist b/sysdeps/unix/sysv/linux/i386/libm.abilist
> index c99c60161d..3413cfdbe7 100644
> --- a/sysdeps/unix/sysv/linux/i386/libm.abilist
> +++ b/sysdeps/unix/sysv/linux/i386/libm.abilist
> @@ -1188,5 +1188,6 @@ GLIBC_2.35 fsqrt F
> GLIBC_2.35 fsqrtl F
> GLIBC_2.35 hypot F
> GLIBC_2.35 hypotf F
> +GLIBC_2.40 exp10 F
> GLIBC_2.40 fmod F
> GLIBC_2.40 fmodf F
> --
> 2.34.1
>
Also need a bug report.
--
H.J.
next prev parent reply other threads:[~2024-03-27 20:15 UTC|newest]
Thread overview: 27+ messages / expand[flat|nested] mbox.gz Atom feed top
2024-03-27 19:40 [PATCH v2 00/10] Fix some libm static issues Adhemerval Zanella
2024-03-27 19:40 ` [PATCH v2 01/10] math: Add support for auto static math tests Adhemerval Zanella
2024-03-27 19:40 ` [PATCH v2 02/10] math: Fix i386 and m68k fmod/fmodf on static build (BZ 31488) Adhemerval Zanella
2024-03-27 19:40 ` [PATCH v2 03/10] i386: Use generic fmod Adhemerval Zanella
2024-03-27 19:55 ` H.J. Lu
2024-03-27 20:37 ` Adhemerval Zanella Netto
2024-03-27 21:38 ` H.J. Lu
2024-03-28 14:11 ` Adhemerval Zanella Netto
2024-03-28 14:51 ` H.J. Lu
2024-03-28 15:14 ` Adhemerval Zanella Netto
[not found] ` <CAMe9rOqhQDA-zk=+oTvdoPpq=rGEhtan0couaZ3Z_fxeFpa=7A@mail.gmail.com>
[not found] ` <9bf7af32-afa3-439c-84c0-76e76b220e44@linaro.org>
[not found] ` <CAMe9rOqTcmRc9mvguQnDRFb=BTjZJ5CptxjtPjtBGoB-a5mc6Q@mail.gmail.com>
[not found] ` <0e4733a3-d569-4a73-b2d3-001cc6c3d751@linaro.org>
2024-03-28 16:00 ` H.J. Lu
2024-03-28 18:22 ` Adhemerval Zanella Netto
2024-03-28 18:38 ` Joseph Myers
2024-03-28 19:37 ` Adhemerval Zanella Netto
2024-03-28 19:57 ` H.J. Lu
2024-03-27 19:40 ` [PATCH v2 04/10] i386: Use generic fmodf Adhemerval Zanella
2024-03-27 19:55 ` H.J. Lu
2024-03-27 19:40 ` [PATCH v2 05/10] math: Fix i386 and m68k exp10 on static build Adhemerval Zanella
2024-03-27 19:57 ` H.J. Lu
2024-03-27 20:39 ` Adhemerval Zanella Netto
2024-03-27 20:55 ` Joseph Myers
2024-03-27 19:40 ` [PATCH v2 06/10] i386: Use generic exp10 Adhemerval Zanella
2024-03-27 20:14 ` H.J. Lu [this message]
2024-03-27 19:40 ` [PATCH v2 07/10] math: Fix isnanf128 static build Adhemerval Zanella
2024-03-27 19:40 ` [PATCH v2 08/10] math: Provided copysignf128 for static libm on alpha, s390, and sparcv9 Adhemerval Zanella
2024-03-27 19:40 ` [PATCH v2 09/10] math: Provide frexpf128 " Adhemerval Zanella
2024-03-27 19:40 ` [PATCH v2 10/10] math: Provide modf128 " Adhemerval Zanella
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
List information: https://www.gnu.org/software/libc/involved.html
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=CAMe9rOo4FQeyjPphhti45rd3wmB5+wD4EaMaTVN+bYHGBAobhg@mail.gmail.com \
--to=hjl.tools@gmail.com \
--cc=adhemerval.zanella@linaro.org \
--cc=fweimer@redhat.com \
--cc=josmyers@redhat.com \
--cc=libc-alpha@sourceware.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).