1274079Sngie/*	$NetBSD: t_fpsetmask.c,v 1.14 2014/11/04 00:20:19 justin Exp $ */
2272343Sngie
3272343Sngie/*-
4272343Sngie * Copyright (c) 1995 The NetBSD Foundation, Inc.
5272343Sngie * All rights reserved.
6272343Sngie *
7272343Sngie * Redistribution and use in source and binary forms, with or without
8272343Sngie * modification, are permitted provided that the following conditions
9272343Sngie * are met:
10272343Sngie * 1. Redistributions of source code must retain the above copyright
11272343Sngie *    notice, this list of conditions and the following disclaimer.
12272343Sngie * 2. Redistributions in binary form must reproduce the above copyright
13272343Sngie *    notice, this list of conditions and the following disclaimer in the
14272343Sngie *    documentation and/or other materials provided with the distribution.
15272343Sngie *
16272343Sngie * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
17272343Sngie * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
18272343Sngie * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
19272343Sngie * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
20272343Sngie * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
21272343Sngie * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
22272343Sngie * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
23272343Sngie * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
24272343Sngie * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
25272343Sngie * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
26272343Sngie * POSSIBILITY OF SUCH DAMAGE.
27272343Sngie */
28272343Sngie
29272343Sngie#include <sys/param.h>
30272343Sngie
31272343Sngie#include <atf-c.h>
32272343Sngie
33272343Sngie#include <stdio.h>
34272343Sngie#include <signal.h>
35272343Sngie#include <float.h>
36272343Sngie#include <setjmp.h>
37272343Sngie#include <stdlib.h>
38272343Sngie#include <string.h>
39272343Sngie
40272343Sngie#include "isqemu.h"
41272343Sngie
42272343Sngie#ifndef _FLOAT_IEEE754
43272343Sngie
44272343SngieATF_TC(no_test);
45272343SngieATF_TC_HEAD(no_test, tc)
46272343Sngie{
47272343Sngie
48272343Sngie	atf_tc_set_md_var(tc, "descr", "Dummy test case");
49272343Sngie}
50272343Sngie
51272343SngieATF_TC_BODY(no_test, tc)
52272343Sngie{
53272343Sngie
54272343Sngie	atf_tc_skip("Test not available on this architecture.");
55272343Sngie}
56272343Sngie
57272343Sngie#else /* defined(_FLOAT_IEEE754) */
58272343Sngie
59272343Sngie#include <ieeefp.h>
60272343Sngie
61272343Sngieconst char *skip_mesg;
62272343Sngieconst char *skip_arch;
63272343Sngie
64272343Sngievoid		sigfpe(int, siginfo_t *, void *);
65272343Sngie
66272343Sngievolatile sig_atomic_t signal_caught;
67272343Sngievolatile int sicode;
68272343Sngie
69272343Sngiestatic volatile const float	f_one   = 1.0;
70272343Sngiestatic volatile const float	f_zero  = 0.0;
71272343Sngiestatic volatile const double	d_one   = 1.0;
72272343Sngiestatic volatile const double	d_zero  = 0.0;
73272343Sngiestatic volatile const long double ld_one  = 1.0;
74272343Sngiestatic volatile const long double ld_zero = 0.0;
75272343Sngie
76272343Sngiestatic volatile const float	f_huge = FLT_MAX;
77272343Sngiestatic volatile const float	f_tiny = FLT_MIN;
78272343Sngiestatic volatile const double	d_huge = DBL_MAX;
79272343Sngiestatic volatile const double	d_tiny = DBL_MIN;
80272343Sngiestatic volatile const long double ld_huge = LDBL_MAX;
81272343Sngiestatic volatile const long double ld_tiny = LDBL_MIN;
82272343Sngie
83272343Sngiestatic volatile float f_x;
84272343Sngiestatic volatile double d_x;
85272343Sngiestatic volatile long double ld_x;
86272343Sngie
87272343Sngie/* trip divide by zero */
88272343Sngiestatic void
89272343Sngief_dz(void)
90272343Sngie{
91272343Sngie
92272343Sngie	f_x = f_one / f_zero;
93272343Sngie}
94272343Sngie
95272343Sngiestatic void
96272343Sngied_dz(void)
97272343Sngie{
98272343Sngie
99272343Sngie	d_x = d_one / d_zero;
100272343Sngie}
101272343Sngie
102272343Sngiestatic void
103272343Sngield_dz(void)
104272343Sngie{
105272343Sngie
106272343Sngie	ld_x = ld_one / ld_zero;
107272343Sngie}
108272343Sngie
109272343Sngie/* trip invalid operation */
110272343Sngiestatic void
111272343Sngied_inv(void)
112272343Sngie{
113272343Sngie
114272343Sngie	d_x = d_zero / d_zero;
115272343Sngie}
116272343Sngie
117272343Sngiestatic void
118272343Sngield_inv(void)
119272343Sngie{
120272343Sngie
121272343Sngie	ld_x = ld_zero / ld_zero;
122272343Sngie}
123272343Sngie
124272343Sngiestatic void
125272343Sngief_inv(void)
126272343Sngie{
127272343Sngie
128272343Sngie	f_x = f_zero / f_zero;
129272343Sngie}
130272343Sngie
131272343Sngie/* trip overflow */
132272343Sngiestatic void
133272343Sngief_ofl(void)
134272343Sngie{
135272343Sngie
136272343Sngie	f_x = f_huge * f_huge;
137272343Sngie}
138272343Sngie
139272343Sngiestatic void
140272343Sngied_ofl(void)
141272343Sngie{
142272343Sngie
143272343Sngie	d_x = d_huge * d_huge;
144272343Sngie}
145272343Sngie
146272343Sngiestatic void
147272343Sngield_ofl(void)
148272343Sngie{
149272343Sngie
150272343Sngie	ld_x = ld_huge * ld_huge;
151272343Sngie}
152272343Sngie
153272343Sngie/* trip underflow */
154272343Sngiestatic void
155272343Sngief_ufl(void)
156272343Sngie{
157272343Sngie
158272343Sngie	f_x = f_tiny * f_tiny;
159272343Sngie}
160272343Sngie
161272343Sngiestatic void
162272343Sngied_ufl(void)
163272343Sngie{
164272343Sngie
165272343Sngie	d_x = d_tiny * d_tiny;
166272343Sngie}
167272343Sngie
168272343Sngiestatic void
169272343Sngield_ufl(void)
170272343Sngie{
171272343Sngie
172272343Sngie	ld_x = ld_tiny * ld_tiny;
173272343Sngie}
174272343Sngie
175272343Sngiestruct ops {
176272343Sngie	void (*op)(void);
177272343Sngie	fp_except mask;
178272343Sngie	int sicode;
179272343Sngie};
180272343Sngie
181272343Sngiestatic const struct ops float_ops[] = {
182272343Sngie	{ f_dz, FP_X_DZ, FPE_FLTDIV },
183272343Sngie	{ f_inv, FP_X_INV, FPE_FLTINV },
184272343Sngie	{ f_ofl, FP_X_OFL, FPE_FLTOVF },
185272343Sngie	{ f_ufl, FP_X_UFL, FPE_FLTUND },
186272343Sngie	{ NULL, 0, 0 }
187272343Sngie};
188272343Sngie
189272343Sngiestatic const struct ops double_ops[] = {
190272343Sngie	{ d_dz, FP_X_DZ, FPE_FLTDIV },
191272343Sngie	{ d_inv, FP_X_INV, FPE_FLTINV },
192272343Sngie	{ d_ofl, FP_X_OFL, FPE_FLTOVF },
193272343Sngie	{ d_ufl, FP_X_UFL, FPE_FLTUND },
194272343Sngie	{ NULL, 0, 0 }
195272343Sngie};
196272343Sngie
197272343Sngiestatic const struct ops long_double_ops[] = {
198272343Sngie	{ ld_dz, FP_X_DZ, FPE_FLTDIV },
199272343Sngie	{ ld_inv, FP_X_INV, FPE_FLTINV },
200272343Sngie	{ ld_ofl, FP_X_OFL, FPE_FLTOVF },
201272343Sngie	{ ld_ufl, FP_X_UFL, FPE_FLTUND },
202272343Sngie	{ NULL, 0, 0 }
203272343Sngie};
204272343Sngie
205272343Sngiestatic sigjmp_buf b;
206272343Sngie
207272343Sngiestatic void
208272343Sngiefpsetmask_masked(const struct ops *test_ops)
209272343Sngie{
210272343Sngie	struct sigaction sa;
211272343Sngie	fp_except ex1, ex2;
212272343Sngie	const struct ops *t;
213272343Sngie
214272343Sngie	/* mask all exceptions, clear history */
215272343Sngie	fpsetmask(0);
216272343Sngie	fpsetsticky(0);
217272343Sngie
218272343Sngie	/* set up signal handler */
219272343Sngie	sa.sa_sigaction = sigfpe;
220272343Sngie	sigemptyset(&sa.sa_mask);
221272343Sngie	sa.sa_flags = SA_SIGINFO;
222272343Sngie	sigaction(SIGFPE, &sa, 0);
223272343Sngie	signal_caught = 0;
224272343Sngie
225272343Sngie	/*
226272343Sngie	 * exceptions masked, check whether "sticky" bits are set correctly
227272343Sngie	 */
228272343Sngie	for (t = test_ops; t->op != NULL; t++) {
229272343Sngie		(*t->op)();
230272343Sngie		ex1 = fpgetsticky();
231272343Sngie		ATF_CHECK_EQ(ex1 & t->mask, t->mask);
232272343Sngie		ATF_CHECK_EQ(signal_caught, 0);
233272343Sngie
234272343Sngie		/* check correct fpsetsticky() behaviour */
235272343Sngie		ex2 = fpsetsticky(0);
236272343Sngie		ATF_CHECK_EQ(fpgetsticky(), 0);
237272343Sngie		ATF_CHECK_EQ(ex1, ex2);
238272343Sngie	}
239272343Sngie}
240272343Sngie
241272343Sngie/* force delayed exceptions to be delivered */
242272343Sngie#define BARRIER() fpsetmask(0); f_x = f_one * f_one
243272343Sngie
244272343Sngiestatic void
245272343Sngiefpsetmask_unmasked(const struct ops *test_ops)
246272343Sngie{
247272343Sngie	struct sigaction sa;
248272343Sngie	int r;
249272343Sngie	const struct ops *volatile t;
250272343Sngie
251272343Sngie	/* mask all exceptions, clear history */
252272343Sngie	fpsetmask(0);
253272343Sngie	fpsetsticky(0);
254272343Sngie
255272343Sngie	/* set up signal handler */
256272343Sngie	sa.sa_sigaction = sigfpe;
257272343Sngie	sigemptyset(&sa.sa_mask);
258272343Sngie	sa.sa_flags = SA_SIGINFO;
259272343Sngie	sigaction(SIGFPE, &sa, 0);
260272343Sngie	signal_caught = 0;
261272343Sngie
262272343Sngie	/*
263272343Sngie	 * exception unmasked, check SIGFPE delivery and correct siginfo
264272343Sngie	 */
265272343Sngie	for (t = test_ops; t->op != NULL; t++) {
266272343Sngie		fpsetmask(t->mask);
267272343Sngie		r = sigsetjmp(b, 1);
268272343Sngie		if (!r) {
269272343Sngie			(*t->op)();
270272343Sngie			BARRIER();
271272343Sngie		}
272272343Sngie		ATF_CHECK_EQ(signal_caught, 1);
273272343Sngie		ATF_CHECK_EQ(sicode, t->sicode);
274272343Sngie		signal_caught = 0;
275272343Sngie	}
276272343Sngie}
277272343Sngie
278272343Sngievoid
279272343Sngiesigfpe(int s, siginfo_t *si, void *c)
280272343Sngie{
281272343Sngie	signal_caught = 1;
282272343Sngie	sicode = si->si_code;
283272343Sngie	siglongjmp(b, 1);
284272343Sngie}
285272343Sngie
286272343Sngie#define TEST(m, t)							\
287272343Sngie	ATF_TC(m##_##t);						\
288272343Sngie									\
289272343Sngie	ATF_TC_HEAD(m##_##t, tc)					\
290272343Sngie	{								\
291272343Sngie									\
292272343Sngie		atf_tc_set_md_var(tc, "descr",				\
293272343Sngie		    "Test " ___STRING(m) " exceptions for "		\
294272343Sngie		    ___STRING(t) "values");				\
295272343Sngie	}								\
296272343Sngie									\
297272343Sngie	ATF_TC_BODY(m##_##t, tc)					\
298272343Sngie	{								\
299272343Sngie		if (strcmp(MACHINE, "macppc") == 0)			\
300272343Sngie			atf_tc_expect_fail("PR port-macppc/46319");	\
301272343Sngie									\
302272343Sngie		if (isQEMU())						\
303272343Sngie			atf_tc_expect_fail("PR misc/44767");		\
304272343Sngie									\
305272343Sngie		m(t##_ops);						\
306272343Sngie	}
307272343Sngie
308272343SngieTEST(fpsetmask_masked, float)
309272343SngieTEST(fpsetmask_masked, double)
310272343SngieTEST(fpsetmask_masked, long_double)
311272343SngieTEST(fpsetmask_unmasked, float)
312272343SngieTEST(fpsetmask_unmasked, double)
313272343SngieTEST(fpsetmask_unmasked, long_double)
314272343Sngie
315272343SngieATF_TC(fpsetmask_basic);
316272343SngieATF_TC_HEAD(fpsetmask_basic, tc)
317272343Sngie{
318272343Sngie	atf_tc_set_md_var(tc, "descr", "A basic test of fpsetmask(3)");
319272343Sngie}
320272343Sngie
321272343SngieATF_TC_BODY(fpsetmask_basic, tc)
322272343Sngie{
323272343Sngie	size_t i;
324272343Sngie	fp_except_t msk, lst[] = { FP_X_INV, FP_X_DZ, FP_X_OFL, FP_X_UFL };
325272343Sngie
326272343Sngie	msk = fpgetmask();
327272343Sngie	for (i = 0; i < __arraycount(lst); i++) {
328272343Sngie		fpsetmask(msk | lst[i]);
329272343Sngie		ATF_CHECK((fpgetmask() & lst[i]) != 0);
330272343Sngie		fpsetmask(msk & lst[i]);
331272343Sngie		ATF_CHECK((fpgetmask() & lst[i]) == 0);
332272343Sngie	}
333272343Sngie
334272343Sngie}
335272343Sngie
336272343Sngie#endif /* defined(_FLOAT_IEEE754) */
337272343Sngie
338272343SngieATF_TP_ADD_TCS(tp)
339272343Sngie{
340272343Sngie
341272343Sngie#ifndef _FLOAT_IEEE754
342272343Sngie	ATF_TP_ADD_TC(tp, no_test);
343272343Sngie#else
344272343Sngie	ATF_TP_ADD_TC(tp, fpsetmask_basic);
345272343Sngie	ATF_TP_ADD_TC(tp, fpsetmask_masked_float);
346272343Sngie	ATF_TP_ADD_TC(tp, fpsetmask_masked_double);
347272343Sngie	ATF_TP_ADD_TC(tp, fpsetmask_masked_long_double);
348272343Sngie	ATF_TP_ADD_TC(tp, fpsetmask_unmasked_float);
349272343Sngie	ATF_TP_ADD_TC(tp, fpsetmask_unmasked_double);
350272343Sngie	ATF_TP_ADD_TC(tp, fpsetmask_unmasked_long_double);
351272343Sngie#endif
352272343Sngie
353272343Sngie	return atf_no_error();
354272343Sngie}
355