Commit 533fd64f authored by Leo Pound Singer's avatar Leo Pound Singer

Use global Numpy random state

I am seeding the global random number generator with
`np.random.seed`. I want MCMC sampling to respect the state of
this random number generator.
parent 971ea9d9
......@@ -100,7 +100,7 @@ def ez_emcee(log_prob_fn, lo, hi, nindep=200,
sampler = Sampler(nwalkers, ndim, log_prob_fn, logp,
ntemps=ntemps, loglargs=args, loglkwargs=kwargs,
logpargs=[lo, hi], **options)
logpargs=[lo, hi], random=np.random, **options)
pos = np.random.uniform(lo, hi, (ntemps, nwalkers, ndim))
# Burn in
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment