If you want to fork this repository to propose merge requests, please send an email to tiago@skewed.de, and your project limit will be raised.

Commit f5a3e59a authored by Tiago Peixoto's avatar Tiago Peixoto

Fix MixedMeasuredBlockState.mcmc_sweep()

parent d29ea506
......@@ -585,9 +585,12 @@ class MixedMeasuredBlockState(UncertainBaseState):
b_dh = [max(0, hs[j] - hstep), hs[j] + hstep]
pb = 1./min(1, hs[j])
Sb = self._state.entropy()
latent_edges = kwargs.get("entropy_args", {}).get("latent_edges", True)
density = False
Sb = self._state.entropy(latent_edges, density)
Sa = self._state.entropy()
Sa = self._state.entropy(latent_edges, density)
if Sa < Sb or numpy.random.random() < exp(-(Sa-Sb) + log(pb) - log(pf)):
return (Sa-Sb, 1, 1)
Markdown is supported
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment