Abstract
Both barycentric Lagrange interpolation and barycentric rational interpolation are thought to be stable and effective methods for approximating a given function on some special point sets. A direct evaluation of these interpolants due to N interpolation points at M sampling points requires \(\mathcal {O}(NM)\) arithmetic operations. In this paper, we introduce two fast multipole methods to reduce the complexity to \(\mathcal {O}(\max \left \{N,M\right \})\). The convergence analysis is also presented in this paper.
Similar content being viewed by others
References
Abramowitz, M., Stegun, I.A.: Handbook of mathematical functions. National Bureau of Standards, Washington DC (1964)
Baker, R.D., Jackson, D.: Statistical application of barycentric rational interpolants: an alternative to splines. Comput. Stat. 29, 1065–1081 (2014)
Baltensperger, R., Berrut, J.P., Noël, B.: Exponential convergence of a linear rational interpolant between transformed Chebyshev points. Math. Comp. 12, 1109–1120 (1999)
Bernstein, S.N.: Sur l’ordre de la meilleure approximation des fonctions continues par les polynômes de degré donné. Mém. Acad. Roy Belg. 4, 1–103 (1912)
Berrut, J.P.: Rational functions for guaranteed and experimentally well-conditioned global interpolation. Comput. Math. Appl. 15, 1–16 (1988)
Berrut, J.P., Baltensperger, R., Mittelmann, H.D.: Recent developments in barycentric rational interpolation. In: de Bruin, M. G., Mache, D. H., Szabados, J. (eds.) Trends and applications in constructive approximation. International Series of Numerical Mathematics, vol. 151, pp. 27–51. Basel, Birkhäuser (2005)
Berrut, J.P., Floater, M.S., Klein, G.: Convergence rates of derivatives of a family of barycentric rational interpolants. Appl. Numer. Math. 61, 989–1000 (2011)
Berrut, J.P., Trefethen, L.N.: Barycentric Lagrange interpolation. SIAM Rev. 46, 501–517 (2004)
Boyd, J.P.: Computing zeros on a real interval through Chebyshev expansion and polynomial rootfinding. SIAM J. Numer. Anal. 40, 1666–1682 (2002)
Boyd, J.P.: Finding the zeros of a univariate equation: proxy rootfinders, Chebyshev interpolation, and the companion matrix. SIAM Rev. 55, 375–396 (2013)
Carrier, J., Greengard, L., Rokhlin, V.: A fast adaptive multipole algorithm for particle simulations. SIAM J. Sci. Stat. Comput. 9, 669–686 (1988)
Dutt, A., Gu, M., Rokhlin, V.: Fast algorithms for polynomial interpolation, integration, and differentiation. SIAM J. Numer. Anal. 33, 1689–1711 (1996)
Floater, M.S., Hormann, K.: Barycentric rational interpolation with no poles and high rates of approximation. Numer. Math. 107, 315–331 (2007)
Gentleman, W.M.: Implementing Clenshaw-Curtis quadrature. Comm. ACM 15, 337–346 (1972)
Gentleman, W.M.: Algorithm 424: Clenshaw-Curtis quadrature [D1]. Comm. ACM 15, 353–355 (1972)
Good, I.J.: The colleague matrix, a Chebyshev analogue of the companion matrix. Quart. J. Math. 12, 61–68 (1961)
Greengard, L., Rokhlin, V.: A fast algorithm for particle simulations. J. Comput. Phys. 73, 325–348 (1987)
Greengard, L., Rokhlin, V.: A new version of the fast multipole method for the Laplace equation in three dimensions. Acta Numer. 6, 229–269 (1997)
Hale, N., Trefethen, L.N.: Chebfun and numerical quadrature. Sci. China 55, 1749–1760 (2012)
Higham, N.J.: The numerical stability of barycentric Lagrange interpolation. IMA J. Numer. Anal. 24, 547–556 (2004)
Klein, G., Berrut, J.: Linear barycentric rational quadrature. BIT 52, 407–424 (2012)
Liu, Y.J.: Fast multipole boundary element method: theory and applications in engineering. Cambridge University Press, Cambridge (2009)
Martinsson, P.G., Rokhlin, V.: An accelerated kernel-independent fast multipole method in one dimension. SIAM J. Sci. Comput. 29, 1160–1178 (2007)
Runge, C.: Über empirische Funktionen und die Interpolation zwischen äquidistanten Ordinaten. Z Math. Phys. 46, 224–243 (1901)
Salzer, H.E.: Lagrange interpolation at the Chebyshev points \(x_{n},\nu =\cos (\pi /n), \nu =0(1)n\); some unnoted advantages. Comput. J. 15, 156–159 (1972)
Schwarz, H.R.: Numerische mathematik, 4th edn. Teubner, Stuttgart (1997)
Specht, W.: Die lage der nullstellen eines polynoms. Math. Nachr. 21, 201–222 (1960)
Taylor, W.J.: Method of Lagrangian curvilinear interpolation. J. Research. Natl. Bur. Stand. 35, 151–159 (1945)
Trefethen, L.N.: Is Gauss quadrature better than Clenshaw-Curtis?. SIAM Rev. 50, 67–87 (2008)
Trefethen, L.N.: Approximation theory and approximation practice. SIAM, Philadelphia (2013)
Trefethen, L.N., et al.: Chebfun Version 4.0, The Chebfun Development Team, http://www.maths.ox.ac.uk/chebfun/ (2011)
Wang, H., Huybrechs, D., Vandewalle, S.: Explicit barycentric weights for polynomial interpolation in the roots or extrema of classical orthogonal polynomials. Math. Comp. 83, 2893–2914 (2014)
Wang, H., Xiang, S.: On the convergence rates of Legendre approximation. Math. Comp. 81, 861–877 (2012)
Xiang, S.: On interpolation approximation: convergence rates for polynomial interpolation for functions of limited regularity. SIAM J. Numer. Anal. 54, 2081–2113 (2016)
Xiang, S., Chen, X., Wang, H.: Error bounds for approximation in Chebyshev points. Numer. Math. 116, 463–491 (2010)
Yarvin, N., Rokhlin, V.: An improved fast multipole algorithm for potential fields on the line. SIAM J. Numer. Anal. 36, 629–666 (1999)
Zhang, B., Huang, J., Pitsianis, N.P., Sun, X.: RECFMM: Recursive parallelization of the adaptive fast multipole method for coulomb and screened coulomb interactions. Commun. Comput. Phys. 20, 534–550 (2016)
Author information
Authors and Affiliations
Corresponding author
Additional information
This work is supported partly by NSF of China (No.11371376), the Innovation-Driven Project and the Mathematics and Interdisciplinary Sciences Project of Central South University.
Rights and permissions
About this article
Cite this article
Liu, G., Xiang, S. Fast multipole methods for approximating a function from sampling values. Numer Algor 76, 727–743 (2017). https://doi.org/10.1007/s11075-017-0279-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-017-0279-z