int-interval-stack 0.3.0

Integer half-open interval stack structures for overlap multiplicity.
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
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75

use rayon::iter::ParallelIterator;

use crate::int_co_stack::test_support::{iv_i32, stack_from_intervals};

use super::*;

fn collect_runs(window: StackWindow<'_, int_interval::I32CO>) -> Vec<((i32, i32), usize)> {
    window
        .iter_height_runs()
        .map(|run| ((run.interval.start(), run.interval.end_excl()), run.height))
        .collect()
}

fn collect_runs_parallel(window: StackWindow<'_, int_interval::I32CO>) -> Vec<((i32, i32), usize)> {
    window
        .par_iter_height_runs()
        .map(|run| ((run.interval.start(), run.interval.end_excl()), run.height))
        .collect()
}

#[test]
fn serial_and_parallel_height_runs_match_for_single_run() {
    let stack = stack_from_intervals(&[(2, 8)]);
    let window = StackWindow::new(&stack, iv_i32(3, 6));

    assert_eq!(collect_runs(window), collect_runs_parallel(window));
}

#[test]
fn serial_and_parallel_height_runs_match_for_multiple_runs() {
    let stack = stack_from_intervals(&[(1, 4), (3, 6)]);
    let window = StackWindow::new(&stack, iv_i32(2, 5));

    assert_eq!(
        collect_runs(window),
        vec![((2, 3), 1), ((3, 4), 2), ((4, 5), 1)],
    );

    assert_eq!(collect_runs(window), collect_runs_parallel(window));
}

#[test]
fn serial_and_parallel_height_runs_match_with_zero_gap() {
    let stack = stack_from_intervals(&[(1, 3), (5, 7)]);
    let window = StackWindow::new(&stack, iv_i32(2, 6));

    assert_eq!(
        collect_runs(window),
        vec![((2, 3), 1), ((3, 5), 0), ((5, 6), 1)],
    );

    assert_eq!(collect_runs(window), collect_runs_parallel(window));
}

#[test]
fn serial_and_parallel_height_runs_match_after_all_change_points() {
    let stack = stack_from_intervals(&[(1, 3)]);
    let window = StackWindow::new(&stack, iv_i32(3, 6));

    assert_eq!(collect_runs(window), vec![((3, 6), 0)]);
    assert_eq!(collect_runs(window), collect_runs_parallel(window));
}

#[test]
fn serial_and_parallel_height_runs_match_when_window_edges_are_change_points() {
    let stack = stack_from_intervals(&[(1, 4), (3, 6)]);
    let window = StackWindow::new(&stack, iv_i32(1, 6));

    assert_eq!(
        collect_runs(window),
        vec![((1, 3), 1), ((3, 4), 2), ((4, 6), 1)],
    );

    assert_eq!(collect_runs(window), collect_runs_parallel(window));
}