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

use Source;

/// The Xorshift+ algorithm [1].
///
/// 1. https://en.wikipedia.org/wiki/Xorshift#Xorshift.2B
pub struct XorshiftPlus {
    state: [u64; 2],
}

impl XorshiftPlus {
    /// Create a generator.
    #[inline(always)]
    pub fn new(seed: [u64; 2]) -> XorshiftPlus {
        XorshiftPlus { state: seed }
    }
}

impl Source for XorshiftPlus {
    #[inline(always)]
    fn next(&mut self) -> u64 {
        let (mut x, y) = (self.state[0], self.state[1]);

        self.state[0] = y;
        x = x ^ (x << 23);
        x = x ^ (x >> 17);
        x = x ^ y ^ (y >> 26);
        self.state[1] = x;

        x.wrapping_add(y)
    }
}