Commit e81cdd98 authored by Tiago Peixoto's avatar Tiago Peixoto

PPBlockState(): add compatibility functions

This fixes issue #668
parent d546bb2b
Pipeline #688 passed with stage
in 100 minutes and 16 seconds
......@@ -130,6 +130,10 @@ class PPBlockState(object):
r"Returns the total number of blocks."
return len(np.unique(self.b.fa))
def get_nonempty_B(self):
r"Alias to :meth:`~PPBlockState.get_B`."
return self.get_B()
def get_Be(self):
r"""Returns the effective number of blocks, defined as :math:`e^{H}`, with
:math:`H=-\sum_r\frac{n_r}{N}\ln \frac{n_r}{N}`, where :math:`n_r` is
......@@ -140,6 +144,18 @@ class PPBlockState(object):
w /= w.sum()
return numpy.exp(-(w*log(w)).sum())
def virtual_vertex_move(self, v, s, **kwargs):
r"""Computes the entropy difference if vertex ``v`` is moved to block ``s``. The
remaining parameters are the same as in
return self._state.virtual_move(int(v), self.b[v], s,
def move_vertex(self, v, s):
r"""Move vertex ``v`` to block ``s``."""
self._state.move_vertex(int(v), int(s))
def entropy(self, uniform=False, degree_dl_kind="distributed", **kwargs):
r"""Return the model entropy (negative log-likelihood).
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