Bivariate Lagrange interpolation at the node points of non-degenerate Lissajous curves

Motivated by an application in Magnetic Particle Imaging, we study bivariate Lagrange interpolation at the node points of Lissajous curves. The resulting theory is a generalization of the polynomial interpolation theory developed for a node set known as Padua points. With appropriately defined polynomial spaces, we will show that the node points of non-degenerate Lissajous curves allow unique interpolation and can be used for quadrature rules in the bivariate setting. An explicit formula for the Lagrange polynomials allows to compute the interpolating polynomial with a simple algorithmic scheme. Compared to the already established schemes of the Padua and Xu points, the numerical results for the proposed scheme show similar approximation errors and a similar growth of the Lebesgue constant.