Infinite looping in hierarchy_minimize
Hello. I'm trying to test the latest version of graph_tool (0.23dev) and building the nested sbm with edge weights. I'm launching the following command: minimize_nested_blockmodel_dl( graph, deg_corr=False, b_min = graph.vp.ikind, B_max = 150, B_min = 5, mcmc_equilibrate_args = dict(wait = 100, nbreaks=2, max_niter=100000), mcmc_args = dict(niter=10), state_args = dict( clabel = graph.vp.ikind, recs = [ln2rev], rec_types = ['real-normal'] ), verbose=True, ) And can see that there is infinite looping: ... level 386 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 387 l=387 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=387 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=387 Bisect at B = 2 with S = 34.5778592783179 l=387 Best result: B = 2, S = 34.5778592783179 level 387 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 388 l=388 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=388 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=388 Bisect at B = 2 with S = 34.5778592783179 l=388 Best result: B = 2, S = 34.5778592783179 level 388 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 389 l=389 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=389 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=389 Bisect at B = 2 with S = 34.5778592783179 l=389 Best result: B = 2, S = 34.5778592783179 level 389 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 390 l=390 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=390 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=390 Bisect at B = 2 with S = 34.5778592783179 l=390 Best result: B = 2, S = 34.5778592783179 level 390 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 391 l=391 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=391 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=391 Bisect at B = 2 with S = 34.5778592783179 l=391 Best result: B = 2, S = 34.5778592783179 level 391 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 392 l=392 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=392 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=392 Bisect at B = 2 with S = 34.5778592783179 l=392 Best result: B = 2, S = 34.5778592783179 level 392 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 393 l=393 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=393 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=393 Bisect at B = 2 with S = 34.5778592783179 l=393 Best result: B = 2, S = 34.5778592783179 level 393 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 394 l=394 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=394 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=394 Bisect at B = 2 with S = 34.5778592783179 l=394 Best result: B = 2, S = 34.5778592783179 level 394 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 395 l=395 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=395 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=395 Bisect at B = 2 with S = 34.5778592783179 l=395 Best result: B = 2, S = 34.5778592783179 level 395 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 396 l=396 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=396 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=396 Bisect at B = 2 with S = 34.5778592783179 l=396 Best result: B = 2, S = 34.5778592783179 level 396 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 397 l=397 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=397 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=397 Bisect at B = 2 with S = 34.5778592783179 l=397 Best result: B = 2, S = 34.5778592783179 level 397 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 398 l=398 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=398 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=398 Bisect at B = 2 with S = 34.5778592783179 l=398 Best result: B = 2, S = 34.5778592783179 level 398 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 399 l=399 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=399 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=399 Bisect at B = 2 with S = 34.5778592783179 l=399 Best result: B = 2, S = 34.5778592783179 level 399 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 400 l=400 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=400 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=400 Bisect at B = 2 with S = 34.5778592783179 l=400 Best result: B = 2, S = 34.5778592783179 What is wrong with the launch of minimization ? Thank you Valeriy
On 06.09.2017 20:21, Valery Topinsky wrote:
Hello.
I'm trying to test the latest version of graph_tool (0.23dev) and building the nested sbm with edge weights.
I'm launching the following command:
minimize_nested_blockmodel_dl( graph, deg_corr=False, b_min = graph.vp.ikind, B_max = 150, B_min = 5, mcmc_equilibrate_args = dict(wait = 100, nbreaks=2, max_niter=100000), mcmc_args = dict(niter=10), state_args = dict( clabel = graph.vp.ikind, recs = [ln2rev], rec_types = ['real-normal'] ), verbose=True, )
And can see that there is infinite looping: ... level 386 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 387 l=387 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=387 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=387 Bisect at B = 2 with S = 34.5778592783179 l=387 Best result: B = 2, S = 34.5778592783179 level 387 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 388 l=388 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=388 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=388 Bisect at B = 2 with S = 34.5778592783179 l=388 Best result: B = 2, S = 34.5778592783179 level 388 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 389 l=389 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=389 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=389 Bisect at B = 2 with S = 34.5778592783179 l=389 Best result: B = 2, S = 34.5778592783179 level 389 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 390 l=390 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=390 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=390 Bisect at B = 2 with S = 34.5778592783179 l=390 Best result: B = 2, S = 34.5778592783179 level 390 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 391 l=391 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=391 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=391 Bisect at B = 2 with S = 34.5778592783179 l=391 Best result: B = 2, S = 34.5778592783179 level 391 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 392 l=392 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=392 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=392 Bisect at B = 2 with S = 34.5778592783179 l=392 Best result: B = 2, S = 34.5778592783179 level 392 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 393 l=393 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=393 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=393 Bisect at B = 2 with S = 34.5778592783179 l=393 Best result: B = 2, S = 34.5778592783179 level 393 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 394 l=394 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=394 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=394 Bisect at B = 2 with S = 34.5778592783179 l=394 Best result: B = 2, S = 34.5778592783179 level 394 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 395 l=395 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=395 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=395 Bisect at B = 2 with S = 34.5778592783179 l=395 Best result: B = 2, S = 34.5778592783179 level 395 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 396 l=396 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=396 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=396 Bisect at B = 2 with S = 34.5778592783179 l=396 Best result: B = 2, S = 34.5778592783179 level 396 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 397 l=397 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=397 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=397 Bisect at B = 2 with S = 34.5778592783179 l=397 Best result: B = 2, S = 34.5778592783179 level 397 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 398 l=398 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=398 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=398 Bisect at B = 2 with S = 34.5778592783179 l=398 Best result: B = 2, S = 34.5778592783179 level 398 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 399 l=399 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=399 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=399 Bisect at B = 2 with S = 34.5778592783179 l=399 Best result: B = 2, S = 34.5778592783179 level 399 : replaced (2, 1) -> (2, 2) , dS: -15.0337981072 400 l=400 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=400 Current bracket: (2, 2, 2) (34.577859278317902, 34.577859278317902, 34.577859278317902) l=400 Bisect at B = 2 with S = 34.5778592783179 l=400 Best result: B = 2, S = 34.5778592783179
What is wrong with the launch of minimization ?
Please provide a minimal and self-contained program that shows the problem (including the network used) so I can investigate. Best, Tiago -- Tiago de Paula Peixoto <tiago@skewed.de>
participants (2)
-
Tiago de Paula Peixoto -
Valery Topinsky