Abstract
We present a stable and efficient Fortran implementation of polynomial interpolation at the Padua points on the square [ − 1,1]
2
. These points are unisolvent and their Lebesgue constant has minimal order of growth (log square of the degree). The algorithm is based on the representation of the Lagrange interpolation formula in a suitable orthogonal basis, and takes advantage of a new matrix formulation together with the machine-specific optimized BLAS subroutine for the matrix-matrix product. Extension to interpolation on rectangles, triangles and ellipses is also described.
Publisher
Association for Computing Machinery (ACM)
Subject
Applied Mathematics,Software
Reference18 articles.
1. AMD 2006. AMD Core Math Library (ACML). Version 3.6.0. Available at http://developer.amd.com/acml.aspx. AMD 2006. AMD Core Math Library (ACML). Version 3.6.0. Available at http://developer.amd.com/acml.aspx.
2. Multivariate simultaneous approximation
3. An updated set of basic linear algebra subprograms (BLAS)
4. On polynomial interpolation of two variables
5. Bivariate interpolation at Xu points: results, extensions and applications;Bos L.;Electron. Trans. Numer. Anal.,2006
Cited by
13 articles.
订阅此论文施引文献
订阅此论文施引文献，注册后可以免费订阅5篇论文的施引文献，订阅后可以查看论文全部施引文献