From 533fd64f3af0431f669e7957023fa864d9a6a337 Mon Sep 17 00:00:00 2001 From: Leo Singer Date: Tue, 11 Sep 2018 22:38:52 -0400 Subject: [PATCH] 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. --- ligo/skymap/bayestar/ez_emcee.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/ligo/skymap/bayestar/ez_emcee.py b/ligo/skymap/bayestar/ez_emcee.py index b91327b..664647e 100644 --- a/ligo/skymap/bayestar/ez_emcee.py +++ b/ligo/skymap/bayestar/ez_emcee.py @@ -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 -- GitLab