1 /* $NetBSD: t_fpsetmask.c,v 1.16 2016/03/12 11:55:14 martin Exp $ */
2
3 /*-
4 * Copyright (c) 1995 The NetBSD Foundation, Inc.
5 * All rights reserved.
6 *
7 * Redistribution and use in source and binary forms, with or without
8 * modification, are permitted provided that the following conditions
9 * are met:
10 * 1. Redistributions of source code must retain the above copyright
11 * notice, this list of conditions and the following disclaimer.
12 * 2. Redistributions in binary form must reproduce the above copyright
13 * notice, this list of conditions and the following disclaimer in the
14 * documentation and/or other materials provided with the distribution.
15 *
16 * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
17 * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
18 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
19 * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
20 * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
21 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
22 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
23 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
24 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
25 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
26 * POSSIBILITY OF SUCH DAMAGE.
27 */
28
29 #include <sys/param.h>
30
31 #include <atf-c.h>
32
33 #include <stdio.h>
34 #include <signal.h>
35 #include <float.h>
36 #include <setjmp.h>
37 #include <stdlib.h>
38 #include <string.h>
39
40 #include "isqemu.h"
41
42 #ifndef _FLOAT_IEEE754
43
44 ATF_TC(no_test);
ATF_TC_HEAD(no_test,tc)45 ATF_TC_HEAD(no_test, tc)
46 {
47
48 atf_tc_set_md_var(tc, "descr", "Dummy test case");
49 }
50
ATF_TC_BODY(no_test,tc)51 ATF_TC_BODY(no_test, tc)
52 {
53
54 atf_tc_skip("Test not available on this architecture.");
55 }
56
57 #else /* defined(_FLOAT_IEEE754) */
58
59 #include <ieeefp.h>
60
61 #if __arm__ && !__SOFTFP__
62 /*
63 * Some NEON fpus do not implement IEEE exception handling,
64 * skip these tests if running on them and compiled for
65 * hard float.
66 */
67 #define FPU_PREREQ() \
68 if (0 == fpsetmask(fpsetmask(FP_X_INV))) \
69 atf_tc_skip("FPU does not implement exception handling");
70 #endif
71
72 #ifndef FPU_PREREQ
73 #define FPU_PREREQ() /* nothing */
74 #endif
75
76 void sigfpe(int, siginfo_t *, void *);
77
78 volatile sig_atomic_t signal_caught;
79 volatile int sicode;
80
81 static volatile const float f_one = 1.0;
82 static volatile const float f_zero = 0.0;
83 static volatile const double d_one = 1.0;
84 static volatile const double d_zero = 0.0;
85 static volatile const long double ld_one = 1.0;
86 static volatile const long double ld_zero = 0.0;
87
88 static volatile const float f_huge = FLT_MAX;
89 static volatile const float f_tiny = FLT_MIN;
90 static volatile const double d_huge = DBL_MAX;
91 static volatile const double d_tiny = DBL_MIN;
92 static volatile const long double ld_huge = LDBL_MAX;
93 static volatile const long double ld_tiny = LDBL_MIN;
94
95 static volatile float f_x;
96 static volatile double d_x;
97 static volatile long double ld_x;
98
99 /* trip divide by zero */
100 static void
f_dz(void)101 f_dz(void)
102 {
103
104 f_x = f_one / f_zero;
105 }
106
107 static void
d_dz(void)108 d_dz(void)
109 {
110
111 d_x = d_one / d_zero;
112 }
113
114 static void
ld_dz(void)115 ld_dz(void)
116 {
117
118 ld_x = ld_one / ld_zero;
119 }
120
121 /* trip invalid operation */
122 static void
d_inv(void)123 d_inv(void)
124 {
125
126 d_x = d_zero / d_zero;
127 }
128
129 static void
ld_inv(void)130 ld_inv(void)
131 {
132
133 ld_x = ld_zero / ld_zero;
134 }
135
136 static void
f_inv(void)137 f_inv(void)
138 {
139
140 f_x = f_zero / f_zero;
141 }
142
143 /* trip overflow */
144 static void
f_ofl(void)145 f_ofl(void)
146 {
147
148 f_x = f_huge * f_huge;
149 }
150
151 static void
d_ofl(void)152 d_ofl(void)
153 {
154
155 d_x = d_huge * d_huge;
156 }
157
158 static void
ld_ofl(void)159 ld_ofl(void)
160 {
161
162 ld_x = ld_huge * ld_huge;
163 }
164
165 /* trip underflow */
166 static void
f_ufl(void)167 f_ufl(void)
168 {
169
170 f_x = f_tiny * f_tiny;
171 }
172
173 static void
d_ufl(void)174 d_ufl(void)
175 {
176
177 d_x = d_tiny * d_tiny;
178 }
179
180 static void
ld_ufl(void)181 ld_ufl(void)
182 {
183
184 ld_x = ld_tiny * ld_tiny;
185 }
186
187 struct ops {
188 void (*op)(void);
189 fp_except mask;
190 int sicode;
191 };
192
193 static const struct ops float_ops[] = {
194 { f_dz, FP_X_DZ, FPE_FLTDIV },
195 { f_inv, FP_X_INV, FPE_FLTINV },
196 { f_ofl, FP_X_OFL, FPE_FLTOVF },
197 { f_ufl, FP_X_UFL, FPE_FLTUND },
198 { NULL, 0, 0 }
199 };
200
201 static const struct ops double_ops[] = {
202 { d_dz, FP_X_DZ, FPE_FLTDIV },
203 { d_inv, FP_X_INV, FPE_FLTINV },
204 { d_ofl, FP_X_OFL, FPE_FLTOVF },
205 { d_ufl, FP_X_UFL, FPE_FLTUND },
206 { NULL, 0, 0 }
207 };
208
209 static const struct ops long_double_ops[] = {
210 { ld_dz, FP_X_DZ, FPE_FLTDIV },
211 { ld_inv, FP_X_INV, FPE_FLTINV },
212 { ld_ofl, FP_X_OFL, FPE_FLTOVF },
213 { ld_ufl, FP_X_UFL, FPE_FLTUND },
214 { NULL, 0, 0 }
215 };
216
217 static sigjmp_buf b;
218
219 static void
fpsetmask_masked(const struct ops * test_ops)220 fpsetmask_masked(const struct ops *test_ops)
221 {
222 struct sigaction sa;
223 fp_except ex1, ex2;
224 const struct ops *t;
225
226 /* mask all exceptions, clear history */
227 fpsetmask(0);
228 fpsetsticky(0);
229
230 /* set up signal handler */
231 sa.sa_sigaction = sigfpe;
232 sigemptyset(&sa.sa_mask);
233 sa.sa_flags = SA_SIGINFO;
234 sigaction(SIGFPE, &sa, 0);
235 signal_caught = 0;
236
237 /*
238 * exceptions masked, check whether "sticky" bits are set correctly
239 */
240 for (t = test_ops; t->op != NULL; t++) {
241 (*t->op)();
242 ex1 = fpgetsticky();
243 ATF_CHECK_EQ(ex1 & t->mask, t->mask);
244 ATF_CHECK_EQ(signal_caught, 0);
245
246 /* check correct fpsetsticky() behaviour */
247 ex2 = fpsetsticky(0);
248 ATF_CHECK_EQ(fpgetsticky(), 0);
249 ATF_CHECK_EQ(ex1, ex2);
250 }
251 }
252
253 /* force delayed exceptions to be delivered */
254 #define BARRIER() fpsetmask(0); f_x = f_one * f_one
255
256 static void
fpsetmask_unmasked(const struct ops * test_ops)257 fpsetmask_unmasked(const struct ops *test_ops)
258 {
259 struct sigaction sa;
260 int r;
261 const struct ops *volatile t;
262
263 /* mask all exceptions, clear history */
264 fpsetmask(0);
265 fpsetsticky(0);
266
267 /* set up signal handler */
268 sa.sa_sigaction = sigfpe;
269 sigemptyset(&sa.sa_mask);
270 sa.sa_flags = SA_SIGINFO;
271 sigaction(SIGFPE, &sa, 0);
272 signal_caught = 0;
273
274 /*
275 * exception unmasked, check SIGFPE delivery and correct siginfo
276 */
277 for (t = test_ops; t->op != NULL; t++) {
278 fpsetmask(t->mask);
279 r = sigsetjmp(b, 1);
280 if (!r) {
281 (*t->op)();
282 BARRIER();
283 }
284 ATF_CHECK_EQ(signal_caught, 1);
285 ATF_CHECK_EQ(sicode, t->sicode);
286 signal_caught = 0;
287 }
288 }
289
290 void
sigfpe(int s,siginfo_t * si,void * c)291 sigfpe(int s, siginfo_t *si, void *c)
292 {
293 signal_caught = 1;
294 sicode = si->si_code;
295 siglongjmp(b, 1);
296 }
297
298 #define TEST(m, t) \
299 ATF_TC(m##_##t); \
300 \
301 ATF_TC_HEAD(m##_##t, tc) \
302 { \
303 \
304 atf_tc_set_md_var(tc, "descr", \
305 "Test " ___STRING(m) " exceptions for " \
306 ___STRING(t) "values"); \
307 } \
308 \
309 ATF_TC_BODY(m##_##t, tc) \
310 { \
311 \
312 FPU_PREREQ(); \
313 \
314 if (strcmp(MACHINE, "macppc") == 0) \
315 atf_tc_expect_fail("PR port-macppc/46319"); \
316 \
317 if (isQEMU()) \
318 atf_tc_expect_fail("PR misc/44767"); \
319 \
320 m(t##_ops); \
321 }
322
323 TEST(fpsetmask_masked, float)
324 TEST(fpsetmask_masked, double)
325 TEST(fpsetmask_masked, long_double)
326 TEST(fpsetmask_unmasked, float)
327 TEST(fpsetmask_unmasked, double)
328 TEST(fpsetmask_unmasked, long_double)
329
330 ATF_TC(fpsetmask_basic);
ATF_TC_HEAD(fpsetmask_basic,tc)331 ATF_TC_HEAD(fpsetmask_basic, tc)
332 {
333 atf_tc_set_md_var(tc, "descr", "A basic test of fpsetmask(3)");
334 }
335
ATF_TC_BODY(fpsetmask_basic,tc)336 ATF_TC_BODY(fpsetmask_basic, tc)
337 {
338 size_t i;
339 fp_except_t msk, lst[] = { FP_X_INV, FP_X_DZ, FP_X_OFL, FP_X_UFL };
340
341 FPU_PREREQ();
342
343 msk = fpgetmask();
344 for (i = 0; i < __arraycount(lst); i++) {
345 fpsetmask(msk | lst[i]);
346 ATF_CHECK((fpgetmask() & lst[i]) != 0);
347 fpsetmask(msk & ~lst[i]);
348 ATF_CHECK((fpgetmask() & lst[i]) == 0);
349 }
350
351 }
352
353 #endif /* defined(_FLOAT_IEEE754) */
354
ATF_TP_ADD_TCS(tp)355 ATF_TP_ADD_TCS(tp)
356 {
357
358 #ifndef _FLOAT_IEEE754
359 ATF_TP_ADD_TC(tp, no_test);
360 #else
361 ATF_TP_ADD_TC(tp, fpsetmask_basic);
362 ATF_TP_ADD_TC(tp, fpsetmask_masked_float);
363 ATF_TP_ADD_TC(tp, fpsetmask_masked_double);
364 ATF_TP_ADD_TC(tp, fpsetmask_masked_long_double);
365 ATF_TP_ADD_TC(tp, fpsetmask_unmasked_float);
366 ATF_TP_ADD_TC(tp, fpsetmask_unmasked_double);
367 ATF_TP_ADD_TC(tp, fpsetmask_unmasked_long_double);
368 #endif
369
370 return atf_no_error();
371 }
372