1 | /* Implement fma using soft-fp. |
2 | Copyright (C) 2013-2022 Free Software Foundation, Inc. |
3 | This file is part of the GNU C Library. |
4 | |
5 | The GNU C Library is free software; you can redistribute it and/or |
6 | modify it under the terms of the GNU Lesser General Public |
7 | License as published by the Free Software Foundation; either |
8 | version 2.1 of the License, or (at your option) any later version. |
9 | |
10 | In addition to the permissions in the GNU Lesser General Public |
11 | License, the Free Software Foundation gives you unlimited |
12 | permission to link the compiled version of this file into |
13 | combinations with other programs, and to distribute those |
14 | combinations without any restriction coming from the use of this |
15 | file. (The Lesser General Public License restrictions do apply in |
16 | other respects; for example, they cover modification of the file, |
17 | and distribution when not linked into a combine executable.) |
18 | |
19 | The GNU C Library is distributed in the hope that it will be useful, |
20 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
21 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
22 | Lesser General Public License for more details. |
23 | |
24 | You should have received a copy of the GNU Lesser General Public |
25 | License along with the GNU C Library; if not, see |
26 | <https://www.gnu.org/licenses/>. */ |
27 | |
28 | #define NO_MATH_REDIRECT |
29 | #define dfmal __hide_dfmal |
30 | #define f32xfmaf64 __hide_f32xfmaf64 |
31 | #include <math.h> |
32 | #undef dfmal |
33 | #undef f32xfmaf64 |
34 | #include <libc-diag.h> |
35 | #include <libm-alias-double.h> |
36 | #include <math-narrow-alias.h> |
37 | |
38 | /* R_e is not set in cases where it is not used in packing, but the |
39 | compiler does not see that it is set in all cases where it is |
40 | used, resulting in warnings that it may be used uninitialized. |
41 | The location of the warning differs in different versions of GCC, |
42 | it may be where R is defined using a macro or it may be where the |
43 | macro is defined. */ |
44 | DIAG_PUSH_NEEDS_COMMENT; |
45 | DIAG_IGNORE_NEEDS_COMMENT (4.9, "-Wmaybe-uninitialized" ); |
46 | |
47 | #include "soft-fp.h" |
48 | #include "double.h" |
49 | |
50 | double |
51 | __fma (double a, double b, double c) |
52 | { |
53 | FP_DECL_EX; |
54 | FP_DECL_D (A); |
55 | FP_DECL_D (B); |
56 | FP_DECL_D (C); |
57 | FP_DECL_D (R); |
58 | double r; |
59 | |
60 | FP_INIT_ROUNDMODE; |
61 | FP_UNPACK_D (A, a); |
62 | FP_UNPACK_D (B, b); |
63 | FP_UNPACK_D (C, c); |
64 | FP_FMA_D (R, A, B, C); |
65 | FP_PACK_D (r, R); |
66 | FP_HANDLE_EXCEPTIONS; |
67 | |
68 | return r; |
69 | } |
70 | DIAG_POP_NEEDS_COMMENT; |
71 | |
72 | #ifndef __fma |
73 | libm_alias_double (__fma, fma) |
74 | libm_alias_double_narrow (__fma, fma) |
75 | #endif |
76 | |