From 628d90c5f97b6f0f8b79a079b682febb1f486a38 Mon Sep 17 00:00:00 2001 From: Vineet Gupta Date: Fri, 29 May 2020 15:58:33 -0700 Subject: ieee754: provide gcc builtins based generic fma functions Reviewed-by: Adhemerval Zanella --- sysdeps/ieee754/ldbl-128/s_fma.c | 5 +++++ 1 file changed, 5 insertions(+) (limited to 'sysdeps/ieee754/ldbl-128/s_fma.c') diff --git a/sysdeps/ieee754/ldbl-128/s_fma.c b/sysdeps/ieee754/ldbl-128/s_fma.c index 13f8ed3b87..0feff0a71c 100644 --- a/sysdeps/ieee754/ldbl-128/s_fma.c +++ b/sysdeps/ieee754/ldbl-128/s_fma.c @@ -21,6 +21,7 @@ #include #include #include +#include /* This implementation relies on long double being more than twice as precise as double and uses rounding to odd in order to avoid problems @@ -31,6 +32,9 @@ double __fma (double x, double y, double z) { +#if USE_FMA_BUILTIN + return __builtin_fma (x, y, z); +#else fenv_t env; /* Multiplication is always exact. */ long double temp = (long double) x * (long double) y; @@ -50,6 +54,7 @@ __fma (double x, double y, double z) feupdateenv (&env); /* And finally truncation with round to nearest. */ return (double) u.d; +#endif /* ! USE_FMA_BUILTIN */ } #ifndef __fma libm_alias_double (__fma, fma) -- cgit 1.4.1