topic Okay, now I have a better in Intel® oneAPI Math Kernel Library

initializing multiple random number generator streams for parallel code

cygnetmama — Sat, 05 Mar 2016 00:07:33 GMT

I am using the set of 6024 Mersenne Twister pseudorandom number generators to generate multiple, independent sequences of random numbers within an OpenMP parallel region. When I initialized the random number generators, I wrote the streams to file so I could be sure they were different. One of the two numbers is unique for each stream, but the second number is identical for all streams. For example, for 4 threads, I get

thread           3 stream -2039061504          49
thread           1 stream -2039060224          49
thread           2 stream -2039060864          49
thread           0 stream -2039057664          49

Is this okay or does it mean I have done something wrong? A stripped down version of the relevant code appears below.

Thank you!

      PROGRAM Rng19_2
      USE ModApril2013
      USE Tmod19
      USE omp_lib
      USE MKL_VSL_TYPE
      USE MKL_VSL

      IMPLICIT NONE
      TYPE(VSL_STREAM_STATE)::stream

      MySd=23
      nruns=200
      CALL omp_set_num_threads(nThrds)

      !$OMP PARALLEL DEFAULT(NONE)
     &PRIVATE (irun,URVs,id,errcode,stream,ThrdSd)
     &SHARED (nruns,MySd)

      id = omp_get_thread_num()
      ThrdSd = MySd + id
      errcode=vslnewstream(stream, VSL_BRNG_MT2203+id, ThrdSd)
      write(79,*) 'thread', id, 'stream', stream

      !$OMP DO SCHEDULE(STATIC,nruns/nThrds)
      NewPopulationLoop: DO irun=1,nruns

        errcode=VDRngUniform(VSL_RNG_METHOD_UNIFORM_STD,stream,nURVs,
     & URVs,0.0_DP,1.0_DP)

      END DO NewPopulationLoop
      !$OMP END DO

      errcode=vsldeletestream( stream )

!$OMP END PARALLEL

END PROGRAM Rng19_2

The numbers you printed if

Ilya_B_Intel — Thu, 10 Mar 2016 11:13:35 GMT

The numbers you printed if combined togather identify the stream, and it is ok if some of them are identical but not all at once.

We run your code (with some dependencies removed) and that should work as intended.

In our run we see the following numbers with streams:

thread id=           3 stream= -1140778368       32645
thread id=           2 stream= -939521408       32645
thread id=           0 stream=    47153152           0
thread id=           1 stream= -1006630272       32645

The numbers generated by those streams in threads are different. Here are first 2 numbers from different threads:

URVs= 0.327342507662252       0.882119875401258
URVs= 0.256897222949192       0.952226318651810
URVs= 0.515935511561111       0.137881062459201
URVs= 0.584926510695368       0.183347462909296

Thank you so much for helping

cygnetmama — Thu, 10 Mar 2016 19:16:02 GMT

Thank you so much for helping me with this!

I have one more question, please. Did you use a different seed or did you use MySd=23 . I would have thought with the same seed that you would have been able to replicate my results.

Thanks again.

Hello cygnetmama,

Eugeny_G_Intel — Fri, 11 Mar 2016 10:21:14 GMT

Hello cygnetmama,

Here is a program we used (it is your test-case simplified so that we can run it). The seed is ThrdSd = MySd + id. And we generate 2 random values. There is no sense to print stream. It is just internal library id, different for each stream. There is sense to print what stream generates, i.e. URVs in our case, and it is reproduceable on same processor, if the stream initialized with same method, seed, etc.

include 'mkl_vsl.f90'

      PROGRAM Rng19_2
      USE omp_lib
      USE MKL_VSL_TYPE
      USE MKL_VSL

      !IMPLICIT NONE
      TYPE(VSL_STREAM_STATE)::stream
      REAL(8)               ::URVs(2)
      INTEGER               ::ThrdSd
      REAL(8)               ::A0,A1

      A0=0.0
      A1=1.0
      nURVs=2
      MySd=23
      nruns=2
      CALL omp_set_num_threads(4) !nThrds

      !$OMP PARALLEL DEFAULT(NONE)
     &PRIVATE (irun,URVs,id,errcode,stream,ThrdSd)
     &SHARED (nruns,nURVs,MySd,A0,A1)

      id = omp_get_thread_num()
      ThrdSd = MySd + id
      errcode=vslnewstream(stream, VSL_BRNG_MT2203+id, ThrdSd)
      write(*,*) 'thread id=', id, 'stream=', stream

        errcode=VDRngUniform(VSL_RNG_METHOD_UNIFORM_STD,stream,
     & nURVs,URVs,A0,A1)
      write(*,*) 'thread id =',id, 'URVs=', URVs

errcode=vsldeletestream( stream )

!$OMP END PARALLEL

END PROGRAM Rng19_2

Currently output is (we run on Haswell):

thread id=           0 stream=    19902464           0
thread id=           3 stream= -1946154368       32553
thread id=           1 stream= -1811936640       32553
thread id=           2 stream= -1744827776       32553
thread id =           1 URVs= 0.515935511561111       0.137881062459201
thread id =           2 URVs= 0.584926510695368       0.183347462909296
thread id =           0 URVs= 0.256897222949192       0.952226318651810
thread id =           3 URVs= 0.327342507662252       0.882119875401258

Okay, now I have a better

cygnetmama — Fri, 11 Mar 2016 16:54:02 GMT

Okay, now I have a better understanding. Thank you very much for helping me!