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
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
use rand_core::impls::fill_bytes_via_next;
use rand_core::le::read_u64_into;
use rand_core::{SeedableRng, RngCore, Error};
use Seed512;
#[derive(Debug, Clone)]
pub struct Xoshiro512Plus {
s: [u64; 8],
}
impl Xoshiro512Plus {
pub fn from_seed_u64(seed: u64) -> Xoshiro512Plus {
from_splitmix!(seed)
}
pub fn jump(&mut self) {
impl_jump!(u64, self, [
0x33ed89b6e7a353f9, 0x760083d7955323be, 0x2837f2fbb5f22fae,
0x4b8c5674d309511c, 0xb11ac47a7ba28c25, 0xf1be7667092bcc1c,
0x53851efdb6df0aaf, 0x1ebbc8b23eaf25db
]);
}
}
impl SeedableRng for Xoshiro512Plus {
type Seed = Seed512;
#[inline]
fn from_seed(seed: Seed512) -> Xoshiro512Plus {
deal_with_zero_seed!(seed, Self);
let mut state = [0; 8];
read_u64_into(&seed.0, &mut state);
Xoshiro512Plus { s: state }
}
}
impl RngCore for Xoshiro512Plus {
#[inline]
fn next_u32(&mut self) -> u32 {
self.next_u64() as u32
}
#[inline]
fn next_u64(&mut self) -> u64 {
let result_plus = self.s[0].wrapping_add(self.s[2]);
impl_xoshiro_large!(self);
result_plus
}
#[inline]
fn fill_bytes(&mut self, dest: &mut [u8]) {
fill_bytes_via_next(self, dest);
}
#[inline]
fn try_fill_bytes(&mut self, dest: &mut [u8]) -> Result<(), Error> {
self.fill_bytes(dest);
Ok(())
}
}
#[cfg(test)]
mod tests {
use super::*;
#[test]
fn reference() {
let mut rng = Xoshiro512Plus::from_seed(Seed512(
[1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0,
3, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0,
5, 0, 0, 0, 0, 0, 0, 0, 6, 0, 0, 0, 0, 0, 0, 0,
7, 0, 0, 0, 0, 0, 0, 0, 8, 0, 0, 0, 0, 0, 0, 0]));
let expected = [
4, 8, 4113, 25169936, 52776585412635, 57174648719367,
9223482039571869716, 9331471677901559830, 9340533895746033672,
14078399799840753678,
];
for &e in &expected {
assert_eq!(rng.next_u64(), e);
}
}
}