1 /*
2  * Written by Maya Rashish <maya@NetBSD.org>
3  * Public domain.
4  *
5  * Testing IEEE-754 rounding modes (and lrint)
6  */
7 
8 #include <atf-c.h>
9 #include <fenv.h>
10 #ifdef __HAVE_FENV
11 #include <math.h>
12 #include <stdio.h>
13 #include <stdlib.h>
14 
15 /*#pragma STDC FENV_ACCESS ON gcc?? */
16 
17 #define INT 9223L
18 
19 #define EPSILON 0.001
20 
21 static const struct {
22 	int round_mode;
23 	double input;
24 	long int expected;
25 } values[] = {
26 	{ FE_DOWNWARD,		3.7,		3},
27 	{ FE_DOWNWARD,		-3.7,		-4},
28 	{ FE_DOWNWARD,		+0,		0},
29 	{ FE_DOWNWARD,		-INT-0.01,	-INT-1},
30 	{ FE_DOWNWARD,		+INT-0.01,	INT-1},
31 	{ FE_DOWNWARD,		-INT+0.01,	-INT},
32 	{ FE_DOWNWARD,		+INT+0.01,	INT},
33 #if 0 /* cpu bugs? */
34 	{ FE_DOWNWARD,		-0,		-1},
35 
36 	{ FE_UPWARD,		+0,		1},
37 #endif
38 	{ FE_UPWARD,		-0,		0},
39 	{ FE_UPWARD,		-123.7,		-123},
40 	{ FE_UPWARD,		123.999,	124},
41 	{ FE_UPWARD,		-INT-0.01,	-INT},
42 	{ FE_UPWARD,		+INT-0.01,	INT},
43 	{ FE_UPWARD,		-INT+0.01,	-INT+1},
44 	{ FE_UPWARD,		+INT+0.01,	INT+1},
45 
46 	{ FE_TOWARDZERO,	1.99,		1},
47 	{ FE_TOWARDZERO,	-1.99,		-1},
48 	{ FE_TOWARDZERO,	0.2,		0},
49 	{ FE_TOWARDZERO,	INT+0.01,	INT},
50 	{ FE_TOWARDZERO,	INT-0.01,	INT - 1},
51 	{ FE_TOWARDZERO,	-INT+0.01,	-INT + 1},
52 	{ FE_TOWARDZERO,	+0,		0},
53 	{ FE_TOWARDZERO,	-0,		0},
54 
55 	{ FE_TONEAREST,		-INT-0.01,	-INT},
56 	{ FE_TONEAREST,		+INT-0.01,	INT},
57 	{ FE_TONEAREST,		-INT+0.01,	-INT},
58 	{ FE_TONEAREST,		+INT+0.01,	INT},
59 	{ FE_TONEAREST,		-INT-0.501,	-INT-1},
60 	{ FE_TONEAREST,		+INT-0.501,	INT-1},
61 	{ FE_TONEAREST,		-INT+0.501,	-INT+1},
62 	{ FE_TONEAREST,		+INT+0.501,	INT+1},
63 	{ FE_TONEAREST,		+0,		0},
64 	{ FE_TONEAREST,		-0,		0},
65 };
66 
67 ATF_TC(fe_round);
68 ATF_TC_HEAD(fe_round, tc)
69 {
70 	atf_tc_set_md_var(tc, "descr","Checking IEEE 754 rounding modes using lrint");
71 }
72 
73 ATF_TC_BODY(fe_round, tc)
74 {
75 	long int received;
76 
77 	for (unsigned int i = 0; i < __arraycount(values); i++) {
78 		fesetround(values[i].round_mode);
79 
80 		received = lrint(values[i].input);
81 		ATF_CHECK_MSG(
82 		    (labs(received - values[i].expected) < EPSILON),
83 		    "lrint rounding wrong, difference too large\n"
84 		    "input: %f (index %d): got %ld, expected %ld\n",
85 		    values[i].input, i, received, values[i].expected);
86 
87 		/* Do we get the same rounding mode out? */
88 		ATF_CHECK_MSG(
89 		    (fegetround() == values[i].round_mode),
90 		    "Didn't get the same rounding mode out!\n"
91 		    "(index %d) fed in %d rounding mode, got %d out\n",
92 		    i, fegetround(), values[i].round_mode);
93 	}
94 }
95 
96 ATF_TP_ADD_TCS(tp)
97 {
98 
99 	ATF_TP_ADD_TC(tp, fe_round);
100 
101 	return atf_no_error();
102 }
103 #else
104 ATF_TC(t_nofe_round);
105 
106 ATF_TC_HEAD(t_nofe_round, tc)
107 {
108 	atf_tc_set_md_var(tc, "descr",
109 	    "dummy test case - no fenv.h support");
110 }
111 
112 
113 ATF_TC_BODY(t_nofe_round, tc)
114 {
115 	atf_tc_skip("no fenv.h support on this architecture");
116 }
117 
118 ATF_TP_ADD_TCS(tp)
119 {
120 	ATF_TP_ADD_TC(tp, t_nofe_round);
121 	return atf_no_error();
122 }
123 
124 #endif
125