Commit 2d315d45 authored by Tiago Peixoto's avatar Tiago Peixoto
Browse files

Forbid creation of empty blocks in mcmc_sweep()

parent 1441551e
......@@ -1197,6 +1197,10 @@ void move_sweep(EMprop mrs, Vprop mrp, Vprop mrm, Vprop wr, Vprop b,
vertex_t v = vertex(vlist[i], g);
vertex_t r = b[v];
// blocks can't become empty
if (wr[r] == 1)
continue;
if (!sequential)
moves[i] = r;
......
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