moved to romu_trio rand

This commit is contained in:
Dominik Maier 2020-12-09 22:30:31 +01:00
parent 0b8308135d
commit c893c053b1
5 changed files with 107 additions and 38 deletions

View File

@ -1,14 +1,22 @@
//! Compare the speed of rand implementations //! Compare the speed of rand implementations
use afl::utils::{Rand, XorShift64Rand, Xoshiro256StarRand}; use afl::utils::{
Lehmer64Rand, Rand, RomuDuoJrRand, RomuTrioRand, XorShift64Rand, Xoshiro256StarRand,
};
use criterion::{black_box, criterion_group, criterion_main, Criterion}; use criterion::{black_box, criterion_group, criterion_main, Criterion};
fn criterion_benchmark(c: &mut Criterion) { fn criterion_benchmark(c: &mut Criterion) {
let mut xorshift = XorShift64Rand::new(0); let mut xorshift = XorShift64Rand::new(1);
let mut xoshiro = Xoshiro256StarRand::new(0); let mut xoshiro = Xoshiro256StarRand::new(1);
let mut romu = RomuDuoJrRand::new(1);
let mut lehmer = Lehmer64Rand::new(1);
let mut romu_trio = RomuTrioRand::new(1);
c.bench_function("xorshift", |b| b.iter(|| black_box(xorshift.next()))); c.bench_function("xorshift", |b| b.iter(|| black_box(xorshift.next())));
c.bench_function("xoshiro", |b| b.iter(|| black_box(xoshiro.next()))); c.bench_function("xoshiro", |b| b.iter(|| black_box(xoshiro.next())));
c.bench_function("romu", |b| b.iter(|| black_box(romu.next())));
c.bench_function("romu_trio", |b| b.iter(|| black_box(romu_trio.next())));
c.bench_function("lehmer", |b| b.iter(|| black_box(lehmer.next())));
} }
criterion_group!(benches, criterion_benchmark); criterion_group!(benches, criterion_benchmark);

View File

@ -301,7 +301,8 @@ where
let cur = current_milliseconds(); let cur = current_milliseconds();
if cur - last > 60 * 100 { if cur - last > 60 * 100 {
last = cur; last = cur;
manager.fire(Event::update_stats(state.executions(), state.executions_over_seconds()), manager.fire(
Event::update_stats(state.executions(), state.executions_over_seconds()),
state, state,
corpus, corpus,
)?; // TODO self.new_execs}); )?; // TODO self.new_execs});

View File

@ -33,9 +33,7 @@ pub enum BrokerEventResult {
Forward, Forward,
} }
pub trait ShowStats { pub trait ShowStats {}
}
/* /*

View File

@ -300,4 +300,3 @@ where
} }
} }
*/ */

View File

@ -9,14 +9,16 @@ use xxhash_rust::xxh3::xxh3_64_with_seed;
#[cfg(feature = "std")] #[cfg(feature = "std")]
use std::time::{SystemTime, UNIX_EPOCH}; use std::time::{SystemTime, UNIX_EPOCH};
pub type StdRand = Xoshiro256StarRand; pub type StdRand = RomuTrioRand;
/// Ways to get random around here /// Ways to get random around here
pub trait Rand: Debug { pub trait Rand: Debug {
// Sets the seed of this Rand // Sets the seed of this Rand
fn set_seed(&mut self, seed: u64); fn set_seed(&mut self, seed: u64);
// Gets the next 64 bit value // Gets the next 64 bit value
fn next(&mut self) -> u64; fn next(&mut self) -> u64;
// Gets a value below the given 64 bit val (inclusive) // Gets a value below the given 64 bit val (inclusive)
fn below(&mut self, upper_bound_excl: u64) -> u64 { fn below(&mut self, upper_bound_excl: u64) -> u64 {
if upper_bound_excl <= 1 { if upper_bound_excl <= 1 {
@ -73,6 +75,15 @@ where
const HASH_CONST: u64 = 0xa5b35705; const HASH_CONST: u64 = 0xa5b35705;
#[cfg(feature = "std")]
/// Gets current nanoseconds since UNIX_EPOCH
pub fn current_nanos() -> u64 {
SystemTime::now()
.duration_since(UNIX_EPOCH)
.unwrap()
.as_nanos() as u64
}
/// XXH3 Based, hopefully speedy, rnd implementation /// XXH3 Based, hopefully speedy, rnd implementation
/// ///
#[derive(Copy, Clone, Debug, Default)] #[derive(Copy, Clone, Debug, Default)]
@ -91,6 +102,7 @@ impl Rand for Xoshiro256StarRand {
self.seeded = true; self.seeded = true;
} }
#[inline]
fn next(&mut self) -> u64 { fn next(&mut self) -> u64 {
let ret: u64 = self.rand_seed[0] let ret: u64 = self.rand_seed[0]
.wrapping_add(self.rand_seed[3]) .wrapping_add(self.rand_seed[3])
@ -125,19 +137,11 @@ impl Xoshiro256StarRand {
ret ret
} }
pub fn to_rc_refcell(self) -> Rc<RefCell<Self>> {
self.into()
}
/// Creates a rand instance, pre-seeded with the current time in nanoseconds. /// Creates a rand instance, pre-seeded with the current time in nanoseconds.
/// Needs stdlib timer /// Needs stdlib timer
#[cfg(feature = "std")] #[cfg(feature = "std")]
pub fn preseeded() -> Self { pub fn preseeded() -> Self {
let seed = SystemTime::now() Self::new(current_nanos())
.duration_since(UNIX_EPOCH)
.unwrap()
.as_nanos() as u64;
Self::new(seed)
} }
} }
@ -155,6 +159,7 @@ impl Rand for XorShift64Rand {
self.seeded = true; self.seeded = true;
} }
#[inline]
fn next(&mut self) -> u64 { fn next(&mut self) -> u64 {
let mut x = self.rand_seed; let mut x = self.rand_seed;
x ^= x << 13; x ^= x << 13;
@ -179,19 +184,11 @@ impl XorShift64Rand {
ret ret
} }
pub fn to_rc_refcell(self) -> Rc<RefCell<Self>> {
self.into()
}
/// Creates a rand instance, pre-seeded with the current time in nanoseconds. /// Creates a rand instance, pre-seeded with the current time in nanoseconds.
/// Needs stdlib timer /// Needs stdlib timer
#[cfg(feature = "std")] #[cfg(feature = "std")]
pub fn preseeded() -> Self { pub fn preseeded() -> Self {
let seed = SystemTime::now() Self::new(current_nanos())
.duration_since(UNIX_EPOCH)
.unwrap()
.as_nanos() as u64;
Self::new(seed)
} }
} }
@ -209,6 +206,7 @@ impl Rand for Lehmer64Rand {
self.seeded = true; self.seeded = true;
} }
#[inline]
fn next(&mut self) -> u64 { fn next(&mut self) -> u64 {
self.rand_seed *= 0xda942042e4dd58b5; self.rand_seed *= 0xda942042e4dd58b5;
return (self.rand_seed >> 64) as u64; return (self.rand_seed >> 64) as u64;
@ -222,26 +220,91 @@ impl Into<Rc<RefCell<Self>>> for Lehmer64Rand {
} }
impl Lehmer64Rand { impl Lehmer64Rand {
/// Creates a new Xoshiro rand with the given seed /// Creates a new Lehmer rand with the given seed
pub fn new(seed: u64) -> Self { pub fn new(seed: u64) -> Self {
let mut ret: Self = Default::default(); let mut ret: Self = Default::default();
ret.set_seed(seed); // TODO: Proper random seed? ret.set_seed(seed);
ret ret
} }
pub fn to_rc_refcell(self) -> Rc<RefCell<Self>> {
self.into()
}
/// Creates a rand instance, pre-seeded with the current time in nanoseconds. /// Creates a rand instance, pre-seeded with the current time in nanoseconds.
/// Needs stdlib timer /// Needs stdlib timer
#[cfg(feature = "std")] #[cfg(feature = "std")]
pub fn preseeded() -> Self { pub fn preseeded() -> Self {
let seed = SystemTime::now() Self::new(current_nanos())
.duration_since(UNIX_EPOCH) }
.unwrap() }
.as_nanos() as u64;
Self::new(seed) /// Extremely quick rand implementation
/// see https://arxiv.org/pdf/2002.11331.pdf
#[derive(Copy, Clone, Debug, Default)]
pub struct RomuTrioRand {
x_state: u64,
y_state: u64,
z_state: u64,
}
impl RomuTrioRand {
pub fn new(seed: u64) -> Self {
let mut rand = Self::default();
rand.set_seed(seed);
rand
}
/// Creates a rand instance, pre-seeded with the current time in nanoseconds.
/// Needs stdlib timer
#[cfg(feature = "std")]
pub fn preseeded() -> Self {
Self::new(current_nanos())
}
}
impl Rand for RomuTrioRand {
fn set_seed(&mut self, seed: u64) {
self.x_state = seed ^ 0x12345;
self.y_state = seed ^ 0x6789A;
self.z_state = seed ^ 0xBCDEF;
}
#[inline]
fn next(&mut self) -> u64 {
let xp = self.x_state;
let yp = self.y_state;
let zp = self.z_state;
self.x_state = 15241094284759029579u64.wrapping_mul(zp);
self.y_state = yp.wrapping_sub(xp).rotate_left(12);
self.z_state = zp.wrapping_sub(yp).rotate_left(44);
xp
}
}
/// see https://arxiv.org/pdf/2002.11331.pdf
#[derive(Copy, Clone, Debug, Default)]
pub struct RomuDuoJrRand {
x_state: u64,
y_state: u64,
}
impl RomuDuoJrRand {
pub fn new(seed: u64) -> Self {
let mut rand = Self::default();
rand.set_seed(seed);
rand
}
}
impl Rand for RomuDuoJrRand {
fn set_seed(&mut self, seed: u64) {
self.x_state = seed ^ 0x12345;
self.y_state = seed ^ 0x6789A;
}
#[inline]
fn next(&mut self) -> u64 {
let xp = self.x_state;
self.x_state = 15241094284759029579u64.wrapping_mul(self.y_state);
self.y_state = self.y_state.wrapping_sub(xp).rotate_left(27);
xp
} }
} }