Polynomial finite element method for domains enclosed by piecewise conics

We consider bivariate piecewise polynomial finite element spaces for curved domains bounded by piecewise conics satisfying homogeneous boundary conditions, construct stable local bases for them using Bernstein–Bézier techniques, prove error bounds and develop optimal assembly algorithms for the finite element system matrices. Numerical experiments confirm the effectiveness of the method.