1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
use crate::bits::{BitBlock, BitSetView, MAX_LEVEL};
use num_traits::{One, Zero};

/// Iterator over the set bits.
pub struct BitIter<B: BitSetView> {
    bitset: B,
    // masked block for each level, consumed bits are cleared
    masks: [B::Bits; MAX_LEVEL],
    // index prefix for each level
    prefixes: [usize; MAX_LEVEL],
}

impl<B: BitSetView> BitIter<B> {
    pub fn new(bitset: B) -> BitIter<B> {
        let mut iter = BitIter {
            bitset,
            masks: [B::Bits::zero(); MAX_LEVEL],
            prefixes: [0; MAX_LEVEL],
        };

        //init mask to perform a full descend when performing the first step
        let last_level = iter.bitset.get_level_count() - 1;
        let top = iter.bitset.get_block(last_level, 0);
        iter.masks[last_level] = top;
        iter
    }

    fn step(&mut self) -> Option<usize> {
        let lc = self.bitset.get_level_count();
        let mut level = 0;
        loop {
            while self.masks[level].is_zero() {
                // no bits in this block, move upward
                level += 1;
                if level >= lc {
                    // top reached and no more bits were found
                    return None;
                }
            }

            // take next set bit
            let prefix = {
                let block = &mut self.masks[level];
                let pos = block.trailing_bit_pos();
                *block = *block ^ (B::Bits::one() << pos);
                (self.prefixes[level] << B::Bits::bit_shift()) | pos
            };

            // move downward
            if level == 0 {
                // bottom reached, prefix is the index of the set bit
                return Some(prefix);
            }
            level -= 1;
            self.masks[level] = self.bitset.get_block(level, prefix);
            self.prefixes[level] = prefix;
        }
    }
}

impl<B: BitSetView> Iterator for BitIter<B> {
    type Item = usize;

    fn next(&mut self) -> Option<Self::Item> {
        self.step()
    }
}