Daniel@0: function fmarginal = add_ev_to_dmarginal(fmarginal, evidence, ns)
Daniel@0: % ADD_EV_TO_DMARGINAL 'pump up' observed nodes back to their original size.
Daniel@0: % fmarginal = add_ev_to_dmarginal(fmarginal, evidence, ns)
Daniel@0: %
Daniel@0: % We introduce 0s into the array in positions which are incompatible with the evidence.
Daniel@0: 
Daniel@0: dom = fmarginal.domain;
Daniel@0: odom = dom(~isemptycell(evidence(dom)));
Daniel@0: vals = cat(1, evidence{odom});
Daniel@0: index = mk_multi_index(length(dom), find_equiv_posns(odom, dom), vals);
Daniel@0: T = 0*myones(ns(dom));
Daniel@0: ens = ns(:)';
Daniel@0: ens(odom) = 1;
Daniel@0: T(index{:}) = myreshape(fmarginal.T, ens(dom));
Daniel@0: fmarginal.T = T;