Infinite looping in hierarchy_minimize

graph-tool
1

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

attachment.html (7.15 KB)

2

Please provide a minimal and self-contained program that shows the problem
(including the network used) so I can investigate.

Best,
Tiago