gtars-core 0.5.7

Core library for gtars: tools for high performance genomic interval analysis
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

#![allow(dead_code)]

use gtars_core::models::{Region, RegionSet};
use std::collections::HashMap;
use std::path::PathBuf;

pub type Iv<'a> = (&'a str, u32, u32);

pub fn make_region(chr: &str, start: u32, end: u32) -> Region {
    Region {
        chr: chr.to_string(),
        start,
        end,
        rest: None,
    }
}

pub fn make_regionset(regions: Vec<Iv>) -> RegionSet {
    RegionSet::from(
        regions
            .into_iter()
            .map(|(chr, start, end)| make_region(chr, start, end))
            .collect::<Vec<Region>>(),
    )
}

pub fn coords(rs: &RegionSet) -> Vec<Iv<'_>> {
    rs.regions
        .iter()
        .map(|r| (r.chr.as_str(), r.start, r.end))
        .collect()
}

pub fn chrom_sizes(sizes: &[(&str, u32)]) -> HashMap<String, u32> {
    sizes.iter().map(|(c, s)| (c.to_string(), *s)).collect()
}

pub fn load_bed(file_name: &str) -> RegionSet {
    let path = PathBuf::from("../tests/data/regionset").join(file_name);
    RegionSet::try_from(path.as_path()).unwrap()
}

pub fn assert_approx(actual: f64, expected: f64) {
    assert!(
        (actual - expected).abs() < 1e-10,
        "expected {expected}, got {actual}"
    );
}