diff options
Diffstat (limited to 'sysdeps/ieee754/ldbl-128')
-rw-r--r-- | sysdeps/ieee754/ldbl-128/s_dfmal.c | 38 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/s_f64xfmaf128.c | 39 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/s_ffmal.c | 34 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/s_fma.c | 6 | ||||
-rw-r--r-- | sysdeps/ieee754/ldbl-128/s_fmal.c | 4 |
5 files changed, 121 insertions, 0 deletions
diff --git a/sysdeps/ieee754/ldbl-128/s_dfmal.c b/sysdeps/ieee754/ldbl-128/s_dfmal.c new file mode 100644 index 0000000..e00ddba --- /dev/null +++ b/sysdeps/ieee754/ldbl-128/s_dfmal.c @@ -0,0 +1,38 @@ +/* Fused multiply-add of long double (ldbl-128) value, narrowing the result + to double. + Copyright (C) 2021 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 + <https://www.gnu.org/licenses/>. */ + +#define f32xfmaf64x __hide_f32xfmaf64x +#define f32xfmaf128 __hide_f32xfmaf128 +#define f64fmaf64x __hide_f64fmaf64x +#define f64fmaf128 __hide_f64fmaf128 +#include <math.h> +#undef f32xfmaf64x +#undef f32xfmaf128 +#undef f64fmaf64x +#undef f64fmaf128 + +#include <math-narrow.h> + +double +__dfmal (_Float128 x, _Float128 y, _Float128 z) +{ + NARROW_FMA_ROUND_TO_ODD (x, y, z, double, union ieee854_long_double, l, + mantissa3, false); +} +libm_alias_double_ldouble (fma) diff --git a/sysdeps/ieee754/ldbl-128/s_f64xfmaf128.c b/sysdeps/ieee754/ldbl-128/s_f64xfmaf128.c new file mode 100644 index 0000000..8616451 --- /dev/null +++ b/sysdeps/ieee754/ldbl-128/s_f64xfmaf128.c @@ -0,0 +1,39 @@ +/* Fused multiply-add of _Float128 value, converting the result to _Float64x. + Copyright (C) 2021 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 + <https://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <math-narrow.h> +#include <tininess.h> + +/* math_ldbl.h defines _Float128 to long double for this directory, + but when they are different, this function must be defined with + _Float128 arguments to avoid defining an alias with an incompatible + type. */ +#undef _Float128 + +#if __HAVE_FLOAT64X_LONG_DOUBLE && __HAVE_DISTINCT_FLOAT128 +_Float64x +__f64xfmaf128 (_Float128 x, _Float128 y, _Float128 z) +{ + NARROW_FMA_ROUND_TO_ODD (x, y, z, _Float64x, union ieee854_long_double, l, + mantissa3, TININESS_AFTER_ROUNDING); +} +libm_alias_float64x_float128 (fma) +#else +/* Defined as an alias of fmal. */ +#endif diff --git a/sysdeps/ieee754/ldbl-128/s_ffmal.c b/sysdeps/ieee754/ldbl-128/s_ffmal.c new file mode 100644 index 0000000..65d5e48 --- /dev/null +++ b/sysdeps/ieee754/ldbl-128/s_ffmal.c @@ -0,0 +1,34 @@ +/* Fused multiply-add of long double (ldbl-128) value, narrowing the result + to float. + Copyright (C) 2021 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 + <https://www.gnu.org/licenses/>. */ + +#define f32fmaf64x __hide_f32fmaf64x +#define f32fmaf128 __hide_f32fmaf128 +#include <math.h> +#undef f32fmaf64x +#undef f32fmaf128 + +#include <math-narrow.h> + +float +__ffmal (_Float128 x, _Float128 y, _Float128 z) +{ + NARROW_FMA_ROUND_TO_ODD (x, y, z, float, union ieee854_long_double, l, + mantissa3, false); +} +libm_alias_float_ldouble (fma) diff --git a/sysdeps/ieee754/ldbl-128/s_fma.c b/sysdeps/ieee754/ldbl-128/s_fma.c index 4795e71..bc90c88 100644 --- a/sysdeps/ieee754/ldbl-128/s_fma.c +++ b/sysdeps/ieee754/ldbl-128/s_fma.c @@ -17,10 +17,15 @@ <https://www.gnu.org/licenses/>. */ #define NO_MATH_REDIRECT +#define dfmal __hide_dfmal +#define f32xfmaf64 __hide_f32xfmaf64 #include <math.h> +#undef dfmal +#undef f32xfmaf64 #include <fenv.h> #include <ieee754.h> #include <libm-alias-double.h> +#include <math-narrow-alias.h> #include <math-use-builtins.h> /* This implementation relies on long double being more than twice as @@ -58,4 +63,5 @@ __fma (double x, double y, double z) } #ifndef __fma libm_alias_double (__fma, fma) +libm_alias_double_narrow (__fma, fma) #endif diff --git a/sysdeps/ieee754/ldbl-128/s_fmal.c b/sysdeps/ieee754/ldbl-128/s_fmal.c index aff9efc..4595ad2 100644 --- a/sysdeps/ieee754/ldbl-128/s_fmal.c +++ b/sysdeps/ieee754/ldbl-128/s_fmal.c @@ -18,12 +18,15 @@ #define NO_MATH_REDIRECT #include <float.h> +#define f64xfmaf128 __hide_f64xfmaf128 #include <math.h> +#undef f64xfmaf128 #include <fenv.h> #include <ieee754.h> #include <math-barriers.h> #include <math_private.h> #include <libm-alias-ldouble.h> +#include <math-narrow-alias.h> #include <tininess.h> #include <math-use-builtins.h> @@ -303,3 +306,4 @@ __fmal (_Float128 x, _Float128 y, _Float128 z) #endif /* ! USE_FMAL_BUILTIN */ } libm_alias_ldouble (__fma, fma) +libm_alias_ldouble_narrow (__fma, fma) |