wolffd@0: function marg2 = shrink_obs_dims_in_gaussian(marg1, dom, evidence, ns)
wolffd@0: % SHRINK_OBS_DIMS_IN_GAUSSIAN Remove observed dimensions from mu/Sigma
wolffd@0: % function marg2 = shrink_obs_dims_in_gaussian(marg1, dom, evidence, ns)
wolffd@0: 
wolffd@0: % This is used by loopy
wolffd@0: 
wolffd@0: hdom = dom(isemptycell(evidence(dom)));
wolffd@0: ndx = find_equiv_posns(hdom, dom);
wolffd@0: b = block(ndx, ns(dom));
wolffd@0: marg2.mu = marg1.mu(b);
wolffd@0: marg2.Sigma = marg1.Sigma(b,b);
wolffd@0: marg2.domain = marg1.domain;