bitmap.h revision 1.2.34.1 1 1.2.34.1 pgoyette /* $NetBSD: bitmap.h,v 1.2.34.1 2018/09/06 06:56:36 pgoyette Exp $ */
2 1.2 riastrad
3 1.2 riastrad /*-
4 1.2.34.1 pgoyette * Copyright (c) 2018 The NetBSD Foundation, Inc.
5 1.2 riastrad * All rights reserved.
6 1.2 riastrad *
7 1.2 riastrad * This code is derived from software contributed to The NetBSD Foundation
8 1.2 riastrad * by Taylor R. Campbell.
9 1.2 riastrad *
10 1.2 riastrad * Redistribution and use in source and binary forms, with or without
11 1.2 riastrad * modification, are permitted provided that the following conditions
12 1.2 riastrad * are met:
13 1.2 riastrad * 1. Redistributions of source code must retain the above copyright
14 1.2 riastrad * notice, this list of conditions and the following disclaimer.
15 1.2 riastrad * 2. Redistributions in binary form must reproduce the above copyright
16 1.2 riastrad * notice, this list of conditions and the following disclaimer in the
17 1.2 riastrad * documentation and/or other materials provided with the distribution.
18 1.2 riastrad *
19 1.2 riastrad * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
20 1.2 riastrad * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
21 1.2 riastrad * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
22 1.2 riastrad * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
23 1.2 riastrad * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
24 1.2 riastrad * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
25 1.2 riastrad * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
26 1.2 riastrad * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
27 1.2 riastrad * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
28 1.2 riastrad * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
29 1.2 riastrad * POSSIBILITY OF SUCH DAMAGE.
30 1.2 riastrad */
31 1.2 riastrad
32 1.2 riastrad #ifndef _LINUX_BITMAP_H_
33 1.2 riastrad #define _LINUX_BITMAP_H_
34 1.2 riastrad
35 1.2.34.1 pgoyette #include <sys/param.h>
36 1.2.34.1 pgoyette #include <sys/types.h>
37 1.2.34.1 pgoyette #include <sys/systm.h>
38 1.2.34.1 pgoyette
39 1.2.34.1 pgoyette /*
40 1.2.34.1 pgoyette * bitmap_zero(bitmap, nbits)
41 1.2.34.1 pgoyette *
42 1.2.34.1 pgoyette * Zero a bitmap that was allocated to have nbits bits. Yes, this
43 1.2.34.1 pgoyette * zeros bits past nbits.
44 1.2.34.1 pgoyette */
45 1.2.34.1 pgoyette static inline void
46 1.2.34.1 pgoyette bitmap_zero(unsigned long *bitmap, size_t nbits)
47 1.2.34.1 pgoyette {
48 1.2.34.1 pgoyette const size_t bpl = NBBY * sizeof(*bitmap);
49 1.2.34.1 pgoyette size_t n = howmany(nbits, bpl);
50 1.2.34.1 pgoyette
51 1.2.34.1 pgoyette memset(bitmap, 0, n * sizeof(*bitmap));
52 1.2.34.1 pgoyette }
53 1.2.34.1 pgoyette
54 1.2.34.1 pgoyette /*
55 1.2.34.1 pgoyette * bitmap_empty(bitmap, nbits)
56 1.2.34.1 pgoyette *
57 1.2.34.1 pgoyette * Return true if all bits at 0, 1, 2, ..., nbits-2, nbits-1 are
58 1.2.34.1 pgoyette * 0, or false if any of them is 1.
59 1.2.34.1 pgoyette */
60 1.2.34.1 pgoyette static inline bool
61 1.2.34.1 pgoyette bitmap_empty(const unsigned long *bitmap, size_t nbits)
62 1.2.34.1 pgoyette {
63 1.2.34.1 pgoyette const size_t bpl = NBBY * sizeof(*bitmap);
64 1.2.34.1 pgoyette
65 1.2.34.1 pgoyette for (; nbits >= bpl; nbits -= bpl) {
66 1.2.34.1 pgoyette if (*bitmap++)
67 1.2.34.1 pgoyette return false;
68 1.2.34.1 pgoyette }
69 1.2.34.1 pgoyette
70 1.2.34.1 pgoyette if (nbits) {
71 1.2.34.1 pgoyette if (*bitmap & ~(~0UL << nbits))
72 1.2.34.1 pgoyette return false;
73 1.2.34.1 pgoyette }
74 1.2.34.1 pgoyette
75 1.2.34.1 pgoyette return true;
76 1.2.34.1 pgoyette }
77 1.2.34.1 pgoyette
78 1.2.34.1 pgoyette /*
79 1.2.34.1 pgoyette * bitmap_weight(bitmap, nbits)
80 1.2.34.1 pgoyette *
81 1.2.34.1 pgoyette * Compute the number of 1 bits at 0, 1, 2, ..., nbits-2, nbits-1.
82 1.2.34.1 pgoyette */
83 1.2.34.1 pgoyette static inline int
84 1.2.34.1 pgoyette bitmap_weight(const unsigned long *bitmap, size_t nbits)
85 1.2.34.1 pgoyette {
86 1.2.34.1 pgoyette const size_t bpl = NBBY * sizeof(*bitmap);
87 1.2.34.1 pgoyette int weight = 0;
88 1.2.34.1 pgoyette
89 1.2.34.1 pgoyette for (; nbits >= bpl; nbits -= bpl)
90 1.2.34.1 pgoyette weight += popcountl(*bitmap++);
91 1.2.34.1 pgoyette if (nbits)
92 1.2.34.1 pgoyette weight += popcountl(*bitmap & ~(~0UL << nbits));
93 1.2.34.1 pgoyette
94 1.2.34.1 pgoyette return weight;
95 1.2.34.1 pgoyette }
96 1.2.34.1 pgoyette
97 1.2.34.1 pgoyette /*
98 1.2.34.1 pgoyette * bitmap_set(bitmap, startbit, nbits)
99 1.2.34.1 pgoyette *
100 1.2.34.1 pgoyette * Set bits at startbit, startbit+1, ..., startbit+nbits-2,
101 1.2.34.1 pgoyette * startbit+nbits-1 to 1.
102 1.2.34.1 pgoyette */
103 1.2.34.1 pgoyette static inline void
104 1.2.34.1 pgoyette bitmap_set(unsigned long *bitmap, size_t startbit, size_t nbits)
105 1.2.34.1 pgoyette {
106 1.2.34.1 pgoyette const size_t bpl = NBBY * sizeof(*bitmap);
107 1.2.34.1 pgoyette unsigned long *p = bitmap + startbit/bpl;
108 1.2.34.1 pgoyette unsigned initial = startbit%bpl;
109 1.2.34.1 pgoyette
110 1.2.34.1 pgoyette /* Handle an initial odd word if any. */
111 1.2.34.1 pgoyette if (initial) {
112 1.2.34.1 pgoyette /* Does the whole thing fit in a single word? */
113 1.2.34.1 pgoyette if (nbits <= bpl - initial) {
114 1.2.34.1 pgoyette /* Yes: just set nbits starting at initial. */
115 1.2.34.1 pgoyette *p |= ~(~0ULL << nbits) << initial;
116 1.2.34.1 pgoyette return;
117 1.2.34.1 pgoyette }
118 1.2.34.1 pgoyette /* Nope: set all bits above initial, and advance. */
119 1.2.34.1 pgoyette *p++ |= ~0ULL << initial;
120 1.2.34.1 pgoyette nbits -= bpl - initial;
121 1.2.34.1 pgoyette }
122 1.2.34.1 pgoyette
123 1.2.34.1 pgoyette /* Set the middle part to all bits 1. */
124 1.2.34.1 pgoyette for (; nbits >= bpl; nbits -= bpl)
125 1.2.34.1 pgoyette *p++ = ~0UL;
126 1.2.34.1 pgoyette
127 1.2.34.1 pgoyette /* Handle a final odd word if any by setting its low nbits. */
128 1.2.34.1 pgoyette if (nbits)
129 1.2.34.1 pgoyette *p |= ~(~0ULL << nbits);
130 1.2.34.1 pgoyette }
131 1.2.34.1 pgoyette
132 1.2.34.1 pgoyette /*
133 1.2.34.1 pgoyette * bitmap_clear(bitmap, startbit, nbits)
134 1.2.34.1 pgoyette *
135 1.2.34.1 pgoyette * Clear bits at startbit, startbit+1, ..., startbit+nbits-2,
136 1.2.34.1 pgoyette * startbit+nbits-1, replacing them by 0.
137 1.2.34.1 pgoyette */
138 1.2.34.1 pgoyette static inline void
139 1.2.34.1 pgoyette bitmap_clear(unsigned long *bitmap, size_t startbit, size_t nbits)
140 1.2.34.1 pgoyette {
141 1.2.34.1 pgoyette const size_t bpl = NBBY * sizeof(*bitmap);
142 1.2.34.1 pgoyette unsigned long *p = bitmap + startbit/bpl;
143 1.2.34.1 pgoyette unsigned initial = startbit%bpl;
144 1.2.34.1 pgoyette
145 1.2.34.1 pgoyette /* Handle an initial odd word if any. */
146 1.2.34.1 pgoyette if (initial) {
147 1.2.34.1 pgoyette /* Does the whole thing fit in a single word? */
148 1.2.34.1 pgoyette if (nbits <= bpl - initial) {
149 1.2.34.1 pgoyette /* Yes: just clear nbits starting at initial. */
150 1.2.34.1 pgoyette *p &= ~(~(~0ULL << nbits) << initial);
151 1.2.34.1 pgoyette return;
152 1.2.34.1 pgoyette }
153 1.2.34.1 pgoyette /* Nope: clear all bits above initial, and advance. */
154 1.2.34.1 pgoyette *p++ &= ~(~0ULL << initial);
155 1.2.34.1 pgoyette nbits -= bpl - initial;
156 1.2.34.1 pgoyette }
157 1.2.34.1 pgoyette
158 1.2.34.1 pgoyette /* Zero the middle part. */
159 1.2.34.1 pgoyette for (; nbits >= bpl; nbits -= bpl)
160 1.2.34.1 pgoyette *p++ = 0UL;
161 1.2.34.1 pgoyette
162 1.2.34.1 pgoyette /* Handle a final odd word if any by clearing its low nbits. */
163 1.2.34.1 pgoyette if (nbits)
164 1.2.34.1 pgoyette *p &= ~0ULL << nbits;
165 1.2.34.1 pgoyette }
166 1.2.34.1 pgoyette
167 1.2.34.1 pgoyette /*
168 1.2.34.1 pgoyette * bitmap_and(dst, src1, src2, nbits)
169 1.2.34.1 pgoyette *
170 1.2.34.1 pgoyette * Set dst to be the bitwise AND of src1 and src2, all bitmaps
171 1.2.34.1 pgoyette * allocated to have nbits bits. Yes, this modifies bits past
172 1.2.34.1 pgoyette * nbits. Any pair of {dst, src1, src2} may be aliases.
173 1.2.34.1 pgoyette */
174 1.2.34.1 pgoyette static inline void
175 1.2.34.1 pgoyette bitmap_and(unsigned long *dst, const unsigned long *src1,
176 1.2.34.1 pgoyette const unsigned long *src2, size_t nbits)
177 1.2.34.1 pgoyette {
178 1.2.34.1 pgoyette const size_t bpl = NBBY * sizeof(unsigned long);
179 1.2.34.1 pgoyette size_t n = howmany(nbits, bpl);
180 1.2.34.1 pgoyette
181 1.2.34.1 pgoyette while (n --> 0)
182 1.2.34.1 pgoyette *dst++ = *src1++ & *src2++;
183 1.2.34.1 pgoyette }
184 1.2.34.1 pgoyette
185 1.2.34.1 pgoyette /*
186 1.2.34.1 pgoyette * bitmap_or(dst, src1, src2, nbits)
187 1.2.34.1 pgoyette *
188 1.2.34.1 pgoyette * Set dst to be the bitwise inclusive-OR of src1 and src2, all
189 1.2.34.1 pgoyette * bitmaps allocated to have nbits bits. Yes, this modifies bits
190 1.2.34.1 pgoyette * past nbits. Any pair of {dst, src1, src2} may be aliases.
191 1.2.34.1 pgoyette */
192 1.2.34.1 pgoyette static inline void
193 1.2.34.1 pgoyette bitmap_or(unsigned long *dst, const unsigned long *src1,
194 1.2.34.1 pgoyette const unsigned long *src2, size_t nbits)
195 1.2.34.1 pgoyette {
196 1.2.34.1 pgoyette const size_t bpl = NBBY * sizeof(unsigned long);
197 1.2.34.1 pgoyette size_t n = howmany(nbits, bpl);
198 1.2.34.1 pgoyette
199 1.2.34.1 pgoyette while (n --> 0)
200 1.2.34.1 pgoyette *dst++ = *src1++ | *src2++;
201 1.2.34.1 pgoyette }
202 1.2.34.1 pgoyette
203 1.2 riastrad #endif /* _LINUX_BITMAP_H_ */
204