function marginal = marginal_nodes(engine, b, nodes, t, add_ev)
% MARGINAL_NODES Compute the marginal on the specified nodes (hmm_2TBN)
% marginal = marginal_nodes(engine, b, nodes, t, add_ev)
%
% nodes must be a singleton set 

assert(length(nodes)==1)
ss = engine.slice_size;

i = nodes(1);
bigT = b.gamma;
dom = i + (t-1)*ss;

%id = engine.marg_singleton_ndx_id(i);
%global SD_NDX
%ndx = SD_NDX{id};
%marginal.T = marg_table_ndxSD(bigT, engine.maximize, ndx);

ns = engine.eff_node_sizes(:);
bigdom = 1:ss;
marginal.T = marg_table(bigT, bigdom + (t-1)*ss, ns(bigdom), dom, engine.maximize);

marginal.domain = dom;
assert(~add_ev);
%if add_ev
%  marginal = add_ev_to_dmarginal(marginal, engine.evidence, engine.node_sizes);
%end