quadrank 0.2.0

Fast rank over binary and size-4 DNA alphabets.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

use crate::binary::BiRanker;

pub use bio::data_structures::rank_select::RankSelect;

impl BiRanker for RankSelect {
    fn new_packed(seq: &[u64]) -> Self {
        let seq_u8 = unsafe { seq.align_to::<u8>().1 };
        let bitvec = bv::BitVec::<u8>::from_bits(seq_u8);
        let k = 512 / 32;
        Self::new(bitvec, k)
    }

    fn size(&self) -> usize {
        self.bits().len() as usize / 8 + std::mem::size_of_val(self.superblocks_1.as_slice())
    }

    unsafe fn rank_unchecked(&self, pos: usize) -> u64 {
        self.rank_1(pos as u64 - 1).unwrap()
    }
}