let bernouilli_parallel prng length (d1,d2) =
  assert (0. <= d1 && d1 <= 1.);
  assert (0. <= d2 && d2 <= 1.);
  Array.init length (fun _ -> Random.State.float prng 1. < d1, Random.State.float prng 1. < d2)