summate.c: sum_with_equalities: extract out sum_over_polytope_0D