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