wolffd@0: function smallpot = marginalize_pot(bigpot, onto, maximize)
wolffd@0: % MARGINALIZE_POT Marginalize a dpot onto a smaller domain.
wolffd@0: % smallpot = marginalize_pot(bigpot, onto, maximize)
wolffd@0: %
wolffd@0: % 'onto' must be in ascending order.
wolffd@0: 
wolffd@0: if nargin < 3, maximize = 0; end
wolffd@0: 
wolffd@0: ns = zeros(1, max(bigpot.domain));
wolffd@0: ns(bigpot.domain) = bigpot.sizes;
wolffd@0: %assert(isequal(bigpot.sizes, mysize(bigpot.T))); % may fail if there are trailing dimensions of size 1
wolffd@0: if issparse(bigpot.T)
wolffd@0:    smallT = marg_sparse_table(bigpot.T, bigpot.domain, bigpot.sizes, onto, maximize);
wolffd@0: else 
wolffd@0:    smallT = marg_table(bigpot.T, bigpot.domain, bigpot.sizes, onto, maximize);
wolffd@0: end
wolffd@0: smallpot = dpot(onto, ns(onto), smallT);