Skip to content
Snippets Groups Projects
Commit f8ef7213 authored by Gregory Ashton's avatar Gregory Ashton
Browse files

Fix missing default kwarg in cpnest

parent ffc85a5d
No related branches found
No related tags found
No related merge requests found
...@@ -40,7 +40,7 @@ class Cpnest(NestedSampler): ...@@ -40,7 +40,7 @@ class Cpnest(NestedSampler):
""" """
default_kwargs = dict(verbose=1, nthreads=1, nlive=500, maxmcmc=1000, default_kwargs = dict(verbose=1, nthreads=1, nlive=500, maxmcmc=1000,
seed=None, poolsize=100, nhamiltonian=0, resume=True, seed=None, poolsize=100, nhamiltonian=0, resume=True,
output=None, proposals=None) output=None, proposals=None, n_periodic_checkpoint=None)
def _translate_kwargs(self, kwargs): def _translate_kwargs(self, kwargs):
if 'nlive' not in kwargs: if 'nlive' not in kwargs:
...@@ -84,17 +84,22 @@ class Cpnest(NestedSampler): ...@@ -84,17 +84,22 @@ class Cpnest(NestedSampler):
self._resolve_proposal_functions() self._resolve_proposal_functions()
model = Model(self.search_parameter_keys, self.priors) model = Model(self.search_parameter_keys, self.priors)
try: out = None
out = CPNest(model, **self.kwargs) remove_kwargs = ["proposals", "n_periodic_checkpoint"]
except TypeError as e: while out is None:
if 'proposals' in self.kwargs.keys(): try:
logger.warning('YOU ARE TRYING TO USE PROPOSALS IN A VERSION OF CPNEST THAT DOES'
'NOT ACCEPT CUSTOM PROPOSALS. SAMPLING WILL COMMENCE WITH THE DEFAULT'
'PROPOSALS.')
del self.kwargs['proposals']
out = CPNest(model, **self.kwargs) out = CPNest(model, **self.kwargs)
else: except TypeError as e:
raise TypeError(e) if len(remove_kwargs) > 0:
kwarg = remove_kwargs.pop(0)
else:
raise TypeError("Unable to initialise cpnest sampler")
logger.info(
"CPNest init. failed with error {}, please update"
.format(e))
logger.info(
"Attempting to rerun with kwarg {} removed".format(kwarg))
self.kwargs.pop(kwarg)
out.run() out.run()
if self.plot: if self.plot:
......
...@@ -112,13 +112,15 @@ class TestCPNest(unittest.TestCase): ...@@ -112,13 +112,15 @@ class TestCPNest(unittest.TestCase):
def test_default_kwargs(self): def test_default_kwargs(self):
expected = dict(verbose=1, nthreads=1, nlive=500, maxmcmc=1000, expected = dict(verbose=1, nthreads=1, nlive=500, maxmcmc=1000,
seed=None, poolsize=100, nhamiltonian=0, resume=True, seed=None, poolsize=100, nhamiltonian=0, resume=True,
output='outdir/cpnest_label/', proposals=None) output='outdir/cpnest_label/', proposals=None,
n_periodic_checkpoint=None)
self.assertDictEqual(expected, self.sampler.kwargs) self.assertDictEqual(expected, self.sampler.kwargs)
def test_translate_kwargs(self): def test_translate_kwargs(self):
expected = dict(verbose=1, nthreads=1, nlive=250, maxmcmc=1000, expected = dict(verbose=1, nthreads=1, nlive=250, maxmcmc=1000,
seed=None, poolsize=100, nhamiltonian=0, resume=True, seed=None, poolsize=100, nhamiltonian=0, resume=True,
output='outdir/cpnest_label/', proposals=None) output='outdir/cpnest_label/', proposals=None,
n_periodic_checkpoint=None)
for equiv in bilby.core.sampler.base_sampler.NestedSampler.npoints_equiv_kwargs: for equiv in bilby.core.sampler.base_sampler.NestedSampler.npoints_equiv_kwargs:
new_kwargs = self.sampler.kwargs.copy() new_kwargs = self.sampler.kwargs.copy()
del new_kwargs['nlive'] del new_kwargs['nlive']
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment