rule30 0.2.1

use rand_core::impls::fill_bytes_via_next;
use rand_core::le::read_u64_into;
use rand_core::{RngCore, SeedableRng};

// Default cells number is 80 x 64 = 5120.
const SIZE: usize = 80;
const SKIP: usize = 4;
const INDEX: [usize; SKIP] = [16, 32, 48, 54];
const MASK: u64 = 0b0001_0001_0001_0001_0001_0001_0001_0001_0001_0001_0001_0001_0001_0001_0001_0001;

/// Pseudo-random number generator using Extended CA.
#[derive(Clone, Eq, PartialEq, Debug)]
pub struct ExtendedCA {
    state: [u64; SIZE],
}

impl ExtendedCA {
    /// Create a new `ExtendedCA` instance with the given state.
    pub fn new(state: [u64; SIZE]) -> Self {
        ExtendedCA { state }
    }

    /// Clone the inner state.
    pub fn state(&self) -> [u64; SIZE] {
        self.state
    }

    /// Reset with the provided state.
    pub fn reset(&mut self, state: [u64; SIZE]) {
        self.state = state;
    }

    /// Extend the state array with 1 more state on both sides to make the evolution code simple.
    #[inline]
    fn extend_state(&self) -> [u64; SIZE + 2] {
        let mut extend = [0u64; SIZE + 2];
        extend[1..(SIZE + 1)].copy_from_slice(&self.state);
        extend[0] = self.state[SIZE - 1];
        extend[SIZE + 1] = self.state[0];
        extend
    }

    /// From the Wolfram reference on Extended CA:
    ///
    /// > "This generator uses a particular five-neighbor rule, so each new
    /// > cell depends on five nonadjacent cells from the previous step."
    ///
    /// It doesn't mention which particular five-neighbor rule is used. So I
    /// would choose 5 neighbors at the same position in the adjacent states.
    /// The evolving rule is similar to the Rule 30 (3 neighbors):
    ///
    /// ```text
    /// next_state[i] =
    ///     (state[i-1] | (state[i-1] <<. 3)) ^ state[i] ^ (state[i+1] | (state[i+1] <<. 3))
    /// ```
    ///
    /// Any odd number step rotation can be used. "3" is chosen since it requires
    /// less round to spread to all the cells. While "1" might break the "nonadjacent"
    /// rule.
    ///
    /// This requires further investigation.
    #[inline]
    fn step(&mut self) {
        let extend = self.extend_state();
        for i in 0..SIZE {
            self.state[i] = (extend[i] | extend[i].rotate_left(3))
                ^ extend[i + 1]
                ^ (extend[i + 2] | extend[i + 2].rotate_left(3));
        }
    }
}

#[derive(Clone, Debug)]
pub struct ExtendedRngSeed(pub [u8; SIZE * 8]);

impl Default for ExtendedRngSeed {
    fn default() -> Self {
        ExtendedRngSeed([0u8; SIZE * 8])
    }
}

impl AsMut<[u8]> for ExtendedRngSeed {
    fn as_mut(&mut self) -> &mut [u8] {
        &mut self.0
    }
}

impl AsRef<[u8]> for ExtendedRngSeed {
    fn as_ref(&self) -> &[u8] {
        &self.0
    }
}

impl SeedableRng for ExtendedCA {
    type Seed = ExtendedRngSeed;

    fn from_seed(seed: Self::Seed) -> Self {
        let mut seed_u64 = [0u64; SIZE];
        read_u64_into(seed.as_ref(), &mut seed_u64);
        ExtendedCA::new(seed_u64)
    }
}

impl RngCore for ExtendedCA {
    #[inline]
    fn next_u32(&mut self) -> u32 {
        self.next_u64() as u32
    }

    #[inline]
    fn next_u64(&mut self) -> u64 {
        let mut num: u64 = 0;
        for i in INDEX.iter() {
            num = (num << 1) | (self.state[*i] & MASK);
        }
        self.step();
        num
    }

    #[inline]
    fn fill_bytes(&mut self, dest: &mut [u8]) {
        fill_bytes_via_next(self, dest)
    }
}

#[cfg(test)]
mod tests {
    use rand_core::{RngCore, SeedableRng};

    use super::{ExtendedCA, ExtendedRngSeed, SIZE};

    #[test]
    fn test_extended_ca_construction() {
        let mut rng_u64 = ExtendedCA::seed_from_u64(42);
        assert_eq!(rng_u64.next_u64(), 12588493861803088380);
        assert_eq!(rng_u64.next_u64(), 13092609834719595470);

        let mut rng_u32 = ExtendedCA::seed_from_u64(42);
        assert_eq!(rng_u32.next_u32(), 2046122492);
        assert_eq!(rng_u32.next_u32(), 3980348366);

        let mut rng_gen = ExtendedCA::from_rng(&mut rng_u64);
        assert_eq!(rng_gen.next_u64(), 9131884500388100982);

        let seed: [u8; SIZE * 8] = core::array::from_fn(|x| x as u8);
        let mut rng_seed = ExtendedCA::from_seed(ExtendedRngSeed(seed));
        assert_eq!(rng_seed.next_u64(), 2238323168904224528);

        let mut rng_double_check = ExtendedCA::seed_from_u64(42);
        assert_eq!(rng_double_check.next_u64(), 12588493861803088380);
    }
}