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 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119
//! Prints a progress meter for number of fastq entries to stderr.
//!
//! # Examples
//!
//! ## fasten_metrics progress
//! While getting read metrics for a large fastq file, print the progress
//! to make the wait a little easier
//! ```bash
//! cat large.fastq | fasten_progress --print | fasten_metrics
//! ```
//! ## fasten_shuffle progress
//! While shuffling a large fastq file, print the progress
//! ```bash
//! cat large_1.fastq large_2.fastq | fasten_progress --print | fasten_shuffle > interleaved.fastq
//! ```
//! ## Two progress bars
//! When there is a halting step in the process like `fasten_sort`, then it might make sense
//! to have two progress bars.
//! However, if there are no halting steps then the progress messages will collide.
//! ```bash
//! cat large_1.fastq large_2.fastq | \
//! fasten_progress --id first --print | \
//! fasten_shuffle | \
//! fasten_sort --paired-end | \
//! fasten_progress --id second --print | \
//! fasten_progress --id collision-with-second --print | \
//! fasten_metrics | column -t
//! ```
//!
//! # Usage
//!
//! ```text
//! Usage: fasten_progress [-h] [-n INT] [-p] [-v] [--id STRING] [--update-every INT]
//!
//! Options:
//! -h, --help Print this help menu.
//! -n, --numcpus INT Number of CPUs (default: 1)
//! -p, --paired-end The input reads are interleaved paired-end
//! -v, --verbose Print more status messages
//! --id STRING Progress identifier. Default: unnamed
//! --update-every INT
//! Update progress every n reads.
//! -p, --print Print the reads back to stdout
//! ```
extern crate getopts;
extern crate fasten;
use std::fs::File;
use std::io::BufReader;
use std::io::BufRead;
use fasten::fasten_base_options;
use fasten::fasten_base_options_matches;
fn main(){
let mut opts = fasten_base_options();
// Options specific to this script
opts.optopt("","id","Progress identifier. Default: unnamed","STRING");
opts.optopt("","update-every","Update progress every n reads.","INT");
opts.optflag("p","print","Print the reads back to stdout");
let matches = fasten_base_options_matches("Prints a progress meter for number of fastq entries", opts);
let print_reads:bool = matches.opt_present("print");
let progress_id:String = {
if matches.opt_present("id") {
matches.opt_str("id")
.expect("ERROR parsing --id")
} else {
String::from("unnamed")
}
};
let my_file = File::open("/dev/stdin").expect("Could not open file");
let my_buffer=BufReader::new(my_file);
let lines_per_read :usize = 4;
let update_every :usize = {
if matches.opt_present("update-every") {
let tmp :usize =
matches.opt_str("update-every")
.expect("ERROR parsing update-every")
.parse()
.expect("ERROR parsing update-every as INT");
tmp * lines_per_read
//tmp
} else {
100
}
};
let mut line_counter = 0;
eprint!("\r{} progress: {}", progress_id, line_counter/lines_per_read);
for res in my_buffer.lines() {
let line=res.expect("ERROR: did not get a line");
if print_reads {
println!("{}", line);
}
line_counter += 1;
match line_counter % update_every {
0=>{
//println!("UPDATE: {}", line_counter);
eprint!("\r{} progress: {}", progress_id, line_counter/lines_per_read);
//eprint!(".");
}
_=>{
}
}
}
eprint!("\n");
let msg = format!("{}: Finished progress on {} reads", progress_id, line_counter);
fasten::logmsg(&msg);
}