Quadrature for Implicitly-defined Finite Element Functions on Curvilinear Polygons

$H^1$-conforming Galerkin methods on polygonal meshes such as VEM, BEM-FEM and Trefftz-FEM employ local finite element functions that are implicitly defined as solutions of Poisson problems having polynomial source and boundary data. Recently, such methods have been extended to allow for mesh cells that are curvilinear polygons. Such extensions present new challenges for determining suitable quadratures. We describe an approach for integrating products of these implicitly defined functions, as well as products of their gradients, that reduces integrals on cells to integrals along their boundaries. Numerical experiments illustrate the practical performance of the proposed methods.