function mdag = enforce_clique( G_m, bvec )
% function mdag = enforce_clique( G_m, bvec )
% takes G_m and forces variables in bvec to be cliques

mdag = G_m;
% connect all of bvec edges
for u=bvec,
  for v=bvec,
    if u ~= v,
      mdag(u,v)=1;
      %disp(sprintf('connected u(%d),v(%d)',u,v));
    end;
  end;
end;
% checking BEL_t-1 created
for u=bvec,
  for v=bvec,
    if u~= v,
      val = mdag(u,v);
      %disp(sprintf('connected u(%d),v(%d)? %d',u,v,val));
    end;
  end;
end;