1 //===-- Unittests for fma ------------------------------------------------===//
2 //
3 // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
4 // See https://llvm.org/LICENSE.txt for license information.
5 // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
6 //
7 //===----------------------------------------------------------------------===//
8 
9 #include "FmaTest.h"
10 
11 #include "src/math/fma.h"
12 
13 struct Inputs {
14   double a, b, c;
15 };
16 
17 struct LlvmLibcFmaTest : public FmaTestTemplate<double> {
18   void test_more_values() {
19     constexpr int N = 236;
20     constexpr Inputs INPUTS[N] = {
21         {0x1p+0, 0x2p+0, 0x3p+0},
22         {0x1.4p+0, 0xcp-4, 0x1p-4},
23         {0x0p+0, 0x0p+0, 0x0p+0},
24         {0x1p+0, 0x0p+0, 0x0p+0},
25         {0x0p+0, 0x1p+0, 0x0p+0},
26         {0x1p+0, 0x1p+0, 0x1p+0},
27         {0x0p+0, 0x0p+0, 0x1p+0},
28         {0x0p+0, 0x0p+0, 0x2p+0},
29         {0x0p+0, 0x0p+0, 0xf.fffffp+124},
30         {0x0p+0, 0x0p+0, 0xf.ffffffffffff8p+1020},
31         {0x0p+0, 0x1p+0, 0x1p+0},
32         {0x1p+0, 0x0p+0, 0x1p+0},
33         {0x0p+0, 0x1p+0, 0x2p+0},
34         {0x1p+0, 0x0p+0, 0x2p+0},
35         {0x0p+0, 0x1p+0, 0xf.fffffp+124},
36         {0x0p+0, 0x1p+0, 0xf.ffffffffffff8p+1020},
37         {0x1p+0, 0x0p+0, 0xf.fffffp+124},
38         {0x1p+0, 0x0p+0, 0xf.ffffffffffff8p+1020},
39         {0x4p-128, 0x4p-128, 0x0p+0},
40         {0x4p-128, 0x4p-1024, 0x0p+0},
41         {0x4p-128, 0x8p-972, 0x0p+0},
42         {0x4p-1024, 0x4p-128, 0x0p+0},
43         {0x4p-1024, 0x4p-1024, 0x0p+0},
44         {0x4p-1024, 0x8p-972, 0x0p+0},
45         {0x8p-972, 0x4p-128, 0x0p+0},
46         {0x8p-972, 0x4p-1024, 0x0p+0},
47         {0x8p-972, 0x8p-972, 0x0p+0},
48         {0x4p-128, 0x4p-128, 0x0p+0},
49         {0x4p-128, 0x4p-1024, 0x0p+0},
50         {0x4p-128, 0x8p-972, 0x0p+0},
51         {0x4p-1024, 0x4p-128, 0x0p+0},
52         {0x4p-1024, 0x4p-1024, 0x0p+0},
53         {0x4p-1024, 0x8p-972, 0x0p+0},
54         {0x8p-972, 0x4p-128, 0x0p+0},
55         {0x8p-972, 0x4p-1024, 0x0p+0},
56         {0x8p-972, 0x8p-972, 0x0p+0},
57         {0x4p-128, 0x4p-128, 0x0p+0},
58         {0x4p-128, 0x4p-1024, 0x0p+0},
59         {0x4p-128, 0x8p-972, 0x0p+0},
60         {0x4p-1024, 0x4p-128, 0x0p+0},
61         {0x4p-1024, 0x4p-1024, 0x0p+0},
62         {0x4p-1024, 0x8p-972, 0x0p+0},
63         {0x8p-972, 0x4p-128, 0x0p+0},
64         {0x8p-972, 0x4p-1024, 0x0p+0},
65         {0x8p-972, 0x8p-972, 0x0p+0},
66         {0x4p-128, 0x4p-128, 0x0p+0},
67         {0x4p-128, 0x4p-1024, 0x0p+0},
68         {0x4p-128, 0x8p-972, 0x0p+0},
69         {0x4p-1024, 0x4p-128, 0x0p+0},
70         {0x4p-1024, 0x4p-1024, 0x0p+0},
71         {0x4p-1024, 0x8p-972, 0x0p+0},
72         {0x8p-972, 0x4p-128, 0x0p+0},
73         {0x8p-972, 0x4p-1024, 0x0p+0},
74         {0x8p-972, 0x8p-972, 0x0p+0},
75         {0x4p-128, 0x4p-128, 0x0p+0},
76         {0x4p-128, 0x4p-1024, 0x0p+0},
77         {0x4p-128, 0x8p-972, 0x0p+0},
78         {0x4p-1024, 0x4p-128, 0x0p+0},
79         {0x4p-1024, 0x4p-1024, 0x0p+0},
80         {0x4p-1024, 0x8p-972, 0x0p+0},
81         {0x8p-972, 0x4p-128, 0x0p+0},
82         {0x8p-972, 0x4p-1024, 0x0p+0},
83         {0x8p-972, 0x8p-972, 0x0p+0},
84         {0x4p-128, 0x4p-128, 0x0p+0},
85         {0x4p-128, 0x4p-1024, 0x0p+0},
86         {0x4p-128, 0x8p-972, 0x0p+0},
87         {0x4p-1024, 0x4p-128, 0x0p+0},
88         {0x4p-1024, 0x4p-1024, 0x0p+0},
89         {0x4p-1024, 0x8p-972, 0x0p+0},
90         {0x8p-972, 0x4p-128, 0x0p+0},
91         {0x8p-972, 0x4p-1024, 0x0p+0},
92         {0x8p-972, 0x8p-972, 0x0p+0},
93         {0x4p-128, 0x4p-128, 0x0p+0},
94         {0x4p-128, 0x4p-1024, 0x0p+0},
95         {0x4p-128, 0x8p-972, 0x0p+0},
96         {0x4p-1024, 0x4p-128, 0x0p+0},
97         {0x4p-1024, 0x4p-1024, 0x0p+0},
98         {0x4p-1024, 0x8p-972, 0x0p+0},
99         {0x8p-972, 0x4p-128, 0x0p+0},
100         {0x8p-972, 0x4p-1024, 0x0p+0},
101         {0x8p-972, 0x8p-972, 0x0p+0},
102         {0x4p-128, 0x4p-128, 0x0p+0},
103         {0x4p-128, 0x4p-1024, 0x0p+0},
104         {0x4p-128, 0x8p-972, 0x0p+0},
105         {0x4p-1024, 0x4p-128, 0x0p+0},
106         {0x4p-1024, 0x4p-1024, 0x0p+0},
107         {0x4p-1024, 0x8p-972, 0x0p+0},
108         {0x8p-972, 0x4p-128, 0x0p+0},
109         {0x8p-972, 0x4p-1024, 0x0p+0},
110         {0x8p-972, 0x8p-972, 0x0p+0},
111         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
112         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
113         {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
114         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
115         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
116         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
117         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
118         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
119         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
120         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
121         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
122         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
123         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
124         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
125         {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
126         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
127         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
128         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
129         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
130         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
131         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
132         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
133         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
134         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
135         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
136         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
137         {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
138         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
139         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
140         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
141         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
142         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
143         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
144         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
145         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
146         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
147         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
148         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
149         {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
150         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
151         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
152         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
153         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
154         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
155         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
156         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
157         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
158         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
159         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
160         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
161         {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
162         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
163         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
164         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
165         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
166         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
167         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
168         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
169         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
170         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
171         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
172         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
173         {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
174         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
175         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
176         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
177         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
178         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
179         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
180         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
181         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
182         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
183         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
184         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
185         {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
186         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
187         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
188         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
189         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
190         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
191         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
192         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
193         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
194         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
195         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
196         {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
197         {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
198         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
199         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
200         {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
201         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
202         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
203         {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
204         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
205         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
206         {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
207         {0x2.fffp+12, 0x1.000002p+0, 0x1.ffffp-24},
208         {0x1.fffp+0, 0x1.00001p+0, 0x1.fffp+0},
209         {0xc.d5e6fp+124, 0x2.6af378p-128, 0x1.f08948p+0},
210         {0x1.9abcdep+100, 0x2.6af378p-128, 0x3.e1129p-28},
211         {0xf.fffffp+124, 0x1.001p+0, 0xf.fffffp+124},
212         {0xf.fffffp+124, 0x1.fffffep+0, 0xf.fffffp+124},
213         {0xf.fffffp+124, 0x2p+0, 0xf.fffffp+124},
214         {0x5p-128, 0x8.00002p-4, 0x1p-128},
215         {0x7.ffffep-128, 0x8.00001p-4, 0x8p-152},
216         {0x8p-152, 0x8p-4, 0x3.fffff8p-128},
217         {0x8p-152, 0x8.8p-4, 0x3.fffff8p-128},
218         {0x8p-152, 0x8p-152, 0x8p+124},
219         {0x8p-152, 0x8p-152, 0x4p-128},
220         {0x8p-152, 0x8p-152, 0x3.fffff8p-128},
221         {0x8p-152, 0x8p-152, 0x8p-152},
222         {0xf.ffp-4, 0xf.ffp-4, 0xf.fep-4},
223         {0x4.000008p-128, 0x4.000008p-28, 0x8p+124},
224         {0x4.000008p-128, 0x4.000008p-28, 0x8p+100},
225         {0x2.fep+12, 0x1.0000000000001p+0, 0x1.ffep-48},
226         {0x1.fffp+0, 0x1.0000000000001p+0, 0x1.fffp+0},
227         {0x1.0000002p+0, 0xf.fffffep-4, 0x1p-300},
228         {0xe.f56df7797f768p+1020, 0x3.7ab6fbbcbfbb4p-1024,
229          0x3.40bf1803497f6p+0},
230         {0x1.deadbeef2feedp+900, 0x3.7ab6fbbcbfbb4p-1024,
231          0x6.817e300692fecp-124},
232         {0xf.ffffffffffff8p+1020, 0x1.001p+0, 0xf.ffffffffffff8p+1020},
233         {0xf.ffffffffffff8p+1020, 0x1.fffffffffffffp+0,
234          0xf.ffffffffffff8p+1020},
235         {0xf.ffffffffffff8p+1020, 0x2p+0, 0xf.ffffffffffff8p+1020},
236         {0x5.a827999fcef3p-540, 0x5.a827999fcef3p-540, 0x0p+0},
237         {0x3.bd5b7dde5fddap-496, 0x3.bd5b7dde5fddap-496, 0xd.fc352bc352bap-992},
238         {0x3.bd5b7dde5fddap-504, 0x3.bd5b7dde5fddap-504,
239          0xd.fc352bc352bap-1008},
240         {0x8p-540, 0x4p-540, 0x4p-1076},
241         {0x1.7fffff8p-968, 0x4p-108, 0x4p-1048},
242         {0x2.8000008p-968, 0x4p-108, 0x4p-1048},
243         {0x2.8p-968, 0x4p-108, 0x4p-1048},
244         {0x2.33956cdae7c2ep-960, 0x3.8e211518bfea2p-108,
245          0x2.02c2b59766d9p-1024},
246         {0x3.a5d5dadd1d3a6p-980, 0x2.9c0cd8c5593bap-64, 0x2.49179ac00d15p-1024},
247         {0x2.2a7aca1773e0cp-908, 0x9.6809186a42038p-128, 0x2.c9e356b3f0fp-1024},
248         {0x3.ffffffffffffep-712, 0x3.ffffffffffffep-276,
249          0x3.fffffc0000ffep-984},
250         {0x5p-1024, 0x8.000000000001p-4, 0x1p-1024},
251         {0x7.ffffffffffffp-1024, 0x8.0000000000008p-4, 0x4p-1076},
252         {0x4p-1076, 0x8p-4, 0x3.ffffffffffffcp-1024},
253         {0x4p-1076, 0x8.8p-4, 0x3.ffffffffffffcp-1024},
254         {0x4p-1076, 0x4p-1076, 0x8p+1020},
255         {0x4p-1076, 0x4p-1076, 0x4p-1024},
256         {0x4p-1076, 0x4p-1076, 0x3.ffffffffffffcp-1024},
257         {0x4p-1076, 0x4p-1076, 0x4p-1076},
258         {0xf.ffffffffffff8p-4, 0xf.ffffffffffff8p-4, 0xf.ffffffffffffp-4},
259         {0x4.0000000000004p-1024, 0x2.0000000000002p-56, 0x8p+1020},
260         {0x4.0000000000004p-1024, 0x2.0000000000002p-56, 0x4p+968},
261         {0x7.fffff8p-128, 0x3.fffffcp+24, 0xf.fffffp+124},
262         {0x7.ffffffffffffcp-1024, 0x7.ffffffffffffcp+52,
263          0xf.ffffffffffff8p+1020},
264     };
265 
266     for (int i = 0; i < N; ++i) {
267       for (int signs = 0; signs < 7; ++signs) {
268         double a = (signs & 4) ? -INPUTS[i].a : INPUTS[i].a;
269         double b = (signs & 2) ? -INPUTS[i].b : INPUTS[i].b;
270         double c = (signs & 1) ? -INPUTS[i].c : INPUTS[i].c;
271         mpfr::TernaryInput<double> input{a, b, c};
272         ASSERT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Fma, input,
273                                        __llvm_libc::fma(a, b, c), 0.5);
274       }
275     }
276   }
277 };
278 
279 TEST_F(LlvmLibcFmaTest, SpecialNumbers) {
280   test_special_numbers(&__llvm_libc::fma);
281 }
282 
283 TEST_F(LlvmLibcFmaTest, SubnormalRange) {
284   test_subnormal_range(&__llvm_libc::fma);
285 }
286 
287 TEST_F(LlvmLibcFmaTest, NormalRange) { test_normal_range(&__llvm_libc::fma); }
288 
289 TEST_F(LlvmLibcFmaTest, ExtraValues) { test_more_values(); }
290