Function sim::fork_rng

pub fn fork_rng(base_rng: &mut XorShiftRng) -> XorShiftRng

Need to explain this trick -- basically keeps consistency between two different simulations when each one might make slightly different sequences of calls to the RNG.