sigalign-core 0.2.1

A core crate for sigalign
Documentation 
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

use super::{Anchor, AnchorTable, TraversedAnchor};

#[inline]
pub fn transform_right_additive_positions_to_traversed_anchor_index(
    anchor_table: &AnchorTable,
    traversed_anchors_buffer: &mut Vec<TraversedAnchor>,
    base_pattern_index: u32,
    base_target_position: u32,
    pattern_size: u32,
) {
    traversed_anchors_buffer.iter_mut().for_each(|tv| {
        let mut pattern_index = base_pattern_index + tv.addt_pattern_index;
        let mut target_position = base_target_position + tv.addt_target_position;
        let anchor_index_in_pattern = {
            loop {
                let anchors_by_pattern = &anchor_table.0[pattern_index as usize];
                match binary_search(
                    anchors_by_pattern,
                    target_position,
                ) {
                    Ok(v) => {
                        break v as u32
                    },
                    Err(_) => {
                        pattern_index -= 1;
                        target_position -= pattern_size;
                    },
                }
            }
        };
        tv.addt_pattern_index = pattern_index;
        tv.addt_target_position = anchor_index_in_pattern;
    });
}

#[inline(always)]
fn binary_search(
    anchors_by_pattern: &Vec<Anchor>,
    target_position: u32,
) -> Result<usize, usize> {
    anchors_by_pattern.binary_search_by_key(&target_position, |anchor| {
        anchor.target_position
    })
}