#![cfg_attr(not(feature = "std"), no_std)]
#[cfg(feature = "std")]
use std::ops::{BitAnd, BitOr};
#[cfg(not(feature = "std"))]
use core::ops::{BitAnd, BitOr};
pub const MAX_DEPTH: usize = 17;
fn sub_mix<T>(seeds: &[T]) -> T where
T: BitAnd<Output = T> + BitOr<Output = T> + Copy
{
(seeds[0] & seeds[1]) | (seeds[1] & seeds[2]) | (seeds[0] & seeds[2])
}
pub fn triplet_mix<T>(seeds: &[T]) -> Result<T, ()> where
T: BitAnd<Output = T> + BitOr<Output = T>,
T: Default + Copy
{
Ok(seeds.iter().cloned().triplet_mix())
}
pub trait TripletMix {
type Item;
fn triplet_mix(self) -> Self::Item;
}
impl<I, T> TripletMix for I where
I: Iterator<Item = T>,
T: BitAnd<Output = T> + BitOr<Output = T> + Default + Copy
{
type Item = T;
fn triplet_mix(self) -> Self::Item {
let mut accum = [[T::default(); 3]; MAX_DEPTH];
let mut result = T::default();
for (i, seed) in self.enumerate() {
accum[0][i % 3] = seed;
let mut index_at_depth = i;
for depth in 0..MAX_DEPTH {
if index_at_depth % 3 != 2 {
break;
}
index_at_depth /= 3;
result = sub_mix(&accum[depth]);
if depth == MAX_DEPTH - 1 {
break;
} else {
accum[depth + 1][index_at_depth % 3] = result;
}
}
}
result
}
}
#[cfg(test)]
mod tests {
use super::*;
#[test]
fn sub_mix_works() {
assert_eq!(sub_mix(&[0, 0, 0][..]), 0);
assert_eq!(sub_mix(&[0, 0, 1][..]), 0);
assert_eq!(sub_mix(&[0, 1, 0][..]), 0);
assert_eq!(sub_mix(&[0, 1, 1][..]), 1);
assert_eq!(sub_mix(&[1, 0, 0][..]), 0);
assert_eq!(sub_mix(&[1, 0, 1][..]), 1);
assert_eq!(sub_mix(&[1, 1, 0][..]), 1);
assert_eq!(sub_mix(&[1, 1, 1][..]), 1);
assert_eq!(sub_mix(&[0, 0, 0][..]), 0);
assert_eq!(sub_mix(&[0, 0, 2][..]), 0);
assert_eq!(sub_mix(&[0, 2, 0][..]), 0);
assert_eq!(sub_mix(&[0, 2, 2][..]), 2);
assert_eq!(sub_mix(&[2, 0, 0][..]), 0);
assert_eq!(sub_mix(&[2, 0, 2][..]), 2);
assert_eq!(sub_mix(&[2, 2, 0][..]), 2);
assert_eq!(sub_mix(&[2, 2, 2][..]), 2);
}
#[test]
fn triplet_mix_works_on_first_level() {
assert_eq!(triplet_mix(&[0, 0, 0][..]).unwrap(), 0);
assert_eq!(triplet_mix(&[0, 0, 1][..]).unwrap(), 0);
assert_eq!(triplet_mix(&[0, 1, 0][..]).unwrap(), 0);
assert_eq!(triplet_mix(&[0, 1, 1][..]).unwrap(), 1);
assert_eq!(triplet_mix(&[1, 0, 0][..]).unwrap(), 0);
assert_eq!(triplet_mix(&[1, 0, 1][..]).unwrap(), 1);
assert_eq!(triplet_mix(&[1, 1, 0][..]).unwrap(), 1);
assert_eq!(triplet_mix(&[1, 1, 1][..]).unwrap(), 1);
assert_eq!(triplet_mix(&[0, 0, 0][..]).unwrap(), 0);
assert_eq!(triplet_mix(&[0, 0, 2][..]).unwrap(), 0);
assert_eq!(triplet_mix(&[0, 2, 0][..]).unwrap(), 0);
assert_eq!(triplet_mix(&[0, 2, 2][..]).unwrap(), 2);
assert_eq!(triplet_mix(&[2, 0, 0][..]).unwrap(), 0);
assert_eq!(triplet_mix(&[2, 0, 2][..]).unwrap(), 2);
assert_eq!(triplet_mix(&[2, 2, 0][..]).unwrap(), 2);
assert_eq!(triplet_mix(&[2, 2, 2][..]).unwrap(), 2);
}
#[test]
fn triplet_mix_works_on_second_level() {
assert_eq!(triplet_mix(&[0, 0, 0, 0, 0, 1, 0, 1, 0][..]).unwrap(), 0);
assert_eq!(triplet_mix(&[0, 1, 1, 1, 0, 0, 1, 0, 1][..]).unwrap(), 1);
assert_eq!(triplet_mix(&[1, 1, 0, 1, 1, 1, 0, 0, 0][..]).unwrap(), 1);
}
#[test]
fn triplet_mix_works_on_third_level() {
assert_eq!(triplet_mix(&[0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 0, 0][..]).unwrap(), 1);
}
}