xref: /freebsd-src/contrib/netbsd-tests/include/sys/t_bitops.c (revision 57718be8fa0bd5edc11ab9a72e68cc71982939a6)
1*57718be8SEnji Cooper /*	$NetBSD: t_bitops.c,v 1.16 2012/12/07 02:28:19 christos Exp $ */
2*57718be8SEnji Cooper 
3*57718be8SEnji Cooper /*-
4*57718be8SEnji Cooper  * Copyright (c) 2011, 2012 The NetBSD Foundation, Inc.
5*57718be8SEnji Cooper  * All rights reserved.
6*57718be8SEnji Cooper  *
7*57718be8SEnji Cooper  * This code is derived from software contributed to The NetBSD Foundation
8*57718be8SEnji Cooper  * by Christos Zoulas and Jukka Ruohonen.
9*57718be8SEnji Cooper  *
10*57718be8SEnji Cooper  * Redistribution and use in source and binary forms, with or without
11*57718be8SEnji Cooper  * modification, are permitted provided that the following conditions
12*57718be8SEnji Cooper  * are met:
13*57718be8SEnji Cooper  * 1. Redistributions of source code must retain the above copyright
14*57718be8SEnji Cooper  *    notice, this list of conditions and the following disclaimer.
15*57718be8SEnji Cooper  * 2. Redistributions in binary form must reproduce the above copyright
16*57718be8SEnji Cooper  *    notice, this list of conditions and the following disclaimer in the
17*57718be8SEnji Cooper  *    documentation and/or other materials provided with the distribution.
18*57718be8SEnji Cooper  *
19*57718be8SEnji Cooper  * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
20*57718be8SEnji Cooper  * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
21*57718be8SEnji Cooper  * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
22*57718be8SEnji Cooper  * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
23*57718be8SEnji Cooper  * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
24*57718be8SEnji Cooper  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
25*57718be8SEnji Cooper  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
26*57718be8SEnji Cooper  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
27*57718be8SEnji Cooper  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
28*57718be8SEnji Cooper  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
29*57718be8SEnji Cooper  * POSSIBILITY OF SUCH DAMAGE.
30*57718be8SEnji Cooper  */
31*57718be8SEnji Cooper #include <sys/cdefs.h>
32*57718be8SEnji Cooper __RCSID("$NetBSD: t_bitops.c,v 1.16 2012/12/07 02:28:19 christos Exp $");
33*57718be8SEnji Cooper 
34*57718be8SEnji Cooper #include <atf-c.h>
35*57718be8SEnji Cooper 
36*57718be8SEnji Cooper #include <sys/cdefs.h>
37*57718be8SEnji Cooper #include <sys/bitops.h>
38*57718be8SEnji Cooper 
39*57718be8SEnji Cooper #include <math.h>
40*57718be8SEnji Cooper #include <inttypes.h>
41*57718be8SEnji Cooper #include <stdlib.h>
42*57718be8SEnji Cooper #include <string.h>
43*57718be8SEnji Cooper 
44*57718be8SEnji Cooper static const struct {
45*57718be8SEnji Cooper 	uint32_t	val;
46*57718be8SEnji Cooper 	int		ffs;
47*57718be8SEnji Cooper 	int		fls;
48*57718be8SEnji Cooper } bits[] = {
49*57718be8SEnji Cooper 
50*57718be8SEnji Cooper 	{ 0x00, 0, 0 }, { 0x01, 1, 1 },	{ 0x02, 2, 2 },	{ 0x03, 1, 2 },
51*57718be8SEnji Cooper 	{ 0x04, 3, 3 }, { 0x05, 1, 3 },	{ 0x06, 2, 3 },	{ 0x07, 1, 3 },
52*57718be8SEnji Cooper 	{ 0x08, 4, 4 }, { 0x09, 1, 4 },	{ 0x0A, 2, 4 },	{ 0x0B, 1, 4 },
53*57718be8SEnji Cooper 	{ 0x0C, 3, 4 }, { 0x0D, 1, 4 },	{ 0x0E, 2, 4 },	{ 0x0F, 1, 4 },
54*57718be8SEnji Cooper 
55*57718be8SEnji Cooper 	{ 0x10, 5, 5 },	{ 0x11, 1, 5 },	{ 0x12, 2, 5 },	{ 0x13, 1, 5 },
56*57718be8SEnji Cooper 	{ 0x14, 3, 5 },	{ 0x15, 1, 5 },	{ 0x16, 2, 5 },	{ 0x17, 1, 5 },
57*57718be8SEnji Cooper 	{ 0x18, 4, 5 },	{ 0x19, 1, 5 },	{ 0x1A, 2, 5 },	{ 0x1B, 1, 5 },
58*57718be8SEnji Cooper 	{ 0x1C, 3, 5 },	{ 0x1D, 1, 5 },	{ 0x1E, 2, 5 },	{ 0x1F, 1, 5 },
59*57718be8SEnji Cooper 
60*57718be8SEnji Cooper 	{ 0xF0, 5, 8 },	{ 0xF1, 1, 8 },	{ 0xF2, 2, 8 },	{ 0xF3, 1, 8 },
61*57718be8SEnji Cooper 	{ 0xF4, 3, 8 },	{ 0xF5, 1, 8 },	{ 0xF6, 2, 8 },	{ 0xF7, 1, 8 },
62*57718be8SEnji Cooper 	{ 0xF8, 4, 8 },	{ 0xF9, 1, 8 },	{ 0xFA, 2, 8 },	{ 0xFB, 1, 8 },
63*57718be8SEnji Cooper 	{ 0xFC, 3, 8 },	{ 0xFD, 1, 8 },	{ 0xFE, 2, 8 },	{ 0xFF, 1, 8 },
64*57718be8SEnji Cooper 
65*57718be8SEnji Cooper };
66*57718be8SEnji Cooper 
67*57718be8SEnji Cooper ATF_TC(bitmap_basic);
68*57718be8SEnji Cooper ATF_TC_HEAD(bitmap_basic, tc)
69*57718be8SEnji Cooper {
70*57718be8SEnji Cooper         atf_tc_set_md_var(tc, "descr", "A basic test of __BITMAP_*");
71*57718be8SEnji Cooper }
72*57718be8SEnji Cooper 
73*57718be8SEnji Cooper ATF_TC_BODY(bitmap_basic, tc)
74*57718be8SEnji Cooper {
75*57718be8SEnji Cooper 	__BITMAP_TYPE(, uint32_t, 65536) bm;
76*57718be8SEnji Cooper 	__BITMAP_ZERO(&bm);
77*57718be8SEnji Cooper 
78*57718be8SEnji Cooper 	ATF_REQUIRE(__BITMAP_SIZE(uint32_t, 65536) == 2048);
79*57718be8SEnji Cooper 
80*57718be8SEnji Cooper 	ATF_REQUIRE(__BITMAP_SHIFT(uint32_t) == 5);
81*57718be8SEnji Cooper 
82*57718be8SEnji Cooper 	ATF_REQUIRE(__BITMAP_MASK(uint32_t) == 31);
83*57718be8SEnji Cooper 
84*57718be8SEnji Cooper 	for (size_t i = 0; i < 65536; i += 2)
85*57718be8SEnji Cooper 		__BITMAP_SET(i, &bm);
86*57718be8SEnji Cooper 
87*57718be8SEnji Cooper 	for (size_t i = 0; i < 2048; i++)
88*57718be8SEnji Cooper 		ATF_REQUIRE(bm._b[i] == 0x55555555);
89*57718be8SEnji Cooper 
90*57718be8SEnji Cooper 	for (size_t i = 0; i < 65536; i++)
91*57718be8SEnji Cooper 		if (i & 1)
92*57718be8SEnji Cooper 			ATF_REQUIRE(!__BITMAP_ISSET(i, &bm));
93*57718be8SEnji Cooper 		else {
94*57718be8SEnji Cooper 			ATF_REQUIRE(__BITMAP_ISSET(i, &bm));
95*57718be8SEnji Cooper 			__BITMAP_CLR(i, &bm);
96*57718be8SEnji Cooper 		}
97*57718be8SEnji Cooper 
98*57718be8SEnji Cooper 	for (size_t i = 0; i < 65536; i += 2)
99*57718be8SEnji Cooper 		ATF_REQUIRE(!__BITMAP_ISSET(i, &bm));
100*57718be8SEnji Cooper }
101*57718be8SEnji Cooper 
102*57718be8SEnji Cooper ATF_TC(fast_divide32);
103*57718be8SEnji Cooper ATF_TC_HEAD(fast_divide32, tc)
104*57718be8SEnji Cooper {
105*57718be8SEnji Cooper 	atf_tc_set_md_var(tc, "descr", "A basic test of fast_divide32(3)");
106*57718be8SEnji Cooper }
107*57718be8SEnji Cooper 
108*57718be8SEnji Cooper ATF_TC_BODY(fast_divide32, tc)
109*57718be8SEnji Cooper {
110*57718be8SEnji Cooper 	uint32_t a, b, q, r, m;
111*57718be8SEnji Cooper 	uint8_t i, s1, s2;
112*57718be8SEnji Cooper 
113*57718be8SEnji Cooper 	a = 0xFFFF;
114*57718be8SEnji Cooper 	b = 0x000F;
115*57718be8SEnji Cooper 
116*57718be8SEnji Cooper 	fast_divide32_prepare(b, &m, &s1, &s2);
117*57718be8SEnji Cooper 
118*57718be8SEnji Cooper 	q = fast_divide32(a, b, m, s1, s2);
119*57718be8SEnji Cooper 	r = fast_remainder32(a, b, m, s1, s2);
120*57718be8SEnji Cooper 
121*57718be8SEnji Cooper 	ATF_REQUIRE(q == 0x1111 && r == 0);
122*57718be8SEnji Cooper 
123*57718be8SEnji Cooper 	for (i = 1; i < __arraycount(bits); i++) {
124*57718be8SEnji Cooper 
125*57718be8SEnji Cooper 		a = bits[i].val;
126*57718be8SEnji Cooper 		b = bits[i].ffs;
127*57718be8SEnji Cooper 
128*57718be8SEnji Cooper 		fast_divide32_prepare(b, &m, &s1, &s2);
129*57718be8SEnji Cooper 
130*57718be8SEnji Cooper 		q = fast_divide32(a, b, m, s1, s2);
131*57718be8SEnji Cooper 		r = fast_remainder32(a, b, m, s1, s2);
132*57718be8SEnji Cooper 
133*57718be8SEnji Cooper 		ATF_REQUIRE(q == a / b);
134*57718be8SEnji Cooper 		ATF_REQUIRE(r == a % b);
135*57718be8SEnji Cooper 	}
136*57718be8SEnji Cooper }
137*57718be8SEnji Cooper 
138*57718be8SEnji Cooper ATF_TC(ffsfls);
139*57718be8SEnji Cooper ATF_TC_HEAD(ffsfls, tc)
140*57718be8SEnji Cooper {
141*57718be8SEnji Cooper 	atf_tc_set_md_var(tc, "descr", "Test ffs32(3)-family for correctness");
142*57718be8SEnji Cooper }
143*57718be8SEnji Cooper 
144*57718be8SEnji Cooper ATF_TC_BODY(ffsfls, tc)
145*57718be8SEnji Cooper {
146*57718be8SEnji Cooper 	uint8_t i;
147*57718be8SEnji Cooper 
148*57718be8SEnji Cooper 	ATF_REQUIRE(ffs32(0) == 0x00);
149*57718be8SEnji Cooper 	ATF_REQUIRE(fls32(0) == 0x00);
150*57718be8SEnji Cooper 	ATF_REQUIRE(ffs64(0) == 0x00);
151*57718be8SEnji Cooper 	ATF_REQUIRE(fls64(0) == 0x00);
152*57718be8SEnji Cooper 
153*57718be8SEnji Cooper 	ATF_REQUIRE(ffs32(UINT32_MAX) == 0x01);
154*57718be8SEnji Cooper 	ATF_REQUIRE(fls32(UINT32_MAX) == 0x20);
155*57718be8SEnji Cooper 	ATF_REQUIRE(ffs64(UINT64_MAX) == 0x01);
156*57718be8SEnji Cooper 	ATF_REQUIRE(fls64(UINT64_MAX) == 0x40);
157*57718be8SEnji Cooper 
158*57718be8SEnji Cooper 	for (i = 1; i < __arraycount(bits); i++) {
159*57718be8SEnji Cooper 
160*57718be8SEnji Cooper 		ATF_REQUIRE(ffs32(bits[i].val) == bits[i].ffs);
161*57718be8SEnji Cooper 		ATF_REQUIRE(fls32(bits[i].val) == bits[i].fls);
162*57718be8SEnji Cooper 		ATF_REQUIRE(ffs64(bits[i].val) == bits[i].ffs);
163*57718be8SEnji Cooper 		ATF_REQUIRE(fls64(bits[i].val) == bits[i].fls);
164*57718be8SEnji Cooper 
165*57718be8SEnji Cooper 		ATF_REQUIRE(ffs32(bits[i].val << 1) == bits[i].ffs + 1);
166*57718be8SEnji Cooper 		ATF_REQUIRE(fls32(bits[i].val << 1) == bits[i].fls + 1);
167*57718be8SEnji Cooper 		ATF_REQUIRE(ffs64(bits[i].val << 1) == bits[i].ffs + 1);
168*57718be8SEnji Cooper 		ATF_REQUIRE(fls64(bits[i].val << 1) == bits[i].fls + 1);
169*57718be8SEnji Cooper 
170*57718be8SEnji Cooper 		ATF_REQUIRE(ffs32(bits[i].val << 9) == bits[i].ffs + 9);
171*57718be8SEnji Cooper 		ATF_REQUIRE(fls32(bits[i].val << 9) == bits[i].fls + 9);
172*57718be8SEnji Cooper 		ATF_REQUIRE(ffs64(bits[i].val << 9) == bits[i].ffs + 9);
173*57718be8SEnji Cooper 		ATF_REQUIRE(fls64(bits[i].val << 9) == bits[i].fls + 9);
174*57718be8SEnji Cooper 	}
175*57718be8SEnji Cooper }
176*57718be8SEnji Cooper 
177*57718be8SEnji Cooper ATF_TC(ilog2_basic);
178*57718be8SEnji Cooper ATF_TC_HEAD(ilog2_basic, tc)
179*57718be8SEnji Cooper {
180*57718be8SEnji Cooper 	atf_tc_set_md_var(tc, "descr", "Test ilog2(3) for correctness");
181*57718be8SEnji Cooper }
182*57718be8SEnji Cooper 
183*57718be8SEnji Cooper ATF_TC_BODY(ilog2_basic, tc)
184*57718be8SEnji Cooper {
185*57718be8SEnji Cooper 	uint64_t i, x;
186*57718be8SEnji Cooper 
187*57718be8SEnji Cooper 	for (i = x = 0; i < 64; i++) {
188*57718be8SEnji Cooper 
189*57718be8SEnji Cooper 		x = (uint64_t)1 << i;
190*57718be8SEnji Cooper 
191*57718be8SEnji Cooper 		ATF_REQUIRE(i == (uint64_t)ilog2(x));
192*57718be8SEnji Cooper 	}
193*57718be8SEnji Cooper }
194*57718be8SEnji Cooper 
195*57718be8SEnji Cooper ATF_TC(ilog2_log2);
196*57718be8SEnji Cooper ATF_TC_HEAD(ilog2_log2, tc)
197*57718be8SEnji Cooper {
198*57718be8SEnji Cooper 	atf_tc_set_md_var(tc, "descr", "Test log2(3) vs. ilog2(3)");
199*57718be8SEnji Cooper }
200*57718be8SEnji Cooper 
201*57718be8SEnji Cooper ATF_TC_BODY(ilog2_log2, tc)
202*57718be8SEnji Cooper {
203*57718be8SEnji Cooper #ifdef __vax__
204*57718be8SEnji Cooper 	atf_tc_skip("Test is unavailable on vax because of lack of log2()");
205*57718be8SEnji Cooper #else
206*57718be8SEnji Cooper 	double  x, y;
207*57718be8SEnji Cooper 	uint64_t i;
208*57718be8SEnji Cooper 
209*57718be8SEnji Cooper 	/*
210*57718be8SEnji Cooper 	 * This may fail under QEMU; see PR misc/44767.
211*57718be8SEnji Cooper 	 */
212*57718be8SEnji Cooper 	for (i = 1; i < UINT32_MAX; i += UINT16_MAX) {
213*57718be8SEnji Cooper 
214*57718be8SEnji Cooper 		x = log2(i);
215*57718be8SEnji Cooper 		y = (double)(ilog2(i));
216*57718be8SEnji Cooper 
217*57718be8SEnji Cooper 		ATF_REQUIRE(ceil(x) >= y);
218*57718be8SEnji Cooper 
219*57718be8SEnji Cooper 		if (fabs(floor(x) - y) > 1.0e-40) {
220*57718be8SEnji Cooper 			atf_tc_expect_fail("PR misc/44767");
221*57718be8SEnji Cooper 			atf_tc_fail("log2(%"PRIu64") != "
222*57718be8SEnji Cooper 			    "ilog2(%"PRIu64")", i, i);
223*57718be8SEnji Cooper 		}
224*57718be8SEnji Cooper 	}
225*57718be8SEnji Cooper #endif
226*57718be8SEnji Cooper }
227*57718be8SEnji Cooper 
228*57718be8SEnji Cooper ATF_TP_ADD_TCS(tp)
229*57718be8SEnji Cooper {
230*57718be8SEnji Cooper 
231*57718be8SEnji Cooper         ATF_TP_ADD_TC(tp, bitmap_basic);
232*57718be8SEnji Cooper 	ATF_TP_ADD_TC(tp, fast_divide32);
233*57718be8SEnji Cooper 	ATF_TP_ADD_TC(tp, ffsfls);
234*57718be8SEnji Cooper 	ATF_TP_ADD_TC(tp, ilog2_basic);
235*57718be8SEnji Cooper 	ATF_TP_ADD_TC(tp, ilog2_log2);
236*57718be8SEnji Cooper 
237*57718be8SEnji Cooper 	return atf_no_error();
238*57718be8SEnji Cooper }
239