wolffd@0: function ens = comp_eff_node_sizes(ns, cnodes, ev, domain)
wolffd@0: 
wolffd@0: dnodes = mysetdiff(1:length(ns), cnodes);
wolffd@0: odom = domain(~isemptycell(evidence(domain)));
wolffd@0: cdom = myintersect(cnodes, domain);
wolffd@0: ddom = myintersect(dnodes, domain);
wolffd@0: cobs = myintersect(cdom, odom);
wolffd@0: dobs = myintersect(ddom, odom);
wolffd@0: ens = ns; 
wolffd@0: ens(cobs) = 0;
wolffd@0: ens(dobs) = 1;