multiflip_mcmc_sweep() with pclabel reaches `nan` entropy at top state
A minimal example:
import graph_tool.all as gt
import numpy as np
g = gt.collection.data["celegansneural"]
N = g.num_vertices()
pclabel = np.array([int(i > N / 2) for i in range(N)])
nbstate = gt.minimize_nested_blockmodel_dl(g, state_args={"pclabel": pclabel})
while not np.isnan(nbstate.entropy()):
nbstate.multiflip_mcmc_sweep(beta=np.inf, niter=10)
level_entropies = {i: nbstate.level_entropy(i) for i in range(len(nbstate.levels))}
print(level_entropies)
Edited by Ale Abdo