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
44ATF_TC(no_test);
45ATF_TC_HEAD(no_test, tc)
46{
47
48	atf_tc_set_md_var(tc, "descr", "Dummy test case");
49}
50
51ATF_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
76void		sigfpe(int, siginfo_t *, void *);
77
78volatile sig_atomic_t signal_caught;
79volatile int sicode;
80
81static volatile const float	f_one   = 1.0;
82static volatile const float	f_zero  = 0.0;
83static volatile const double	d_one   = 1.0;
84static volatile const double	d_zero  = 0.0;
85static volatile const long double ld_one  = 1.0;
86static volatile const long double ld_zero = 0.0;
87
88static volatile const float	f_huge = FLT_MAX;
89static volatile const float	f_tiny = FLT_MIN;
90static volatile const double	d_huge = DBL_MAX;
91static volatile const double	d_tiny = DBL_MIN;
92static volatile const long double ld_huge = LDBL_MAX;
93static volatile const long double ld_tiny = LDBL_MIN;
94
95static volatile float f_x;
96static volatile double d_x;
97static volatile long double ld_x;
98
99/* trip divide by zero */
100static void
101f_dz(void)
102{
103
104	f_x = f_one / f_zero;
105}
106
107static void
108d_dz(void)
109{
110
111	d_x = d_one / d_zero;
112}
113
114static void
115ld_dz(void)
116{
117
118	ld_x = ld_one / ld_zero;
119}
120
121/* trip invalid operation */
122static void
123d_inv(void)
124{
125
126	d_x = d_zero / d_zero;
127}
128
129static void
130ld_inv(void)
131{
132
133	ld_x = ld_zero / ld_zero;
134}
135
136static void
137f_inv(void)
138{
139
140	f_x = f_zero / f_zero;
141}
142
143/* trip overflow */
144static void
145f_ofl(void)
146{
147
148	f_x = f_huge * f_huge;
149}
150
151static void
152d_ofl(void)
153{
154
155	d_x = d_huge * d_huge;
156}
157
158static void
159ld_ofl(void)
160{
161
162	ld_x = ld_huge * ld_huge;
163}
164
165/* trip underflow */
166static void
167f_ufl(void)
168{
169
170	f_x = f_tiny * f_tiny;
171}
172
173static void
174d_ufl(void)
175{
176
177	d_x = d_tiny * d_tiny;
178}
179
180static void
181ld_ufl(void)
182{
183
184	ld_x = ld_tiny * ld_tiny;
185}
186
187struct ops {
188	void (*op)(void);
189	fp_except mask;
190	int sicode;
191};
192
193static 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
201static 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
209static 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
217static sigjmp_buf b;
218
219static void
220fpsetmask_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
256static void
257fpsetmask_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
290void
291sigfpe(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
323TEST(fpsetmask_masked, float)
324TEST(fpsetmask_masked, double)
325TEST(fpsetmask_masked, long_double)
326TEST(fpsetmask_unmasked, float)
327TEST(fpsetmask_unmasked, double)
328TEST(fpsetmask_unmasked, long_double)
329
330ATF_TC(fpsetmask_basic);
331ATF_TC_HEAD(fpsetmask_basic, tc)
332{
333	atf_tc_set_md_var(tc, "descr", "A basic test of fpsetmask(3)");
334}
335
336ATF_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
355ATF_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