REVIEW 4 minor 14 references
Simple Formula for Integration of Polynomials on a Simplex
T0 review · 0 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Integrating a degree-t polynomial over a simplex reduces to t point evaluations of modified homogeneous forms.
desk verdict Correct, clean, and honest about its own scope; the formula is a repackaging of a classical identity, but a useful one. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central mechanism is a Laplace-transform identity: for a positively homogeneous function $f$ of degree $t > -(n+1)$, $$\int_{\Delta_z} f(x)\,dx = \frac{1}{\Gamma(1+n+t)}\int_{\mathbb{R}_+^n} f(x)\exp(-z^T x)\,dx,$$ obtained by writing the simplex integral as a one-dimensional integral over levels $\{z^T x \le y\}$ and applying Fubini. This reduces monomial integrals to closed-form Gamma-factor expressions, which for polynomials become factorial products. The Bombieri polynomial $\hat{f}(x) = \sum_\alpha \alpha! f_\alpha x^\alpha$ is the object that absorbs these factorial weights, and the evaluation points $\xi_j$ lie on the line segment from the origin to $e/n$ inside the simplex.
What would settle it
Take the standard 2-simplex and the monomial $f(x)=x_1 x_2^3$, whose exact integral is $1/120$; evaluating the paper's formula gives $\frac{1}{2}\hat{f}_4(e/360^{1/4}) = \frac{1}{2}\cdot\frac{6}{360} = \frac{1}{120}$, so any deviation from this value in a direct numerical quadrature check would indicate a misapplication or error.
Extended reading notes
Core claim
The paper proves that for the canonical simplex $\Delta = \{x \ge 0, e^T x \le 1\} \subset \mathbb{R}^n$ and any polynomial $f$ of total degree $t$, written as $f = \sum_{j=0}^t f_j$ with $f_j$ homogeneous of degree $j$, the integral satisfies $$\int_\$\Delta$ f\,dx = \frac{1}{n!}\left(\hat{f}_0 + \sum_{j=1}^t \hat{f}_j(\xi_j)\right), \qquad \xi_j = \frac{e}{((n+1)\cdots(n+j))^{1/j}},$$ where $\hat{f}_j$ is the Bombieri (factorial-weighted) form of $f_j$, i.e. each monomial coefficient $f_\alpha$ is replaced by $\alpha_1!\cdots\alpha_n! f_\alpha$. For a single homogeneous polynomial of degree $t$ the formula simplifies to $\int_\Delta f\,dx = \hat{f}(e)/(n+t)!$, and an affine change of variables transfers the formula to an arbitrary simplex. The paper also extends the identity to positively homogeneous functions of the form $\sum_\alpha f_\alpha x^\alpha$ with all $\alpha_i > -1$, using Gamma factors in place of factorials.
Load-bearing premise
The proof leans on interchanging the simplex integral with an exponential (Laplace) integral, which is valid only when the integrand is absolutely integrable and its homogeneity degree exceeds $-(n+1)$; ordinary polynomials always satisfy this, but the extended formula for non-polynomial functions breaks down when any exponent is at or below $-1$.
Editorial extensions
If this is right
- Integrating any degree-$t$ polynomial over any full-dimensional simplex requires exactly $t$ point evaluations of the modified homogeneous parts, independent of the dimension in terms of the number of points.
- For a power of a linear form, $\int_\Delta (\ell\cdot x)^t dx = t!\, E_t(\ell)/(n+t)!$, where $E_t$ is the complete homogeneous symmetric polynomial of degree $t$; this gives an immediate closed form for such integrals.
- In finite element and extended finite element methods, elementwise polynomial integrals can be computed exactly by this formula, potentially replacing cubature rules that sample the original polynomial at many points.
- The evaluation points $\xi_j$ are easy to generate recursively, since $\xi_j = e/\theta_j^{1/j}$ with $\theta_j = (n+1)\cdots(n+j)$.
- The extension to positively homogeneous functions with real exponents $\alpha_i > -1$ supplies exact integrals for a class of non-polynomial integrands that still arise in applications.
Reading between the lines
- The author stops short of saying so, but the formula implies that classical cubature rules with many points are unnecessary for exact polynomial integration over a simplex: the entire information of the integral is carried by the $t$ modified forms at $t$ aligned points.
- The Laplace-transform identity may be useful beyond the stated class: any function whose integral over the simplex can be computed through the same homogeneity argument yields a similar one-point evaluation, though the absolute-integrability condition $\alpha_i > -1$ is essential.
- The Waring-decomposition subsection suggests a practical acceleration: if a form has low Waring rank, its simplex integral can be computed from the values $E_k(\ell_i)$ for the few linear forms in the decomposition, rather than from the full monomial expansion; this could be tested against existing implementations.
- Because the point set $\{\xi_j\}$ is one-dimensional, the formula might extend to finite-element mass and stiffness matrices by differentiating with respect to parameters inside the polynomial, a route the paper leaves unexplored.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that integrating a polynomial of total degree t on an arbitrary full-dimensional simplex reduces, after an affine change of variables, to t point evaluations: on the canonical simplex ∆ = {x ≥ 0, e^T x ≤ 1}, the integral is 1/n! times the sum over j = 0,...,t of the Bombieri (factorial-weighted) forms \hat f_j evaluated at ξ_j = e/((n+1)...(n+j))^{1/j}. The proof uses homogeneity, a Laplace-transform identity (Theorem 2.1), and closed-form monomial exponential integrals; the result is then transferred to arbitrary simplices. A similar formula is stated for a class of positively homogeneous functions with exponents α_i > -1, and a two-dimensional example is checked.
Significance. If correct, the formula is a clean and practically useful reformulation of classical simplex integration: it avoids cubature weights and vertex summations, requiring only evaluations of related homogeneous polynomials at t aligned points. The derivation is self-contained, has no fitted parameters, and is verified on a concrete example. The author appropriately acknowledges that the identity can be obtained from earlier results; the contribution is therefore primarily an elegant and explicit simplification rather than a fundamentally new method. The non-polynomial extension is a modest additional result. I find the central claim sound.
minor comments (4)
- [Section 2.2, Eq. (2.4)] The factor written as "1/z_e" is ambiguous and not defined in the notation; it should be written as z^{-e} = 1/(z_1...z_n) (or explicitly as ∏_i 1/z_i) to avoid confusion with the vector e.
- [Section 2.3, Eq. (2.11)] The change of variables y = A(x-a) gives dx = dy/|det A|, so Eq. (2.11) should use 1/|det A|, or the text should state that the vertex ordering is chosen so that det A > 0.
- [Section 2.2, Proposition 2.2 proof] The notation z^{-e} appears without a definition; since e already denotes the all-ones vector, the proof should define z^{-e} explicitly and use it consistently with the notation in Eq. (2.4).
- [Corollary 2.1] In the display after Eq. (2.8), the summation is written as f = ∑_{j0} f_j; it should be ∑_{j=0}^t f_j.
Circularity Check
No circularity: Theorem 1.1 is derived self-containedly from the classical monomial integral and standard Laplace/Fubini identities.
full rationale
The paper's central formula (1.2), restated as (2.8), is not an input in disguise. The proof decomposes f into homogeneous parts, applies Theorem 2.1 to each part, and evaluates the monomial integral ∫_{R_+^n} y^α exp(-z^T y) dy = α! / z^{α+e} in closed form. Summing these evaluations gives exactly the Bombieri-type weighted form, so the factorial weights are a bookkeeping device rather than fitted parameters. The key identity (2.3) is derived in the paper from the homogeneity of h(y), Fubini-Tonelli, and the Laplace transform of y^{n+t}, all standard external facts; no result of the author is used to force the conclusion. The cited earlier work [9,10,11] is mentioned as background for the Laplace-transform technique, but the proof does not depend on those papers. The claimed novelty is qualified by the author's own admission that the formula could be obtained from previous results, which is an honesty about provenance rather than circular reasoning. The formula is benchmarked against the classical Dirichlet/monomial integral ∫_Δ x^α dx = α!/(n+|α|)!, and the worked Example 2.1 checks out numerically. Two small issues—the missing absolute value of det(A) in (2.11) for orientation-reversing maps and the reliance on standard convergence conditions for the non-polynomial extension—are correctness or presentation concerns, not circularity. No fitted input is renamed as a prediction, and no self-citation chain bears the derivation's weight. The overall circularity score is therefore 0.
Assumptions & free parameters
assumptions (4)
- standard math Fubini-Tonelli theorem applies to the iterated integrals in the proof of Theorem 2.1.
- standard math Gamma integral identity: integral from 0 to infinity of y^k exp(-z y) dy = k!/z^(k+1) for integers k.
- standard math Any full-dimensional simplex is affinely equivalent to the canonical simplex Delta via y = A(x-a).
- standard math Every polynomial decomposes uniquely as a finite sum of homogeneous forms.
Cite this review
Pith. "Pith review of Simple Formula for Integration of Polynomials on a Simplex." pith.science (2026). https://pith.science/paper/Z4CGVVGA
@misc{pith2026190806736,
author = {Pith},
title = {Pith review of: Simple Formula for Integration of Polynomials on a Simplex},
year = {2026},
howpublished = {\url{https://pith.science/paper/Z4CGVVGA}},
note = {Machine review of arXiv:1908.06736}
}
abstract
We show that integrating a polynomial of degree t on an arbitrary simplex (with respect to Lebesgue measure) reduces to evaluating t homogeneous polynomials of degree j = 1, 2,. .. , t, each at a unique point $\xi$ j of the simplex. This new and very simple formula can be exploited in finite (and extended finite) element methods, as well as in other applications where such integrals are needed.
Figures
Reference graph
Works this paper leans on
-
[1]
Antonietti P . F., Houston P ., Pennesi G., Fast Numerical Integration on Polytopic Meshes with Appli- cations to Discontinuous Galerkin Finite Element Methods, J. Sci. Comp., 77, pp. 339–1370 (2018)
work page 2018
-
[2]
How to integrate a polynomial over a simplex, Math
Baldoni V ., Berline N., De Loera J.A., K¨ oppe M., V ergne M. How to integrate a polynomial over a simplex, Math. Comput., 80, pp. 297–325 (2011)
work page 2011
-
[3]
B¨ ueler B., Enge A., Fukuda K., Exact volume computation f or polytopes: A practical study. In Gil Kalai and G¨ unter M. Ziegler, editors, Polytopes – Combinatorics and Computation , volume 29 of DMV Seminar, pages 131–154. Birkh¨ auser V erlag, Basel (2000)
work page 2000
-
[4]
Chin E.B., Lasserre J.B., Sukumar N., Numerical integrat ion of homogeneous functions on convex and nonconvex polygons and polyhedra, Comput. Mech., 56, pp. 967–981 (2015)
work page 2015
-
[5]
Chin E. B., Sukumar N., Modeling curved interfaces withou t element-partitioning in the extended fi- nite element method, Int. J. Num. Methods Eng. to appear. https://doi.org/10.1002/nme.6150
-
[6]
Frenning G., Efficient V oronoi volume estimation for DEM s imulations of granular materials under confined conditions, MethodsX 2, pp. e79-e90 (2015). Simple formula for integration of polynomials on a simplex 11
work page 2015
-
[7]
Kozhasov K., Michalek M., Sturmfels B., Positivity certi ficates via integral representations, arXiv:1908.04191 (2019)
work page Pith review arXiv 2019
-
[8]
Lasserre J.B., Avrachenkov, K., The multi-dimensional v ersion of ∫ b a xpdx. Amer. Math. Monthly, 108, pp. 151–154 (2001)
work page 2001
Show all 14 references
-
[9]
Lasserre J.B., Zeron, E.S., Solving a class of multivaria ble integration problems via Laplace Lech- niques, Appl. Math. (Warsaw), 28, pp. 391–405 (2001)
2001
-
[10]
1126–1140 (2001)
Lasserre J.B., Zeron E.S., A Laplace transform algorith m for the volume of a convex polytope Journal of the ACM, 48, pp. 1126–1140 (2001)
2001
-
[11]
Lasserre, J.B., Integration on a convex polytope, Proc. Am. Math. Soc., 126, pp. 2433–2441(1998)
1998
-
[12]
Leoni F., Y air Shokef Y ., Attraction Controls the Inversion of Order by Disorder in Buckled Colloidal Monolayers, Phys. Rev. Letters (2017)
2017
-
[13]
Methods Appl
Nagy A.P ., Benson D.J., On the numerical integration of t rimmed isogeometric elements, Comput. Methods Appl. Mech. Eng., 284, pp. 165–185 (2015)
2015
-
[14]
Quantitative Spec- troscopy & Radiative Transfer, 129, pp
Y ong Zhang, Y u Ma, Hong-Liang Yin, He-Ping Tan, Natural e lement method for solving radiative transfer with or without conduction in three-dimensional c omplex geometries, J. Quantitative Spec- troscopy & Radiative Transfer, 129, pp. 118–130 (2013)
2013
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.