wolffd@0: function smallpot = marginalize_pot(bigpot, keepdom, sumoverdom, nodesizes)
wolffd@0: % MARGINALIZE_POT Marginalize a mpot onto a smaller domain.
wolffd@0: % smallpot = marginalize_pot(bigpot, keep)
wolffd@0: 
wolffd@0: keepsize = sum(nodesizes(keepdom));
wolffd@0: [A1, A2, B1, B2, C11, C12, C21, C22] = partition_matrix_vec_3(bigpot.A, bigpot.B, bigpot.C, keepdom, sumoverdom, nodesizes);
wolffd@0: smallpot = scgcpot(keepsize, bigpot.ctailsize, bigpot.p, A1, B1, C11);