pub struct MultiZip<I: Iterator>(Vec<I>);
impl<I: Iterator> MultiZip<I> {
pub fn new(vec_of_iters: Vec<I>) -> Self {
Self(vec_of_iters)
}
}
impl<I: Iterator> Iterator for MultiZip<I> {
type Item = Vec<I::Item>;
fn next(&mut self) -> Option<Self::Item> {
self.0.iter_mut().map(Iterator::next).collect()
}
}
pub trait WrappingGet<T> {
fn wrapping_next(&self, idx: usize) -> &T;
fn wrapping_prev(&self, idx: usize) -> &T;
}
impl<T> WrappingGet<T> for &[T] {
fn wrapping_next(&self, idx: usize) -> &T {
&self[match idx {
_ if idx == self.len() - 1 => 0,
_ => idx + 1,
}]
}
fn wrapping_prev(&self, idx: usize) -> &T {
&self[match idx {
_ if idx == 0 => self.len() - 1,
_ => idx - 1,
}]
}
}
#[derive(Clone, Copy, Debug)]
pub struct FloatComparator {
rel_tol: f64,
abs_tol: f64,
}
impl FloatComparator {
pub fn new(rel_tol: f64, abs_tol: f64) -> Self {
Self { rel_tol, abs_tol }
}
#[inline]
pub fn isclose(self, a: f64, b: f64) -> bool {
if a == b {
return true;
}
if !a.is_finite() || !b.is_finite() {
return false;
}
let diff = (a - b).abs();
diff <= (self.rel_tol * b).abs() || diff <= (self.rel_tol * a).abs() || diff <= self.abs_tol
}
}
impl Default for FloatComparator {
fn default() -> Self {
Self::new(1e-9, 0.0)
}
}
pub fn isclose(a: f64, b: f64) -> bool {
FloatComparator::default().isclose(a, b)
}
#[derive(Clone, Copy, Debug)]
pub struct SearchRange {
pub search_range: u16,
pub entry_selector: u16,
pub range_shift: u16,
}
impl SearchRange {
pub fn compute(n_items: usize, item_size: usize) -> Self {
let entry_selector = (n_items as f64).log2().floor() as usize;
let search_range = (2.0_f64.powi(entry_selector as i32) * item_size as f64) as usize;
let range_shift = (n_items * item_size).saturating_sub(search_range);
SearchRange {
search_range: search_range.try_into().unwrap(),
entry_selector: entry_selector.try_into().unwrap(),
range_shift: range_shift.try_into().unwrap(),
}
}
}
#[cfg(test)]
mod tests {
use super::*;
#[test]
fn simple_search_range() {
let SearchRange {
search_range,
entry_selector,
range_shift,
} = SearchRange::compute(39, 2);
assert_eq!((search_range, entry_selector, range_shift), (64, 5, 14));
}
#[test]
fn search_range_no_crashy() {
let foo = SearchRange::compute(0, 0);
assert_eq!(
(foo.search_range, foo.entry_selector, foo.range_shift),
(0, 0, 0)
)
}
}