vmagnin / exploring_coarrays Goto Github PK

Hi,
Thanks for these very useful examples. Just a remark on the seeding of the PRNG. instead of:

 ! Each image will have its own RNG seed:
 call rng%seed([ -1337_i8, 9812374_i8 ] + 10*this_image())

I would suggest something like:

call rng%seed([ -1337_i8, 9812374_i8 ])
if (this_image() /= 1) then
  do i = 2, this_image()
    call rng%jump()
  end do
end if

The PRNG sequences on the different images are then guaranteed not to overlap.

Best,

Christophe

I suggest adding the topics coarray-fortran and openmp in the About section.