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 $h$-refinement with $q$ fixed, the quadrature schemes yield an order of accuracy of $2q$, where $q$ is the one-dimensional node count; numerical experiments demonstrate up to 22nd order. Under the action of $q$-refinement with the geometry fixed, the convergence is approximately exponential, i.e., doubling $q$ 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.
翻译:高阶二次曲线算法用于计算曲线表面和体积的积分,其几何由(一个或一个以上)多变多元数数的等级组合暗含地基。 算法将隐含地定义的几何转换为隐含地定义的多值高度函数的图形, 并应用一个仅需要一维方位方形的减少维度方法。 特别是, 我们探索高- Legendre 和 tanh- sinh 方法的使用, 并显示二次方程算法继承了它们的高度趋同率。 在以美元固定的( 一个或一个以上)多变量组合中, 二次方程计划产生的精确度为2q$, 其中美元为一维值的高度函数计数; 数字实验显示到22度。 在以美元计方位定的精度动作下, 汇合大约是美元双倍的。 复杂的几何直方位直方位值由三个运算法( 包括: e. g. drelexal deal deal deal deal) 等数、 多度的内径、 多立度的内径、 多立度的内径、 多立、 多立、 多立、 多立度的内径的内径、 多立、 多立、 多立、多立、多立、多立的内径的内径的内径的内径、多立、多立、多立、多立、多立、多立、多立、多立、多立、多立、多立、多立、多立、多的内径的内径的内径、多立、多立、多立、多的内径的内径的内径的内径的内径的内径的内径、多的内、多的内径、多等的内径、多的内、多的内、多的内、多的内、多的内、多的内、多的内、多的内、多的内径、多的内径、多的内、多的内、多的内、多形、多的内、多的内、多的内、多的内、多的内、多形、多的内、多形的内、多的内、多的内、多的内、多的内、多