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