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. Mathematics Subject Classification

In this article, we study bivariate polynomial interpolation on the node points of degenerate Lissajous figures. These node points form Chebyshev lattices of rank 1 and are generalizations of the well-known Padua points. We show that these node points allow unique interpolation in appropriately defined spaces of polynomials and give explicit formulas for the Lagrange basis polynomials. Further, we prove mean and uniform convergence of the interpolating schemes. For the uniform convergence the growth of the Lebesgue constant has to be taken into consideration. It turns out that this growth is of logarithmic nature.

No abstract

Accurate interpolation and approximation techniques for functions with discontinuities are key tools in many applications as, for instance, medical imaging. In this paper, we study an RBF type method for scattered data interpolation that incorporates discontinuities via a variable scaling function. For the construction of the discontinuous basis of kernel functions, information on the edges of the interpolated function is necessary. We characterize the native space spanned by these kernel functions and study error bounds in terms of the fill distance of the node set. To extract the location of the discontinuities, we use a segmentation method based on a classification algorithm from machine learning. The conducted numerical experiments confirm the theoretically derived convergence rates in case that the discontinuities are a priori known. Further, an application to interpolation in magnetic particle imaging shows that the presented method is very promising.

No abstract