English

High-Order Quadrature on Multi-Component Domains Implicitly Defined by Multivariate Polynomials

Numerical Analysis 2021-11-24 v2 Numerical Analysis Algebraic Geometry

Abstract

A high-order quadrature algorithm is presented for computing integrals over curved surfaces and volumes whose geometry is implicitly defined by the level sets of (one or more) multivariate polynomials. The algorithm recasts the implicitly defined geometry as the graph of an implicitly defined, multi-valued height function, and applies a dimension reduction approach needing only one-dimensional quadrature. In particular, we explore the use of Gauss-Legendre and tanh-sinh methods and demonstrate that the quadrature algorithm inherits their high-order convergence rates. Under the action of hh-refinement with qq fixed, the quadrature schemes yield an order of accuracy of 2q2q, where qq is the one-dimensional node count; numerical experiments demonstrate up to 22nd order. Under the action of qq-refinement with the geometry fixed, the convergence is approximately exponential, i.e., doubling qq approximately doubles the number of accurate digits of the computed integral. Complex geometry is automatically handled by the algorithm, including, e.g., multi-component domains, tunnels, and junctions arising from multiple polynomial level sets, as well as self-intersections, cusps, and other kinds of singularities. A variety of numerical experiments demonstrates the quadrature algorithm on two- and three-dimensional problems, including: randomly generated geometry involving multiple high-curvature pieces; challenging examples involving high degree singularities such as cusps; adaptation to simplex constraint cells in addition to hyperrectangular constraint cells; and boolean operations to compute integrals on overlapping domains.

Keywords

Cite

@article{arxiv.2105.08857,
  title  = {High-Order Quadrature on Multi-Component Domains Implicitly Defined by Multivariate Polynomials},
  author = {Robert I. Saye},
  journal= {arXiv preprint arXiv:2105.08857},
  year   = {2021}
}

Comments

44 pages, 18 figures, 4 algorithms, 1 table

R2 v1 2026-06-24T02:14:40.252Z