aboutsummaryrefslogtreecommitdiff
path: root/libc/test/src/math/fma_test.cpp
blob: dd761382631d5e7f77a3a24aa4abb56bda89331c (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
//===-- Unittests for fma ------------------------------------------------===//
//
// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
// See https://llvm.org/LICENSE.txt for license information.
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
//
//===----------------------------------------------------------------------===//

#include "FmaTest.h"

#include "src/math/fma.h"

struct Inputs {
  double a, b, c;
};

struct LlvmLibcFmaTest : public FmaTestTemplate<double> {
  void test_more_values() {
    constexpr int N = 236;
    constexpr Inputs INPUTS[N] = {
        {0x1p+0, 0x2p+0, 0x3p+0},
        {0x1.4p+0, 0xcp-4, 0x1p-4},
        {0x0p+0, 0x0p+0, 0x0p+0},
        {0x1p+0, 0x0p+0, 0x0p+0},
        {0x0p+0, 0x1p+0, 0x0p+0},
        {0x1p+0, 0x1p+0, 0x1p+0},
        {0x0p+0, 0x0p+0, 0x1p+0},
        {0x0p+0, 0x0p+0, 0x2p+0},
        {0x0p+0, 0x0p+0, 0xf.fffffp+124},
        {0x0p+0, 0x0p+0, 0xf.ffffffffffff8p+1020},
        {0x0p+0, 0x1p+0, 0x1p+0},
        {0x1p+0, 0x0p+0, 0x1p+0},
        {0x0p+0, 0x1p+0, 0x2p+0},
        {0x1p+0, 0x0p+0, 0x2p+0},
        {0x0p+0, 0x1p+0, 0xf.fffffp+124},
        {0x0p+0, 0x1p+0, 0xf.ffffffffffff8p+1020},
        {0x1p+0, 0x0p+0, 0xf.fffffp+124},
        {0x1p+0, 0x0p+0, 0xf.ffffffffffff8p+1020},
        {0x4p-128, 0x4p-128, 0x0p+0},
        {0x4p-128, 0x4p-1024, 0x0p+0},
        {0x4p-128, 0x8p-972, 0x0p+0},
        {0x4p-1024, 0x4p-128, 0x0p+0},
        {0x4p-1024, 0x4p-1024, 0x0p+0},
        {0x4p-1024, 0x8p-972, 0x0p+0},
        {0x8p-972, 0x4p-128, 0x0p+0},
        {0x8p-972, 0x4p-1024, 0x0p+0},
        {0x8p-972, 0x8p-972, 0x0p+0},
        {0x4p-128, 0x4p-128, 0x0p+0},
        {0x4p-128, 0x4p-1024, 0x0p+0},
        {0x4p-128, 0x8p-972, 0x0p+0},
        {0x4p-1024, 0x4p-128, 0x0p+0},
        {0x4p-1024, 0x4p-1024, 0x0p+0},
        {0x4p-1024, 0x8p-972, 0x0p+0},
        {0x8p-972, 0x4p-128, 0x0p+0},
        {0x8p-972, 0x4p-1024, 0x0p+0},
        {0x8p-972, 0x8p-972, 0x0p+0},
        {0x4p-128, 0x4p-128, 0x0p+0},
        {0x4p-128, 0x4p-1024, 0x0p+0},
        {0x4p-128, 0x8p-972, 0x0p+0},
        {0x4p-1024, 0x4p-128, 0x0p+0},
        {0x4p-1024, 0x4p-1024, 0x0p+0},
        {0x4p-1024, 0x8p-972, 0x0p+0},
        {0x8p-972, 0x4p-128, 0x0p+0},
        {0x8p-972, 0x4p-1024, 0x0p+0},
        {0x8p-972, 0x8p-972, 0x0p+0},
        {0x4p-128, 0x4p-128, 0x0p+0},
        {0x4p-128, 0x4p-1024, 0x0p+0},
        {0x4p-128, 0x8p-972, 0x0p+0},
        {0x4p-1024, 0x4p-128, 0x0p+0},
        {0x4p-1024, 0x4p-1024, 0x0p+0},
        {0x4p-1024, 0x8p-972, 0x0p+0},
        {0x8p-972, 0x4p-128, 0x0p+0},
        {0x8p-972, 0x4p-1024, 0x0p+0},
        {0x8p-972, 0x8p-972, 0x0p+0},
        {0x4p-128, 0x4p-128, 0x0p+0},
        {0x4p-128, 0x4p-1024, 0x0p+0},
        {0x4p-128, 0x8p-972, 0x0p+0},
        {0x4p-1024, 0x4p-128, 0x0p+0},
        {0x4p-1024, 0x4p-1024, 0x0p+0},
        {0x4p-1024, 0x8p-972, 0x0p+0},
        {0x8p-972, 0x4p-128, 0x0p+0},
        {0x8p-972, 0x4p-1024, 0x0p+0},
        {0x8p-972, 0x8p-972, 0x0p+0},
        {0x4p-128, 0x4p-128, 0x0p+0},
        {0x4p-128, 0x4p-1024, 0x0p+0},
        {0x4p-128, 0x8p-972, 0x0p+0},
        {0x4p-1024, 0x4p-128, 0x0p+0},
        {0x4p-1024, 0x4p-1024, 0x0p+0},
        {0x4p-1024, 0x8p-972, 0x0p+0},
        {0x8p-972, 0x4p-128, 0x0p+0},
        {0x8p-972, 0x4p-1024, 0x0p+0},
        {0x8p-972, 0x8p-972, 0x0p+0},
        {0x4p-128, 0x4p-128, 0x0p+0},
        {0x4p-128, 0x4p-1024, 0x0p+0},
        {0x4p-128, 0x8p-972, 0x0p+0},
        {0x4p-1024, 0x4p-128, 0x0p+0},
        {0x4p-1024, 0x4p-1024, 0x0p+0},
        {0x4p-1024, 0x8p-972, 0x0p+0},
        {0x8p-972, 0x4p-128, 0x0p+0},
        {0x8p-972, 0x4p-1024, 0x0p+0},
        {0x8p-972, 0x8p-972, 0x0p+0},
        {0x4p-128, 0x4p-128, 0x0p+0},
        {0x4p-128, 0x4p-1024, 0x0p+0},
        {0x4p-128, 0x8p-972, 0x0p+0},
        {0x4p-1024, 0x4p-128, 0x0p+0},
        {0x4p-1024, 0x4p-1024, 0x0p+0},
        {0x4p-1024, 0x8p-972, 0x0p+0},
        {0x8p-972, 0x4p-128, 0x0p+0},
        {0x8p-972, 0x4p-1024, 0x0p+0},
        {0x8p-972, 0x8p-972, 0x0p+0},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-128},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x4p-1024},
        {0xf.fffffp+124, 0xf.fffffp+124, 0x8p-972},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.fffffp+124, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.fffffp+124, 0x8p-972},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-128},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x4p-1024},
        {0xf.ffffffffffff8p+1020, 0xf.ffffffffffff8p+1020, 0x8p-972},
        {0x2.fffp+12, 0x1.000002p+0, 0x1.ffffp-24},
        {0x1.fffp+0, 0x1.00001p+0, 0x1.fffp+0},
        {0xc.d5e6fp+124, 0x2.6af378p-128, 0x1.f08948p+0},
        {0x1.9abcdep+100, 0x2.6af378p-128, 0x3.e1129p-28},
        {0xf.fffffp+124, 0x1.001p+0, 0xf.fffffp+124},
        {0xf.fffffp+124, 0x1.fffffep+0, 0xf.fffffp+124},
        {0xf.fffffp+124, 0x2p+0, 0xf.fffffp+124},
        {0x5p-128, 0x8.00002p-4, 0x1p-128},
        {0x7.ffffep-128, 0x8.00001p-4, 0x8p-152},
        {0x8p-152, 0x8p-4, 0x3.fffff8p-128},
        {0x8p-152, 0x8.8p-4, 0x3.fffff8p-128},
        {0x8p-152, 0x8p-152, 0x8p+124},
        {0x8p-152, 0x8p-152, 0x4p-128},
        {0x8p-152, 0x8p-152, 0x3.fffff8p-128},
        {0x8p-152, 0x8p-152, 0x8p-152},
        {0xf.ffp-4, 0xf.ffp-4, 0xf.fep-4},
        {0x4.000008p-128, 0x4.000008p-28, 0x8p+124},
        {0x4.000008p-128, 0x4.000008p-28, 0x8p+100},
        {0x2.fep+12, 0x1.0000000000001p+0, 0x1.ffep-48},
        {0x1.fffp+0, 0x1.0000000000001p+0, 0x1.fffp+0},
        {0x1.0000002p+0, 0xf.fffffep-4, 0x1p-300},
        {0xe.f56df7797f768p+1020, 0x3.7ab6fbbcbfbb4p-1024,
         0x3.40bf1803497f6p+0},
        {0x1.deadbeef2feedp+900, 0x3.7ab6fbbcbfbb4p-1024,
         0x6.817e300692fecp-124},
        {0xf.ffffffffffff8p+1020, 0x1.001p+0, 0xf.ffffffffffff8p+1020},
        {0xf.ffffffffffff8p+1020, 0x1.fffffffffffffp+0,
         0xf.ffffffffffff8p+1020},
        {0xf.ffffffffffff8p+1020, 0x2p+0, 0xf.ffffffffffff8p+1020},
        {0x5.a827999fcef3p-540, 0x5.a827999fcef3p-540, 0x0p+0},
        {0x3.bd5b7dde5fddap-496, 0x3.bd5b7dde5fddap-496, 0xd.fc352bc352bap-992},
        {0x3.bd5b7dde5fddap-504, 0x3.bd5b7dde5fddap-504,
         0xd.fc352bc352bap-1008},
        {0x8p-540, 0x4p-540, 0x4p-1076},
        {0x1.7fffff8p-968, 0x4p-108, 0x4p-1048},
        {0x2.8000008p-968, 0x4p-108, 0x4p-1048},
        {0x2.8p-968, 0x4p-108, 0x4p-1048},
        {0x2.33956cdae7c2ep-960, 0x3.8e211518bfea2p-108,
         0x2.02c2b59766d9p-1024},
        {0x3.a5d5dadd1d3a6p-980, 0x2.9c0cd8c5593bap-64, 0x2.49179ac00d15p-1024},
        {0x2.2a7aca1773e0cp-908, 0x9.6809186a42038p-128, 0x2.c9e356b3f0fp-1024},
        {0x3.ffffffffffffep-712, 0x3.ffffffffffffep-276,
         0x3.fffffc0000ffep-984},
        {0x5p-1024, 0x8.000000000001p-4, 0x1p-1024},
        {0x7.ffffffffffffp-1024, 0x8.0000000000008p-4, 0x4p-1076},
        {0x4p-1076, 0x8p-4, 0x3.ffffffffffffcp-1024},
        {0x4p-1076, 0x8.8p-4, 0x3.ffffffffffffcp-1024},
        {0x4p-1076, 0x4p-1076, 0x8p+1020},
        {0x4p-1076, 0x4p-1076, 0x4p-1024},
        {0x4p-1076, 0x4p-1076, 0x3.ffffffffffffcp-1024},
        {0x4p-1076, 0x4p-1076, 0x4p-1076},
        {0xf.ffffffffffff8p-4, 0xf.ffffffffffff8p-4, 0xf.ffffffffffffp-4},
        {0x4.0000000000004p-1024, 0x2.0000000000002p-56, 0x8p+1020},
        {0x4.0000000000004p-1024, 0x2.0000000000002p-56, 0x4p+968},
        {0x7.fffff8p-128, 0x3.fffffcp+24, 0xf.fffffp+124},
        {0x7.ffffffffffffcp-1024, 0x7.ffffffffffffcp+52,
         0xf.ffffffffffff8p+1020},
    };

    for (int i = 0; i < N; ++i) {
      for (int signs = 0; signs < 7; ++signs) {
        double a = (signs & 4) ? -INPUTS[i].a : INPUTS[i].a;
        double b = (signs & 2) ? -INPUTS[i].b : INPUTS[i].b;
        double c = (signs & 1) ? -INPUTS[i].c : INPUTS[i].c;
        mpfr::TernaryInput<double> input{a, b, c};
        ASSERT_MPFR_MATCH_ALL_ROUNDING(mpfr::Operation::Fma, input,
                                       LIBC_NAMESPACE::fma(a, b, c), 0.5);
      }
    }
  }
};

TEST_F(LlvmLibcFmaTest, SubnormalRange) {
  test_subnormal_range(&LIBC_NAMESPACE::fma);
}

TEST_F(LlvmLibcFmaTest, NormalRange) {
  test_normal_range(&LIBC_NAMESPACE::fma);
}

TEST_F(LlvmLibcFmaTest, ExtraValues) { test_more_values(); }