diff options
Diffstat (limited to 'sysdeps/ieee754/soft-fp')
| -rw-r--r-- | sysdeps/ieee754/soft-fp/Makefile | 3 | ||||
| -rw-r--r-- | sysdeps/ieee754/soft-fp/s_fma.c | 68 | ||||
| -rw-r--r-- | sysdeps/ieee754/soft-fp/s_fmaf.c | 68 | ||||
| -rw-r--r-- | sysdeps/ieee754/soft-fp/s_fmal.c | 66 |
4 files changed, 205 insertions, 0 deletions
diff --git a/sysdeps/ieee754/soft-fp/Makefile b/sysdeps/ieee754/soft-fp/Makefile new file mode 100644 index 0000000000..ada13e8b70 --- /dev/null +++ b/sysdeps/ieee754/soft-fp/Makefile @@ -0,0 +1,3 @@ +ifeq ($(subdir),math) +CPPFLAGS += -I../soft-fp +endif diff --git a/sysdeps/ieee754/soft-fp/s_fma.c b/sysdeps/ieee754/soft-fp/s_fma.c new file mode 100644 index 0000000000..313547d061 --- /dev/null +++ b/sysdeps/ieee754/soft-fp/s_fma.c @@ -0,0 +1,68 @@ +/* Implement fma using soft-fp. + Copyright (C) 2013-2017 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. + + In addition to the permissions in the GNU Lesser General Public + License, the Free Software Foundation gives you unlimited + permission to link the compiled version of this file into + combinations with other programs, and to distribute those + combinations without any restriction coming from the use of this + file. (The Lesser General Public License restrictions do apply in + other respects; for example, they cover modification of the file, + and distribution when not linked into a combine executable.) + + 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 + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <libc-diag.h> +#include <libm-alias-double.h> + +/* R_e is not set in cases where it is not used in packing, but the + compiler does not see that it is set in all cases where it is + used, resulting in warnings that it may be used uninitialized. + The location of the warning differs in different versions of GCC, + it may be where R is defined using a macro or it may be where the + macro is defined. */ +DIAG_PUSH_NEEDS_COMMENT; +DIAG_IGNORE_NEEDS_COMMENT (4.9, "-Wmaybe-uninitialized"); + +#include "soft-fp.h" +#include "double.h" + +double +__fma (double a, double b, double c) +{ + FP_DECL_EX; + FP_DECL_D (A); + FP_DECL_D (B); + FP_DECL_D (C); + FP_DECL_D (R); + double r; + + FP_INIT_ROUNDMODE; + FP_UNPACK_D (A, a); + FP_UNPACK_D (B, b); + FP_UNPACK_D (C, c); + FP_FMA_D (R, A, B, C); + FP_PACK_D (r, R); + FP_HANDLE_EXCEPTIONS; + + return r; +} +DIAG_POP_NEEDS_COMMENT; + +#ifndef __fma +libm_alias_double (__fma, fma) +#endif diff --git a/sysdeps/ieee754/soft-fp/s_fmaf.c b/sysdeps/ieee754/soft-fp/s_fmaf.c new file mode 100644 index 0000000000..295223da50 --- /dev/null +++ b/sysdeps/ieee754/soft-fp/s_fmaf.c @@ -0,0 +1,68 @@ +/* Implement fmaf using soft-fp. + Copyright (C) 2013-2017 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. + + In addition to the permissions in the GNU Lesser General Public + License, the Free Software Foundation gives you unlimited + permission to link the compiled version of this file into + combinations with other programs, and to distribute those + combinations without any restriction coming from the use of this + file. (The Lesser General Public License restrictions do apply in + other respects; for example, they cover modification of the file, + and distribution when not linked into a combine executable.) + + 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 + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <libc-diag.h> +#include <libm-alias-float.h> + +/* R_e is not set in cases where it is not used in packing, but the + compiler does not see that it is set in all cases where it is + used, resulting in warnings that it may be used uninitialized. + The location of the warning differs in different versions of GCC, + it may be where R is defined using a macro or it may be where the + macro is defined. */ +DIAG_PUSH_NEEDS_COMMENT; +DIAG_IGNORE_NEEDS_COMMENT (4.9, "-Wmaybe-uninitialized"); + +#include "soft-fp.h" +#include "single.h" + +float +__fmaf (float a, float b, float c) +{ + FP_DECL_EX; + FP_DECL_S (A); + FP_DECL_S (B); + FP_DECL_S (C); + FP_DECL_S (R); + float r; + + FP_INIT_ROUNDMODE; + FP_UNPACK_S (A, a); + FP_UNPACK_S (B, b); + FP_UNPACK_S (C, c); + FP_FMA_S (R, A, B, C); + FP_PACK_S (r, R); + FP_HANDLE_EXCEPTIONS; + + return r; +} +DIAG_POP_NEEDS_COMMENT; + +#ifndef __fmaf +libm_alias_float (__fma, fma) +#endif diff --git a/sysdeps/ieee754/soft-fp/s_fmal.c b/sysdeps/ieee754/soft-fp/s_fmal.c new file mode 100644 index 0000000000..79be52c4f1 --- /dev/null +++ b/sysdeps/ieee754/soft-fp/s_fmal.c @@ -0,0 +1,66 @@ +/* Implement fmal using soft-fp. + Copyright (C) 2013-2017 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. + + In addition to the permissions in the GNU Lesser General Public + License, the Free Software Foundation gives you unlimited + permission to link the compiled version of this file into + combinations with other programs, and to distribute those + combinations without any restriction coming from the use of this + file. (The Lesser General Public License restrictions do apply in + other respects; for example, they cover modification of the file, + and distribution when not linked into a combine executable.) + + 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 + <http://www.gnu.org/licenses/>. */ + +#include <math.h> +#include <libc-diag.h> +#include <libm-alias-ldouble.h> + +/* R_e is not set in cases where it is not used in packing, but the + compiler does not see that it is set in all cases where it is + used, resulting in warnings that it may be used uninitialized. + The location of the warning differs in different versions of GCC, + it may be where R is defined using a macro or it may be where the + macro is defined. */ +DIAG_PUSH_NEEDS_COMMENT; +DIAG_IGNORE_NEEDS_COMMENT (4.9, "-Wmaybe-uninitialized"); + +#include "soft-fp.h" +#include "quad.h" + +long double +__fmal (long double a, long double b, long double c) +{ + FP_DECL_EX; + FP_DECL_Q (A); + FP_DECL_Q (B); + FP_DECL_Q (C); + FP_DECL_Q (R); + long double r; + + FP_INIT_ROUNDMODE; + FP_UNPACK_Q (A, a); + FP_UNPACK_Q (B, b); + FP_UNPACK_Q (C, c); + FP_FMA_Q (R, A, B, C); + FP_PACK_Q (r, R); + FP_HANDLE_EXCEPTIONS; + + return r; +} +DIAG_POP_NEEDS_COMMENT; + +libm_alias_ldouble (__fma, fma) |