fibertools-rs 0.12.0

Fiber-seq toolkit in rust
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

use crate::cli::BenchmarkOptions;
use anyhow::Result;
use std::time::Instant;

pub fn run_benchmark(opts: &mut BenchmarkOptions) -> Result<()> {
    log::info!("Starting fiber iterator benchmark");

    let mut bam = opts.input.bam_reader();

    // Benchmark input.fibers()
    let start_time = Instant::now();
    let mut fiber_count = 0;

    for fiber in opts.input.fibers(&mut bam) {
        fiber_count += 1;
        // Force materialization but do nothing with the fiber
        std::hint::black_box(&fiber);
    }

    let elapsed = start_time.elapsed();

    log::info!(
        "Processed {} fibers in {:.3} seconds",
        fiber_count,
        elapsed.as_secs_f64()
    );
    log::info!(
        "Rate: {:.0} fibers/second",
        fiber_count as f64 / elapsed.as_secs_f64()
    );

    Ok(())
}