High-order integration on regular triangulated manifolds reaches superalgebraic approximation rates through cubical reparametrizations (Q6877530)

This interesting paper studies high-order integration on regular triangulated manifolds. More precisely, let \(S\) be a compact, orientable, \(d\)-dimensional \(C^{r+1},\ r>0\) manifold embedded into some \(0\leq d\leq m\)-dimensional Euclidean space \(\mathbb R^d\) and let \(f:S\to \mathbb R\) be integrable with \(f\in C^r\). The authors are interested in algorithms to approximate integrals of the form \(\int_{S}fdS\). Such integrals arise in many applications for example in geometric processing, surface-interface and colloidal sciences, various optimization of production processes, spectral methods to solve PDEs on curved surfaces and many others. The authors provide methods for deriving high-order volume elements for such integrals. To do this, they introduce a homeomorphic multilinear hypercube-simplex transformation which reparametrizes an initial flat triangulation of the manifold to a certain cubical mesh. By using square-squeezing, the authors approximate the integrand and the volume element for each hypercube domain of the reparametrized mesh through interpolation in Chebyshev-Lobatto grids. The authors obtain bounds for the integration error in terms of the \(r\)th-order total variation of the integrand and the surface parametrization, predicting high algebraic approximation rates that scale solely with the interpolation degree. For integrands whose total variation is constantly bounded with increasing \(r\), the estimates show that the integration error is decreasing exponentially, while mesh refinements are limited to achieve algebraic rates. Numerous numerics are given to support the theoretical findings. \medskip\N\NThe paper is very well written with a very good set of references.