[graph-tool] Infinite looping in hierarchy_minimize

6 Sep 2017

      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

[graph-tool] Infinite looping in hierarchy_minimize

Valery Topinsky