wolffd@0: function marginal = marginal_nodes(engine, nodes, t)
wolffd@0: % MARGINAL_NODES Compute the marginal on the specified query nodes (bk_ff_hmm)
wolffd@0: % marginal = marginal_nodes(engine, i, t)
wolffd@0: 
wolffd@0: assert(length(nodes)==1);
wolffd@0: i = nodes(end);
wolffd@0: %assert(myismember(i, engine.hnodes));
wolffd@0: marginal = engine.marginals{i,t};
wolffd@0: bnet = bnet_from_engine(engine);
wolffd@0: ss = length(bnet.intra);
wolffd@0: marginal.domain = i + (t-1)*ss;