add_niels_to_pt
add_niels_to_pt(p, pn->n, before_double);
add_niels_to_pt(out, ni, j == n - 1 && i != 1);
add_niels_to_pt(combo,