def compute_monomials(R, nverts, dresverts):
"""Eq. A.8 from arXiv:1411.1418"""
X = np.array(R.gens())
Delta = Polyhedron(vertices=nverts)
npoints = np.array(Delta.integral_points())
P_expon = (npoints @ dresverts.T) + 1
P_monoms = np.power(X, P_expon).prod(axis=1)
return P_monoms
