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
pub mod alignment; pub use alignment::*; use std::cmp::Reverse; use std::collections::BinaryHeap; pub struct DepthIter<'a, R: AlignmentReader<'a>> { iter: AlignmentIter<'a, R>, cur_chrom: i32, cur_pos: usize, heap: BinaryHeap<Reverse<usize>>, next_read: Option<(i32, usize, usize)>, filter: Option<Box<dyn Fn(&Alignment) -> bool + 'a>>, } impl<'a, R: AlignmentReader<'a>> DepthIter<'a, R> { pub fn with_filter<F: Fn(&Alignment) -> bool + 'a>(reader: R, filter: F) -> Self { let (chrom, pos) = reader.start(); let iter = reader.to_alignment_iter(); let mut ret = Self { iter, next_read: None, cur_chrom: chrom as i32, cur_pos: pos as usize, heap: BinaryHeap::new(), filter: Some(Box::new(filter)), }; ret.load_next(); ret } fn load_next(&mut self) { self.next_read = loop { if let Some(Ok(read)) = self.iter.next() { if self.filter.as_ref().map_or(true, |predict| predict(&read)) { break Some(read); } } else { break None; } } .map(|read| (read.ref_id(), read.ref_begin(), read.ref_end())); } } impl<'a, R: AlignmentReader<'a>> Iterator for DepthIter<'a, R> { type Item = (i32, usize, u32); fn next(&mut self) -> Option<Self::Item> { if self.next_read.is_none() && self.heap.len() == 0 { return None; } let ret = (self.cur_chrom, self.cur_pos, self.heap.len() as u32); self.cur_pos += 1; while let Some((tid, left, right)) = self.next_read { if tid != self.cur_chrom { if self.heap.len() == 0 { self.cur_chrom = tid; self.cur_pos = 0; } break; } if left > self.cur_pos { break; } self.heap.push(Reverse(right)); self.next_read = self.iter.next().map(|x| { ( x.as_ref().unwrap().ref_id(), x.as_ref().unwrap().ref_begin(), x.as_ref().unwrap().ref_end(), ) }); } while self.heap.peek().map_or(false, |x| x.0 < self.cur_pos) { self.heap.pop(); } return Some(ret); } } #[cfg(test)] mod tests { #[test] fn it_works() { assert_eq!(2 + 2, 4); } }