Barycentric Hermite interpolation

From MaRDI portal
Publication:2847712




Abstract: Let z1,ldots,zK be distinct grid points. If fk,0 is the prescribed value of a function at the grid point zk, and fk,r the prescribed value of the rforeignlanguage{american}{-th} derivative, for 1leqrleqnk1, the Hermite interpolant is the unique polynomial of degree N1 (N=n1+cdots+nK) which interpolates the prescribed function values and function derivatives. We obtain another derivation of a method for Hermite interpolation recently proposed by Butcher et al. {[}emph{Numerical Algorithms, vol. 56 (2011), p. 319-347}{]}. One advantage of our derivation is that it leads to an efficient method for updating the barycentric weights. If an additional derivative is prescribed at one of the interpolation points, we show how to update the barycentric coefficients using only mathcalOleft(Night) operations. Even in the context of confluent Newton series, a comparably efficient and general method to update the coefficients appears not to be known. If the method is properly implemented, it computes the barycentric weights with fewer operations than other methods and has very good numerical stability even when derivatives of high order are involved. We give a partial explanation of its numerical stability.









This page was built for publication: Barycentric Hermite interpolation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2847712)