Skip to content

Commit

Permalink
Merge remote-tracking branch 'upstream/master'
Browse files Browse the repository at this point in the history
  • Loading branch information
vhaasteren committed Mar 31, 2016
2 parents 5831c50 + 170157f commit eb86802
Showing 1 changed file with 7 additions and 5 deletions.
12 changes: 7 additions & 5 deletions PTMCMCSampler/PTMCMCSampler.py
Original file line number Diff line number Diff line change
Expand Up @@ -60,8 +60,7 @@ class PTSampler(object):

def __init__(self, ndim, logl, logp, cov, groups=None, loglargs=[], loglkwargs={},
logpargs=[], logpkwargs={}, logl_grad=None, logp_grad=None,
comm=MPI.COMM_WORLD,
outDir='./chains', verbose=True, resume=False):
comm=MPI.COMM_WORLD, outDir='./chains', verbose=True, resume=False):

# MPI initialization
self.comm = comm
Expand Down Expand Up @@ -191,13 +190,16 @@ def initialize(self, Niter, ladder=None, Tmin=1, Tmax=None, Tskip=100,
self.addProposalToCycle(nutsjump, NUTSweight)

# add SCAM
self.addProposalToCycle(
self.covarianceJumpProposalSCAM,
self.SCAMweight)
self.addProposalToCycle(self.covarianceJumpProposalSCAM,
self.SCAMweight)

# add AM
self.addProposalToCycle(self.covarianceJumpProposalAM, self.AMweight)

# check length of jump cycle
if len(self.propCycle) == 0:
raise ValueError('No jump proposals specified!')

# randomize cycle
self.randomizeProposalCycle()

Expand Down

0 comments on commit eb86802

Please sign in to comment.