Questions about Call Random_number in gfortran

[msz59]

@RonShepard

Gfortran uses the xoshiro256** PRNG algorithm which is described here: Xorshift - Wikipedia . It has a very long cycle, which is why the internal state requires 33 integers. It has the remarkable property that it is possible to skip forward an arbitrary number of steps. This allows multiple PRNG streams to be generated in independent threads and/or on separate compute nodes that are guaranteed to not overlap with each other.

Very interesting. Two questions/comments:

1. those 32 integers (I guess your 33 is a typo) are probably counted in 8-bit values, as
   call random_seed(size = n)
   returns n==8 (gfortran 13.2.0 x86_64-pc-linux-gnu)
2. do you know whether (and how) can one initialize such independent, non-overlapping PRNG streams? I guess it should be done with proper seeding, but how?