Support dtor sections on Mach-O.
[libfirm] / ir / adt / bitset.h
index f7efa15..5eb3490 100644 (file)
-/**
- * Bitsets.
+/*
+ * Copyright (C) 1995-2008 University of Karlsruhe.  All right reserved.
+ *
+ * This file is part of libFirm.
+ *
+ * This file may be distributed and/or modified under the terms of the
+ * GNU General Public License version 2 as published by the Free Software
+ * Foundation and appearing in the file LICENSE.GPL included in the
+ * packaging of this file.
+ *
+ * Licensees holding valid libFirm Professional Edition licenses may use
+ * this file in accordance with the libFirm Commercial License.
+ * Agreement provided with the Software.
+ *
+ * This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
+ * WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+ * PURPOSE.
  */
 
-#ifndef __FIRM_BITSET_H
-#define __FIRM_BITSET_H
+/**
+ * @file
+ * @brief   A bitset implementation.
+ * @author  Sebastian Hack
+ * @date    15.10.2004
+ * @version $Id$
+ */
+#ifndef FIRM_ADT_BITSET_H
+#define FIRM_ADT_BITSET_H
 
 #include <stdlib.h>
 #include <stdio.h>
 #include <assert.h>
 #include <string.h>
 
+#include "xmalloc.h"
+#include "bitfiddle.h"
 
-#define INLINE inline
-
-#define BITSET_USE_STD
+typedef unsigned int bitset_pos_t;
 
 #include "bitset_std.h"
 
+#if defined(__GNUC__) && defined(__i386__)
+#include "bitset_ia32.h"
+#endif
+
 typedef struct _bitset_t {
-       unsigned units;
-       bitset_unit_t *data;
+       bitset_pos_t units;
+       bitset_pos_t size;
 } bitset_t;
 
-#define BS_UNIT_SIZE sizeof(bitset_unit_t)
-#define BS_UNIT_SIZE_BITS (BS_UNIT_SIZE * 8)
-#define BS_UNIT_MASK (BS_UNIT_SIZE_BITS - 1)
+#define BS_UNIT_SIZE         sizeof(bitset_unit_t)
+#define BS_UNIT_SIZE_BITS    (BS_UNIT_SIZE * 8)
+#define BS_UNIT_MASK         (BS_UNIT_SIZE_BITS - 1)
 
-
-/**
- * Units needed for a given highest bit.
- * @param highest_bit The highest bit that should be storable.
- * @return The number of units needed.
- */
-#define _bitset_units(highest_bit) (round_up2(highest_bit, BS_UNIT_SIZE_BITS) / BS_UNIT_SIZE_BITS)
-
-/**
- * Compute the size in bytes needed for a bitset.
- * This also include the size for the bitset data structure.
- * @param highest_bit The highest bit that shall be storable.
- * @return The overall amount of bytes needed for that bitset.
- */
-#define _bitset_overall_size(highest_bit) \
-       (sizeof(bitset_t) + _bitset_units(highest_bit) * BS_UNIT_SIZE)
+#define BS_DATA(bs)          ((bitset_unit_t *) ((char *) (bs) + sizeof(bitset_t)))
+#define BS_UNITS(bits)       (round_up2(bits, BS_UNIT_SIZE_BITS) / BS_UNIT_SIZE_BITS)
+#define BS_TOTAL_SIZE(bits)  (sizeof(bitset_t) + BS_UNITS(bits) * BS_UNIT_SIZE)
 
 /**
  * Initialize a bitset.
  * This functions should not be called.
+ *
+ * Note that this function needs three macros which must be provided by the
+ * bitfield implementor:
+ * - _bitset_overall_size(size) The overall size that must be
+ *   allocated for the bitfield in bytes.
+ * - _bitset_units(size) The number of units that will be
+ *   present in the bitfield for a given highest bit.
+ * - _bitset_data_ptr(data, size) This produces as pointer to the
+ *   first unit in the allocated memory area. The main reason for this
+ *   macro is, that some bitset implementors want control over memory
+ *   alignment.
+ *
  * @param area A pointer to memory reserved for the bitset.
- * @param units The number of units that are allocated for the bitset.
+ * @param size The size of the bitset in bits.
  * @return A pointer to the initialized bitset.
  */
-static INLINE bitset_t *_bitset_prepare(void *area, unsigned units)
+static inline bitset_t *_bitset_prepare(void *area, bitset_pos_t size)
 {
        bitset_t *ptr = area;
-       ptr->units = units;
-       ptr->data = (bitset_unit_t *) (ptr + 1);
-       memset(ptr->data, 0, BS_UNIT_SIZE * units);
+       memset(ptr, 0, BS_TOTAL_SIZE(size));
+       ptr->units = BS_UNITS(size);
+       ptr->size  = size;
        return ptr;
 }
 
+/**
+ * Mask out all bits, which are only there, because the number
+ * of bits in the set didn't match a unit size boundary.
+ * @param bs The bitset.
+ * @return The masked bitset.
+ */
+static inline bitset_t *_bitset_mask_highest(bitset_t *bs)
+{
+       bitset_pos_t rest = bs->size & BS_UNIT_MASK;
+       if (rest)
+               BS_DATA(bs)[bs->units - 1] &= (1 << rest) - 1;
+       return bs;
+}
+
+/**
+ * Get the capacity of the bitset in bits.
+ * @param bs The bitset.
+ * @return The capacity in bits of the bitset.
+ */
+#define bitset_capacity(bs) ((bs)->units * BS_UNIT_SIZE_BITS)
+
+/**
+ * Get the size of the bitset in bits.
+ * @note Note the difference between capacity and size.
+ * @param bs The bitset.
+ * @return The highest bit which can be set or cleared plus 1.
+ */
+#define bitset_size(bs)  ((bs)->size)
+
 /**
  * Allocate a bitset on an obstack.
  * @param obst The obstack.
- * @param highest_bit The greatest bit that shall be stored in the set.
+ * @param size The greatest bit that shall be stored in the set.
  * @return A pointer to an empty initialized bitset.
  */
-#define bitset_obstack_alloc(obst,highest_bit) \
-  _bitset_prepare(obstack_alloc(obst, _bitset_overall_size(highest_bit)), _bitset_units(highest_bit))
+#define bitset_obstack_alloc(obst,size) \
+       _bitset_prepare(obstack_alloc(obst, BS_TOTAL_SIZE(size)), size)
 
 /**
  * Allocate a bitset via malloc.
- * @param highest_bit The greatest bit that shall be stored in the set.
+ * @param size The greatest bit that shall be stored in the set.
  * @return A pointer to an empty initialized bitset.
  */
-#define bitset_malloc(highest_bit) \
-       _bitset_prepare(malloc(_bitset_overall_size(highest_bit)), _bitset_units(highest_bit))
+#define bitset_malloc(size) \
+       _bitset_prepare(xmalloc(BS_TOTAL_SIZE(size)), size)
 
 /**
  * Free a bitset allocated with bitset_malloc().
@@ -84,43 +138,12 @@ static INLINE bitset_t *_bitset_prepare(void *area, unsigned units)
 
 /**
  * Allocate a bitset on the stack via alloca.
- * @param highest_bit The greatest bit that shall be stored in the set.
+ * @param size The greatest bit that shall be stored in the set.
  * @return A pointer to an empty initialized bitset.
  */
-#define bitset_alloca(highest_bit) \
-       _bitset_prepare(alloca(_bitset_overall_size(highest_bit)), _bitset_units(highest_bit))
-
-/**
- * Print a bitset to a stream.
- * The bitset is printed as a comma seperated list of bits set.
- * @param file The stream.
- * @param bs The bitset.
- */
-static INLINE void bitset_fprint(FILE *file, bitset_t *bs)
-{
-       const char *prefix = "";
-       int i;
-       unsigned k = 0;
+#define bitset_alloca(size) \
+       _bitset_prepare(alloca(BS_TOTAL_SIZE(size)), size)
 
-       putc('[', file);
-       for(i = 0; i < bs->units; i++) {
-               bitset_unit_t j;
-               bitset_unit_t unit = bs->data[i];
-
-#if 0
-               printf("%s%08x", prefix, unit);
-               prefix=":";
-               continue;
-#endif
-               for(j = 1; j != 0; j <<= 1, k++) {
-                       if(unit & j) {
-                               fprintf(file, "%s%u", prefix, k);
-                               prefix = ",";
-                       }
-               }
-       }
-       putc(']', file);
-}
 
 /**
  * Get the unit which contains a specific bit.
@@ -129,10 +152,10 @@ static INLINE void bitset_fprint(FILE *file, bitset_t *bs)
  * @param bit The bit.
  * @return A pointer to the unit containing the bit.
  */
-static INLINE bitset_unit_t *_bitset_get_unit(const bitset_t *bs, unsigned bit)
+static inline bitset_unit_t *_bitset_get_unit(const bitset_t *bs, bitset_pos_t bit)
 {
-       assert(bit < bs->units * BS_UNIT_SIZE_BITS && "Bit too large");
-       return bs->data + bit / BS_UNIT_SIZE_BITS;
+       assert(bit <= bs->size && "Bit to large");
+       return BS_DATA(bs) + bit / BS_UNIT_SIZE_BITS;
 }
 
 /**
@@ -140,7 +163,7 @@ static INLINE bitset_unit_t *_bitset_get_unit(const bitset_t *bs, unsigned bit)
  * @param bs The bitset.
  * @param bit The bit to set.
  */
-static INLINE void bitset_set(bitset_t *bs, unsigned bit)
+static inline void bitset_set(bitset_t *bs, bitset_pos_t bit)
 {
        bitset_unit_t *unit = _bitset_get_unit(bs, bit);
        _bitset_inside_set(unit, bit & BS_UNIT_MASK);
@@ -151,102 +174,428 @@ static INLINE void bitset_set(bitset_t *bs, unsigned bit)
  * @param bs The bitset.
  * @param bit The bit to clear.
  */
-static INLINE void bitset_clear(bitset_t *bs, unsigned bit)
+static inline void bitset_clear(bitset_t *bs, bitset_pos_t bit)
 {
        bitset_unit_t *unit = _bitset_get_unit(bs, bit);
        _bitset_inside_clear(unit, bit & BS_UNIT_MASK);
 }
 
+/**
+ * Check, if a bit is set.
+ * @param bs The bitset.
+ * @param bit The bit to check for.
+ * @return 1, if the bit was set, 0 if not.
+ */
+static inline int bitset_is_set(const bitset_t *bs, bitset_pos_t bit)
+{
+       bitset_unit_t *unit = _bitset_get_unit(bs, bit);
+       return _bitset_inside_is_set(unit, bit & BS_UNIT_MASK);
+}
+
 /**
  * Flip a bit in a bitset.
  * @param bs The bitset.
  * @param bit The bit to flip.
  */
-static INLINE void bitset_flip(bitset_t *bs, unsigned bit)
+static inline void bitset_flip(bitset_t *bs, bitset_pos_t bit)
 {
        bitset_unit_t *unit = _bitset_get_unit(bs, bit);
        _bitset_inside_flip(unit, bit & BS_UNIT_MASK);
 }
 
 /**
- * Find the smallest bit set in the bitset.
+ * Flip the whole bitset.
  * @param bs The bitset.
- * @return The smallest bit set in the bitset.
  */
-static INLINE unsigned bitset_min(bitset_t *bs)
+static inline void bitset_flip_all(bitset_t *bs)
 {
-       unsigned i, ofs = 0;
-       bitset_unit_t *unit;
-
-       for(i = 0, unit = bs->data; i < bs->units; ++i, ++unit) {
-               unsigned pos = _bitset_inside_ntz(unit);
-               if(pos > 0)
-                       return ofs + pos;
-               ofs += BS_UNIT_SIZE_BITS;
-       }
+       bitset_pos_t i;
+       for(i = 0; i < bs->units; i++)
+               _bitset_inside_flip_unit(&BS_DATA(bs)[i]);
+       _bitset_mask_highest(bs);
+}
 
-       return 0;
+/**
+ * Copy a bitset to another.
+ * @param tgt The target bitset.
+ * @param src The source bitset.
+ * @return The target bitset.
+ */
+static inline bitset_t *bitset_copy(bitset_t *tgt, const bitset_t *src)
+{
+       bitset_pos_t tu = tgt->units;
+       bitset_pos_t su = src->units;
+       bitset_pos_t min_units = tu < su ? tu : su;
+       memcpy(BS_DATA(tgt), BS_DATA(src), min_units * BS_UNIT_SIZE);
+       if(tu > min_units)
+               memset(BS_DATA(tgt) + min_units, 0, BS_UNIT_SIZE * (tu - min_units));
+       return _bitset_mask_highest(tgt);
 }
 
 /**
- * Find the greatest bit set in the bitset.
+ * Find the next set bit from a given bit.
+ * @note Note that if pos is set, pos is returned.
  * @param bs The bitset.
- * @return The greatest bit set in the bitset.
+ * @param pos The bit from which to search for the next set bit.
+ * @return The next set bit from pos on, or -1, if no set bit was found
+ * after pos.
  */
-static INLINE unsigned bitset_max(bitset_t *bs)
+static inline bitset_pos_t _bitset_next(const bitset_t *bs,
+               bitset_pos_t pos, int set)
 {
-       unsigned i, max = 0, ofs = 0;
-       bitset_unit_t *unit;
+       bitset_pos_t unit_number = pos / BS_UNIT_SIZE_BITS;
+       bitset_pos_t res;
+
+       if(pos >= bs->size)
+               return -1;
+
+       {
+               bitset_pos_t bit_in_unit = pos & BS_UNIT_MASK;
+               bitset_pos_t in_unit_mask = (1 << bit_in_unit) - 1;
+
+               /*
+                * Mask out the bits smaller than pos in the current unit.
+                * We are only interested in bits set higher than pos.
+                */
+               bitset_unit_t curr_unit = BS_DATA(bs)[unit_number];
+
+               /*
+                * Find the next bit set in the unit.
+                * Mind that this function returns 0, if the unit is -1 and
+                * counts the bits from 1 on.
+                */
+               bitset_pos_t next_in_this_unit =
+                       _bitset_inside_ntz_value((set ? curr_unit : ~curr_unit) & ~in_unit_mask);
+
+               /* If there is a bit set in the current unit, exit. */
+               if (next_in_this_unit < BS_UNIT_SIZE_BITS) {
+                       res = next_in_this_unit + unit_number * BS_UNIT_SIZE_BITS;
+                       return res < bs->size ? res : (bitset_pos_t) -1;
+               }
 
-       for(i = 0, unit = bs->data; i < bs->units; ++i, ++unit) {
-               unsigned pos = _bitset_inside_nlz(unit);
-               if(pos > 0)
-                       max = ofs + pos;
-               ofs += BS_UNIT_SIZE_BITS;
+               /* Else search for set bits in the next units. */
+               else {
+                       bitset_pos_t i;
+                       for(i = unit_number + 1; i < bs->units; ++i) {
+                               bitset_unit_t data = BS_DATA(bs)[i];
+                               bitset_pos_t first_set =
+                                       _bitset_inside_ntz_value(set ? data : ~data);
+
+                               if (first_set < BS_UNIT_SIZE_BITS) {
+                                       res = first_set + i * BS_UNIT_SIZE_BITS;
+                                       return res < bs->size ? res : (bitset_pos_t) -1;
+                               }
+                       }
+               }
        }
 
-       return max;
+       return -1;
 }
 
+#define bitset_next_clear(bs,pos) _bitset_next((bs), (pos), 0)
+#define bitset_next_set(bs,pos) _bitset_next((bs), (pos), 1)
+
+/**
+ * Convenience macro for bitset iteration.
+ * @param bitset The bitset.
+ * @param elm A unsigned long variable.
+ */
+#define bitset_foreach(bitset,elm) \
+       for(elm = bitset_next_set(bitset,0); elm != (bitset_pos_t) -1; elm = bitset_next_set(bitset,elm+1))
+
+
+#define bitset_foreach_clear(bitset,elm) \
+       for(elm = bitset_next_clear(bitset,0); elm != (bitset_pos_t) -1; elm = bitset_next_clear(bitset,elm+1))
+
 /**
  * Count the bits set.
  * This can also be seen as the cardinality of the set.
  * @param bs The bitset.
  * @return The number of bits set in the bitset.
  */
-static INLINE unsigned bitset_pop(bitset_t *bs)
+static inline unsigned bitset_popcnt(const bitset_t *bs)
 {
-       unsigned i, pop = 0;
+       bitset_pos_t  i;
        bitset_unit_t *unit;
+       unsigned      pop = 0;
 
-       for(i = 0, unit = bs->data; i < bs->units; ++i, ++unit)
+       for (i = 0, unit = BS_DATA(bs); i < bs->units; ++i, ++unit)
                pop += _bitset_inside_pop(unit);
 
        return pop;
 }
 
+/**
+ * Clear the bitset.
+ * This sets all bits to zero.
+ * @param bs The bitset.
+ */
+static inline bitset_t *bitset_clear_all(bitset_t *bs)
+{
+       memset(BS_DATA(bs), 0, BS_UNIT_SIZE * bs->units);
+       return bs;
+}
+
+/**
+ * Set the bitset.
+ * This sets all bits to one.
+ * @param bs The bitset.
+ */
+static inline bitset_t *bitset_set_all(bitset_t *bs)
+{
+       memset(BS_DATA(bs), -1, bs->units * BS_UNIT_SIZE);
+       return _bitset_mask_highest(bs);
+}
+
+/**
+ * Check, if one bitset is contained by another.
+ * That is, each bit set in lhs is also set in rhs.
+ * @param lhs A bitset.
+ * @param rhs Another bitset.
+ * @return 1, if all bits in lhs are also set in rhs, 0 otherwise.
+ */
+static inline int bitset_contains(const bitset_t *lhs, const bitset_t *rhs)
+{
+       bitset_pos_t n = lhs->units < rhs->units ? lhs->units : rhs->units;
+       bitset_pos_t i;
+
+       for(i = 0; i < n; ++i) {
+               bitset_unit_t lu = BS_DATA(lhs)[i];
+               bitset_unit_t ru = BS_DATA(rhs)[i];
+
+               if((lu | ru) & ~ru)
+                       return 0;
+       }
+
+       /*
+        * If the left hand sinde is a larger bitset than rhs,
+        * we have to check, that all extra bits in lhs are 0
+        */
+       if(lhs->units > n) {
+               for(i = n; i < lhs->units; ++i) {
+                       if(BS_DATA(lhs)[i] != 0)
+                               return 0;
+               }
+       }
+
+       return 1;
+}
+
+/**
+ * Treat the bitset as a number and subtract 1.
+ * @param bs The bitset.
+ * @return The same bitset.
+ */
+static inline void bitset_minus1(bitset_t *bs)
+{
+#define _SH (sizeof(bitset_unit_t) * 8 - 1)
+
+       bitset_pos_t i;
+
+       for(i = 0; i < bs->units; ++i) {
+               bitset_unit_t unit = BS_DATA(bs)[i];
+               bitset_unit_t um1  = unit - 1;
+
+               BS_DATA(bs)[i] = um1;
+
+               if(((unit >> _SH) ^ (um1 >> _SH)) == 0)
+                       break;
+       }
+#undef _SH
+}
+
+/**
+ * Check if two bitsets intersect.
+ * @param a The first bitset.
+ * @param b The second bitset.
+ * @return 1 if they have a bit in common, 0 if not.
+ */
+static inline int bitset_intersect(const bitset_t *a, const bitset_t *b)
+{
+       bitset_pos_t n = a->units < b->units ? a->units : b->units;
+       bitset_pos_t i;
+
+       for (i = 0; i < n; ++i)
+               if (BS_DATA(a)[i] & BS_DATA(b)[i])
+                       return 1;
+
+       return 0;
+}
+
+/**
+ * set or clear all bits in the range [from;to[.
+ * @param a      The bitset.
+ * @param from   The first index to set to one.
+ * @param to     The last index plus one to set to one.
+ * @param do_set If 1 the bits are set, if 0, they are cleared.
+ */
+static inline void bitset_mod_range(bitset_t *a, bitset_pos_t from, bitset_pos_t to, int do_set)
+{
+       bitset_pos_t from_unit, to_unit, i;
+       bitset_unit_t from_unit_mask, to_unit_mask;
+
+       if (from == to)
+           return;
+
+       if (to < from) {
+               bitset_pos_t tmp = from;
+               from = to;
+               to = tmp;
+       }
+
+       if (to > a->size)
+               to = a->size;
+
+       /*
+        * A small example (for cleaning bits in the same unit).
+        * from   = 7
+        * to     = 19
+        * do_set = 0
+        * result:         xxxxxxx000000000000xxxxxxxxxxxxx
+        * from_unit_mask: 00000001111111111111111111111111
+        * to_unit_mask:   11111111111111111110000000000000
+        * scale:          01234567890123456789012345678901
+        *                           1         2         3
+        */
+
+       from_unit      = from / BS_UNIT_SIZE_BITS;
+       from_unit_mask = ~((1 << from) - 1);
+       from           = from & BS_UNIT_MASK;
+
+       to_unit        = to / BS_UNIT_SIZE_BITS;
+       to_unit_mask   = (1 << to) - 1;
+       to             = to & BS_UNIT_MASK;
+
+       /* do we want to set the bits in the range? */
+       if (do_set) {
+               /* If from and to are in the same unit: */
+               if (from_unit == to_unit)
+                       BS_DATA(a)[from_unit] |= from_unit_mask & to_unit_mask;
+
+               /* Else, we have to treat the from and to units seperately */
+               else {
+                       BS_DATA(a)[from_unit] |= from_unit_mask;
+                       BS_DATA(a)[to_unit]   |= to_unit_mask;
+                       for (i = from_unit + 1; i < to_unit; i++)
+                               BS_DATA(a)[i] = BITSET_UNIT_ALL_ONE;
+               }
+       }
+
+       /* ... or clear them? */
+       else {
+               if (from_unit == to_unit)
+                       BS_DATA(a)[from_unit] &= ~(from_unit_mask & to_unit_mask);
+
+               else {
+                       BS_DATA(a)[from_unit] &= ~from_unit_mask;
+                       BS_DATA(a)[to_unit]   &= ~to_unit_mask;
+                       for (i = from_unit + 1; i < to_unit; i++)
+                               BS_DATA(a)[i] = 0;
+               }
+       }
+}
+
+#define bitset_set_range(bs, from, to)   bitset_mod_range((bs), (from), (to), 1)
+#define bitset_clear_range(bs, from, to) bitset_mod_range((bs), (from), (to), 0)
+
+/**
+ * Check, if a bitset is empty.
+ * @param a The bitset.
+ * @return 1, if the bitset is empty, 0 if not.
+ */
+static inline int bitset_is_empty(const bitset_t *a)
+{
+       bitset_pos_t i;
+       for (i = 0; i < a->units; ++i)
+               if (BS_DATA(a)[i] != 0)
+                       return 0;
+       return 1;
+}
+
+/**
+ * Print a bitset to a stream.
+ * The bitset is printed as a comma separated list of bits set.
+ * @param file The stream.
+ * @param bs The bitset.
+ */
+static inline void bitset_fprint(FILE *file, const bitset_t *bs)
+{
+       const char *prefix = "";
+       int i;
+
+       putc('{', file);
+       for(i = bitset_next_set(bs, 0); i != -1; i = bitset_next_set(bs, i + 1)) {
+               fprintf(file, "%s%u", prefix, i);
+               prefix = ",";
+       }
+       putc('}', file);
+}
+
+static inline void bitset_debug_fprint(FILE *file, const bitset_t *bs)
+{
+       bitset_pos_t i;
+
+       fprintf(file, "%u:", bs->units);
+       for(i = 0; i < bs->units; ++i)
+               fprintf(file, " " BITSET_UNIT_FMT, BS_DATA(bs)[i]);
+}
+
+/**
+ * Perform tgt = tgt \ src operation.
+ * @param tgt  The target bitset.
+ * @param src  The source bitset.
+ * @return the tgt set.
+ */
+static inline bitset_t *bitset_andnot(bitset_t *tgt, const bitset_t *src);
+
+/**
+ * Perform Union, tgt = tgt u src operation.
+ * @param tgt  The target bitset.
+ * @param src  The source bitset.
+ * @return the tgt set.
+ */
+static inline bitset_t *bitset_or(bitset_t *tgt, const bitset_t *src);
+
+/**
+ * Perform tgt = tgt ^ ~src operation.
+ * @param tgt  The target bitset.
+ * @param src  The source bitset.
+ * @return the tgt set.
+ */
+static inline bitset_t *bitset_xor(bitset_t *tgt, const bitset_t *src);
+
 /*
  * Here, the binary operations follow.
  * And, Or, And Not, Xor are available.
  */
-
 #define BINARY_OP(op) \
-static INLINE bitset_t *bitset_ ## op(bitset_t *tgt, bitset_t *src) \
+static inline bitset_t *bitset_ ## op(bitset_t *tgt, const bitset_t *src) \
 { \
-       int i; \
-       int n = tgt->units > src->units ? src->units : tgt->units; \
-       bitset_unit_t *tgt_unit, *src_unit; \
-       src_unit = src->data; \
-       tgt_unit = tgt->data; \
-       for(i = 0; i < n; ++i) \
-               _bitset_inside_ ## op(tgt_unit++, src_unit++); \
-       return tgt; \
+       bitset_pos_t i; \
+       bitset_pos_t n = tgt->units > src->units ? src->units : tgt->units; \
+       for(i = 0; i < n; i += _BITSET_BINOP_UNITS_INC) \
+               _bitset_inside_binop_ ## op(&BS_DATA(tgt)[i], &BS_DATA(src)[i]); \
+       if(n < tgt->units) \
+               _bitset_clear_rest(&BS_DATA(tgt)[i], tgt->units - i); \
+       return _bitset_mask_highest(tgt); \
 }
 
+/*
+ * Define the clear rest macro for the and, since it is the only case,
+ * were non existed (treated as 0) units in the src must be handled.
+ * For all other operations holds: x Op 0 = x for Op in { Andnot, Or, Xor }
+ *
+ * For and, each bitset implementer has to provide the macro
+ * _bitset_clear_units(data, n), which clears n units from the pointer
+ * data on.
+ */
+#define _bitset_clear_rest(data,n) _bitset_inside_clear_units(data, n)
 BINARY_OP(and)
+#undef _bitset_clear_rest
+#define _bitset_clear_rest(data,n) do { } while(0)
+
 BINARY_OP(andnot)
 BINARY_OP(or)
 BINARY_OP(xor)
 
-
 #endif