refactor: optimize dump partition iteration and add progress tracking #20

Merged
coissac merged 1 commits from push-xqswlxlvmyrq into main 2026-06-09 09:34:13 +00:00
2 changed files with 51 additions and 28 deletions
+45 -25
View File
@@ -21,13 +21,14 @@ impl KmerIndex {
///
/// The caller must have set the global kmer length (`obikseq::set_k`) before
/// calling this method.
pub fn dump<W: Write>(
pub fn dump<W: Write, F: Fn() + Send + Sync>(
&self,
out: &mut W,
force_presence: bool,
debug: bool,
head: Option<usize>,
filters: &[Box<dyn KmerFilter>],
on_partition: F,
) -> OKIResult<()> {
let genomes = &self.meta.genomes;
let use_counts = self.meta.config.with_counts && !force_presence;
@@ -46,51 +47,70 @@ impl KmerIndex {
// ── Rows — parallel over partitions ───────────────────────────────────
let n = self.n_partitions();
let remaining = AtomicUsize::new(head.unwrap_or(usize::MAX));
let chunks: Vec<OKIResult<Vec<u8>>> = (0..n)
.into_par_iter()
.map(|i| {
if remaining.load(Ordering::Relaxed) == 0 {
return Ok(vec![]);
}
let write_row = |buf: &mut Vec<u8>, row: &[u32], prefix: &str| {
let _ = buf.write_all(prefix.as_bytes());
for &v in row { let _ = write!(buf, ",{v}"); }
let _ = buf.write_all(b"\n");
};
let chunks: Vec<OKIResult<Vec<u8>>> = if let Some(limit) = head {
// ── Bounded: atomic counter, early exit when limit reached ────────
let remaining = AtomicUsize::new(limit);
(0..n).into_par_iter().map(|i| {
if remaining.load(Ordering::Relaxed) == 0 { return Ok(vec![]); }
let mut buf = Vec::<u8>::new();
let write_kmer = |buf: &mut Vec<u8>, row: &[u32], prefix: &str| -> bool {
let try_write = |buf: &mut Vec<u8>, row: &[u32], prefix: &str| -> bool {
match remaining.fetch_update(Ordering::SeqCst, Ordering::SeqCst, |cur| {
if cur > 0 { Some(cur - 1) } else { None }
}) {
Err(_) => false,
Ok(_) => {
let _ = buf.write_all(prefix.as_bytes());
for &v in row { let _ = write!(buf, ",{v}"); }
let _ = buf.write_all(b"\n");
true
}
Ok(_) => { write_row(buf, row, prefix); true }
}
};
if debug {
self.partition
.iter_partition_kmers_located(i, use_counts, n_genomes, filters, |part, layer, kmer, row| {
let seq = String::from_utf8(kmer.to_ascii())
.unwrap_or_else(|_| "?".repeat(kmer_size));
write_kmer(&mut buf, &row, &format!("{part},{layer},{seq}"))
let seq = String::from_utf8(kmer.to_ascii()).unwrap_or_else(|_| "?".repeat(kmer_size));
try_write(&mut buf, &row, &format!("{part},{layer},{seq}"))
})
.map_err(OKIError::Partition)?;
} else {
self.partition
.iter_partition_kmers(i, use_counts, n_genomes, filters, |kmer, row| {
let seq = String::from_utf8(kmer.to_ascii())
.unwrap_or_else(|_| "?".repeat(kmer_size));
write_kmer(&mut buf, &row, &seq)
let seq = String::from_utf8(kmer.to_ascii()).unwrap_or_else(|_| "?".repeat(kmer_size));
try_write(&mut buf, &row, &seq)
})
.map_err(OKIError::Partition)?;
}
on_partition();
Ok(buf)
}).collect()
} else {
// ── Unbounded: no atomic, no contention ───────────────────────────
(0..n).into_par_iter().map(|i| {
let mut buf = Vec::<u8>::new();
if debug {
self.partition
.iter_partition_kmers_located(i, use_counts, n_genomes, filters, |part, layer, kmer, row| {
let seq = String::from_utf8(kmer.to_ascii()).unwrap_or_else(|_| "?".repeat(kmer_size));
write_row(&mut buf, &row, &format!("{part},{layer},{seq}"));
true
})
.collect();
.map_err(OKIError::Partition)?;
} else {
self.partition
.iter_partition_kmers(i, use_counts, n_genomes, filters, |kmer, row| {
let seq = String::from_utf8(kmer.to_ascii()).unwrap_or_else(|_| "?".repeat(kmer_size));
write_row(&mut buf, &row, &seq);
true
})
.map_err(OKIError::Partition)?;
}
on_partition();
Ok(buf)
}).collect()
};
// ── Sequential write ──────────────────────────────────────────────────
for chunk in chunks {
+4 -1
View File
@@ -3,6 +3,7 @@ use std::path::PathBuf;
use clap::Args;
use obikindex::KmerIndex;
use obisys::progress_bar;
use tracing::info;
use super::predicate::FilterArgs;
@@ -41,12 +42,14 @@ pub fn run(args: DumpArgs) {
);
let filters = args.filter.build_filters(&idx.meta().genomes);
let pb = progress_bar("dump", idx.n_partitions() as u64, "partitions");
let stdout = io::stdout();
let mut out = BufWriter::new(stdout.lock());
idx.dump(&mut out, args.force_presence, args.debug, args.head, &filters).unwrap_or_else(|e| {
idx.dump(&mut out, args.force_presence, args.debug, args.head, &filters, || pb.inc(1)).unwrap_or_else(|e| {
eprintln!("dump error: {e}");
std::process::exit(1);
});
pb.finish_and_clear();
}