Low cardinality positive interior cubature on NURBS-shaped domains

DOI: https://doi.org/10.1007/s10543-023-00958-y
IF: 1.856
2023-03-18
BIT Numerical Mathematics
Abstract:In this paper we propose an algorithm that computes a low cardinality PI-type (Positive weights and Interior nodes) algebraic cubature rule of degree n , with at most nodes, over curvilinear polygons defined by piecewise rational functions. Typical examples are domains whose boundary is defined piecewise by NURBS curves or by composite Bezier curves. The key tools are a relevant but overlooked theorem of 1976 by Wilhelmsen on Tchakaloff sets , a specific in-domain algorithm for such curvilinear polygons and the sparse nonnegative solution of underdetermined moment matching systems by the Lawson–Hanson Nonnegative least squares solver. Many numerical tests are performed to show the flexibility of this approach and the implemented MATLAB toolbox is freely available to the users, in particular for possible applications within FEM/VEM with NURBS-shaped curvilinear elements.
mathematics, applied,computer science, software engineering
What problem does this paper attempt to address?