Problem with random generator

WschW · April 10, 2020, 5:19am

Each time you call DPP the DeterminantalPointProcesses it creates contains a MersenneTwister created with the same seed.

github.com

alshedivat/DeterminantalPointProcesses.jl/blob/80ffff8c9e33b1cfd43fcfe7379f8ffd80ba7240/src/types.jl#L12


      
          abstract type PointProcess end;
          
          

          
# specific types
          mutable struct DeterminantalPointProcess <: PointProcess
              L::Symmetric
              Lfact::Eigen
              size::Int
              rng::AbstractRNG
          
          
    function DeterminantalPointProcess(L::Symmetric; seed::Int=42)
                  Lfact = eigen(L)
                  new(L, Lfact, length(Lfact.values), MersenneTwister(seed))
              end
          
          
    function DeterminantalPointProcess(Lfact::Eigen; seed::Int=42)
                  L = Symmetric((Lfact.vectors .* Lfact.values') * Lfact.vectors')
                  new(L, Lfact, length(Lfact.values), MersenneTwister(seed))
              end
          end

Since the rand method provided by the package uses that MersenneTwister you will always get the same sequence if you create the DeterminantalPointProcesses without setting the seed to a different number each time you create one.

Topic		Replies	Views
Randomly generate a sequence of numbers one by one using the rand function in Distributions module General Usage	17	2920	April 12, 2019
Comparing random sequences between Julia, Numpy, MATLAB General Usage statistics	3	1453	December 2, 2017
`rand()` vs. `rand.()`? General Usage question , proposal	13	2520	June 18, 2017
Parallel random number generation Statistics	11	2493	February 10, 2017
Identical code produces different results. sometimes General Usage question	12	2681	May 12, 2017

Problem with random generator

Related topics