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