1*309059cfSdjm /* $OpenBSD: bitmap.c,v 1.9 2017/10/20 01:56:39 djm Exp $ */
278e9a853Sdjm /*
378e9a853Sdjm * Copyright (c) 2015 Damien Miller <djm@mindrot.org>
478e9a853Sdjm *
578e9a853Sdjm * Permission to use, copy, modify, and distribute this software for any
678e9a853Sdjm * purpose with or without fee is hereby granted, provided that the above
778e9a853Sdjm * copyright notice and this permission notice appear in all copies.
878e9a853Sdjm *
978e9a853Sdjm * THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL WARRANTIES
1078e9a853Sdjm * WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED WARRANTIES OF
1178e9a853Sdjm * MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR
1278e9a853Sdjm * ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL DAMAGES OR ANY DAMAGES
1378e9a853Sdjm * WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN
1478e9a853Sdjm * ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF
1578e9a853Sdjm * OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
1678e9a853Sdjm */
1778e9a853Sdjm
1878e9a853Sdjm #include <sys/types.h>
1978e9a853Sdjm #include <string.h>
2078e9a853Sdjm #include <stdlib.h>
2178e9a853Sdjm
2278e9a853Sdjm #include "bitmap.h"
2378e9a853Sdjm
2478e9a853Sdjm #define BITMAP_WTYPE u_int
2578e9a853Sdjm #define BITMAP_MAX (1<<24)
2678e9a853Sdjm #define BITMAP_BYTES (sizeof(BITMAP_WTYPE))
2778e9a853Sdjm #define BITMAP_BITS (sizeof(BITMAP_WTYPE) * 8)
2878e9a853Sdjm #define BITMAP_WMASK ((BITMAP_WTYPE)BITMAP_BITS - 1)
2978e9a853Sdjm struct bitmap {
3078e9a853Sdjm BITMAP_WTYPE *d;
3178e9a853Sdjm size_t len; /* number of words allocated */
3278e9a853Sdjm size_t top; /* index of top word allocated */
3378e9a853Sdjm };
3478e9a853Sdjm
3578e9a853Sdjm struct bitmap *
bitmap_new(void)3678e9a853Sdjm bitmap_new(void)
3778e9a853Sdjm {
3878e9a853Sdjm struct bitmap *ret;
3978e9a853Sdjm
4078e9a853Sdjm if ((ret = calloc(1, sizeof(*ret))) == NULL)
4178e9a853Sdjm return NULL;
4278e9a853Sdjm if ((ret->d = calloc(1, BITMAP_BYTES)) == NULL) {
4378e9a853Sdjm free(ret);
4478e9a853Sdjm return NULL;
4578e9a853Sdjm }
4678e9a853Sdjm ret->len = 1;
4778e9a853Sdjm ret->top = 0;
4878e9a853Sdjm return ret;
4978e9a853Sdjm }
5078e9a853Sdjm
5178e9a853Sdjm void
bitmap_free(struct bitmap * b)5278e9a853Sdjm bitmap_free(struct bitmap *b)
5378e9a853Sdjm {
5478e9a853Sdjm if (b != NULL && b->d != NULL) {
5557dd3233Sdjm bitmap_zero(b);
5678e9a853Sdjm free(b->d);
5757dd3233Sdjm b->d = NULL;
5878e9a853Sdjm }
5978e9a853Sdjm free(b);
6078e9a853Sdjm }
6178e9a853Sdjm
6278e9a853Sdjm void
bitmap_zero(struct bitmap * b)6378e9a853Sdjm bitmap_zero(struct bitmap *b)
6478e9a853Sdjm {
6578e9a853Sdjm memset(b->d, 0, b->len * BITMAP_BYTES);
6678e9a853Sdjm b->top = 0;
6778e9a853Sdjm }
6878e9a853Sdjm
6978e9a853Sdjm int
bitmap_test_bit(struct bitmap * b,u_int n)7078e9a853Sdjm bitmap_test_bit(struct bitmap *b, u_int n)
7178e9a853Sdjm {
7278e9a853Sdjm if (b->top >= b->len)
7378e9a853Sdjm return 0; /* invalid */
7478e9a853Sdjm if (b->len == 0 || (n / BITMAP_BITS) > b->top)
7578e9a853Sdjm return 0;
7678e9a853Sdjm return (b->d[n / BITMAP_BITS] >> (n & BITMAP_WMASK)) & 1;
7778e9a853Sdjm }
7878e9a853Sdjm
7978e9a853Sdjm static int
reserve(struct bitmap * b,u_int n)8078e9a853Sdjm reserve(struct bitmap *b, u_int n)
8178e9a853Sdjm {
8278e9a853Sdjm BITMAP_WTYPE *tmp;
8378e9a853Sdjm size_t nlen;
8478e9a853Sdjm
8578e9a853Sdjm if (b->top >= b->len || n > BITMAP_MAX)
8678e9a853Sdjm return -1; /* invalid */
8778e9a853Sdjm nlen = (n / BITMAP_BITS) + 1;
8878e9a853Sdjm if (b->len < nlen) {
89bea1fb3cSdjm if ((tmp = recallocarray(b->d, b->len,
90bea1fb3cSdjm nlen, BITMAP_BYTES)) == NULL)
9178e9a853Sdjm return -1;
9278e9a853Sdjm b->d = tmp;
9378e9a853Sdjm b->len = nlen;
9478e9a853Sdjm }
9578e9a853Sdjm return 0;
9678e9a853Sdjm }
9778e9a853Sdjm
9878e9a853Sdjm int
bitmap_set_bit(struct bitmap * b,u_int n)9978e9a853Sdjm bitmap_set_bit(struct bitmap *b, u_int n)
10078e9a853Sdjm {
10178e9a853Sdjm int r;
10278e9a853Sdjm size_t offset;
10378e9a853Sdjm
10478e9a853Sdjm if ((r = reserve(b, n)) != 0)
10578e9a853Sdjm return r;
10678e9a853Sdjm offset = n / BITMAP_BITS;
10778e9a853Sdjm if (offset > b->top)
10878e9a853Sdjm b->top = offset;
10978e9a853Sdjm b->d[offset] |= (BITMAP_WTYPE)1 << (n & BITMAP_WMASK);
11078e9a853Sdjm return 0;
11178e9a853Sdjm }
11278e9a853Sdjm
11378e9a853Sdjm /* Resets b->top to point to the most significant bit set in b->d */
11478e9a853Sdjm static void
retop(struct bitmap * b)11578e9a853Sdjm retop(struct bitmap *b)
11678e9a853Sdjm {
11778e9a853Sdjm if (b->top >= b->len)
11878e9a853Sdjm return;
11978e9a853Sdjm while (b->top > 0 && b->d[b->top] == 0)
12078e9a853Sdjm b->top--;
12178e9a853Sdjm }
12278e9a853Sdjm
12378e9a853Sdjm void
bitmap_clear_bit(struct bitmap * b,u_int n)12478e9a853Sdjm bitmap_clear_bit(struct bitmap *b, u_int n)
12578e9a853Sdjm {
12678e9a853Sdjm size_t offset;
12778e9a853Sdjm
12878e9a853Sdjm if (b->top >= b->len || n > BITMAP_MAX)
12978e9a853Sdjm return; /* invalid */
13078e9a853Sdjm offset = n / BITMAP_BITS;
13178e9a853Sdjm if (offset > b->top)
13278e9a853Sdjm return;
13378e9a853Sdjm b->d[offset] &= ~((BITMAP_WTYPE)1 << (n & BITMAP_WMASK));
13478e9a853Sdjm /* The top may have changed as a result of the clear */
13578e9a853Sdjm retop(b);
13678e9a853Sdjm }
13778e9a853Sdjm
13878e9a853Sdjm size_t
bitmap_nbits(struct bitmap * b)13978e9a853Sdjm bitmap_nbits(struct bitmap *b)
14078e9a853Sdjm {
14178e9a853Sdjm size_t bits;
14278e9a853Sdjm BITMAP_WTYPE w;
14378e9a853Sdjm
14478e9a853Sdjm retop(b);
14578e9a853Sdjm if (b->top >= b->len)
14678e9a853Sdjm return 0; /* invalid */
14778e9a853Sdjm if (b->len == 0 || (b->top == 0 && b->d[0] == 0))
14878e9a853Sdjm return 0;
14978e9a853Sdjm /* Find MSB set */
15078e9a853Sdjm w = b->d[b->top];
15178e9a853Sdjm bits = (b->top + 1) * BITMAP_BITS;
15278e9a853Sdjm while (!(w & ((BITMAP_WTYPE)1 << (BITMAP_BITS - 1)))) {
15378e9a853Sdjm w <<= 1;
15478e9a853Sdjm bits--;
15578e9a853Sdjm }
15678e9a853Sdjm return bits;
15778e9a853Sdjm }
15878e9a853Sdjm
15978e9a853Sdjm size_t
bitmap_nbytes(struct bitmap * b)16078e9a853Sdjm bitmap_nbytes(struct bitmap *b)
16178e9a853Sdjm {
16278e9a853Sdjm return (bitmap_nbits(b) + 7) / 8;
16378e9a853Sdjm }
16478e9a853Sdjm
16578e9a853Sdjm int
bitmap_to_string(struct bitmap * b,void * p,size_t l)16678e9a853Sdjm bitmap_to_string(struct bitmap *b, void *p, size_t l)
16778e9a853Sdjm {
16878e9a853Sdjm u_char *s = (u_char *)p;
16978e9a853Sdjm size_t i, j, k, need = bitmap_nbytes(b);
17078e9a853Sdjm
17178e9a853Sdjm if (l < need || b->top >= b->len)
17278e9a853Sdjm return -1;
17378e9a853Sdjm if (l > need)
17478e9a853Sdjm l = need;
17578e9a853Sdjm /* Put the bytes from LSB backwards */
17678e9a853Sdjm for (i = k = 0; i < b->top + 1; i++) {
17778e9a853Sdjm for (j = 0; j < BITMAP_BYTES; j++) {
17878e9a853Sdjm if (k >= l)
17978e9a853Sdjm break;
18078e9a853Sdjm s[need - 1 - k++] = (b->d[i] >> (j * 8)) & 0xff;
18178e9a853Sdjm }
18278e9a853Sdjm }
18378e9a853Sdjm return 0;
18478e9a853Sdjm }
18578e9a853Sdjm
18678e9a853Sdjm int
bitmap_from_string(struct bitmap * b,const void * p,size_t l)18778e9a853Sdjm bitmap_from_string(struct bitmap *b, const void *p, size_t l)
18878e9a853Sdjm {
18978e9a853Sdjm int r;
19078e9a853Sdjm size_t i, offset, shift;
1917154f74eSdjm const u_char *s = (const u_char *)p;
19278e9a853Sdjm
19378e9a853Sdjm if (l > BITMAP_MAX / 8)
19478e9a853Sdjm return -1;
19578e9a853Sdjm if ((r = reserve(b, l * 8)) != 0)
19678e9a853Sdjm return r;
19778e9a853Sdjm bitmap_zero(b);
19878e9a853Sdjm if (l == 0)
19978e9a853Sdjm return 0;
20078e9a853Sdjm b->top = offset = ((l + (BITMAP_BYTES - 1)) / BITMAP_BYTES) - 1;
20178e9a853Sdjm shift = ((l + (BITMAP_BYTES - 1)) % BITMAP_BYTES) * 8;
20278e9a853Sdjm for (i = 0; i < l; i++) {
20378e9a853Sdjm b->d[offset] |= (BITMAP_WTYPE)s[i] << shift;
20478e9a853Sdjm if (shift == 0) {
20578e9a853Sdjm offset--;
20678e9a853Sdjm shift = BITMAP_BITS - 8;
20778e9a853Sdjm } else
20878e9a853Sdjm shift -= 8;
20978e9a853Sdjm }
21078e9a853Sdjm retop(b);
21178e9a853Sdjm return 0;
21278e9a853Sdjm }
213