wolffd@0: function smallpot = marginalize_pot(bigpot, onto, maximize)
wolffd@0: % MARGINALIZE_POT Marginalize a upot onto a smaller domain.
wolffd@0: % smallpot = marginalize_pot(bigpot, onto, maximize)
wolffd@0: %
wolffd@0: % The maximize argument is ignored
wolffd@0: 
wolffd@0: numer = marg_table(bigpot.p .* bigpot.u, bigpot.domain, bigpot.sizes, onto);
wolffd@0: denom = marg_table(bigpot.p, bigpot.domain, bigpot.sizes, onto);
wolffd@0: 
wolffd@0: p = denom;
wolffd@0: % replace 0s by 1s before dividing. This is valid since demon(i) = 0 => numer(i) = 0
wolffd@0: denom = denom + (denom == 0); 
wolffd@0: u = numer ./ denom;
wolffd@0: 
wolffd@0: ns = zeros(1, max(bigpot.domain));
wolffd@0: ns(bigpot.domain) = bigpot.sizes;
wolffd@0: 
wolffd@0: smallpot = upot(onto, ns(onto), p, u);