Questions about Call Random_number in gfortran

[Beliavsky]

@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.

If it is possible to skip ahead with gfortran’s random_number(), how would you do so? I know you can skip n random variates with

do i=1,n
   call random_number(x)
end do