Also, I think we should first decide how to specify default rng. This would everything much simpler. To test whether ``use_default_rng`` in a tight inner loop is performance killer and not a very elegant solution.