A rational radial basis function method for accurately resolving discontinuities and steep gradients. (English) Zbl 1397.65209

Summary: Radial Basis Function (RBF) methods have become important tools for scattered data interpolation and for solving partial differential equations (PDEs) in complexly shaped domains. When the underlying function is sufficiently smooth, RBF methods can produce exceptional accuracy. However, like other high order numerical methods, if the underlying function has steep gradients or discontinuities the RBF method may/will produce solutions with non-physical oscillations. In this work, a rational RBF method is used to approximate derivatives of functions with steep gradients and discontinuities and to solve PDEs with such solutions. The method is non-linear and is more computationally expensive than the standard RBF method. A modified partition of unity method is discussed as an way to implement the rational RBF method in higher dimensions.


65M70 Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs
65D05 Numerical interpolation
41A05 Interpolation in approximation theory
41A20 Approximation by rational functions
65F15 Numerical computation of eigenvalues and eigenvectors of matrices
Full Text: DOI


[1] Babuka, I.; Melenk, J. M., The partition of unity method, Int. J. Numer. Methods Eng., 40, 4, 727-758, (1997) · Zbl 0949.65117
[2] Buhmann, M. D., Radial basis functions, (2003), Cambridge University Press · Zbl 1038.41001
[3] Clenshaw, C. W.; Lord, K., Rational approximations from Chebyshev series, (Scaife, B. K.P., Studies in Numerical Analysis, (1974), Academic Press), 95-113 · Zbl 0308.65010
[4] Deparis, S.; Forti, D.; Quarteroni, A., A rescaled localized radial basis function interpolation on non-Cartesian and nonconforming grids, SIAM J. Sci. Comput., 36, 6, A2745-A2762, (Jan. 2014)
[5] Fasshauer, G.; McCourt, M., Stable evaluation of Gaussian RBF interpolants, SIAM J. Sci. Comput., 34, 737-762, (2012)
[6] Fasshauer, G. E., Meshfree approximation methods with Matlab, (2007), World Scientific · Zbl 1123.65001
[7] Fornberg, B.; Larsson, E.; Flyer, N., Stable computations with Gaussian radial basis functions, SIAM J. Sci. Comput., 33, 869-892, (2011) · Zbl 1227.65018
[8] Franke, R., Locally determined smooth interpolation at irregularly spaced points in several variables, IMA J. Appl. Math., 19, 4, 471-482, (1977) · Zbl 0361.65008
[9] Hesthaven, J.; Kaber, S.; Lurati, L., Padé-Legendre interpolants for Gibbs reconstruction, J. Sci. Comput., 28, 2-3, 337-359, (2006) · Zbl 1102.65011
[10] Jakobsson, S.; Andersson, B.; Edelvik, F., Rational radial basis function interpolation with applications to antenna design, J. Comput. Appl. Math., 233, 889-904, (2009) · Zbl 1178.65009
[11] Lehoucq, R. B.; Sorensen, D. C.; Yang, C., ARPACK users guide, (Jan. 1998), Society for Industrial & Applied Mathematics (SIAM)
[12] De Marchi, S.; Idda, A.; Santin, G., A rescaled method for RBF approximation, (Approximation Theory XV, San Antonio 2016, (2017), Springer International Publishing), 39-59 · Zbl 1385.65017
[13] MATLAB, (2016), The MathWorks Inc. Natick, Massachusetts, version 9.1 (R2016b)
[14] Min, M. S.; Kaber, S. M.; Don, W. S., Fourier-Padé approximations and filtering for spectral simulations of an incompressible Boussinesq convection problem, Math. Comput., 76, 1275-1290, (2007) · Zbl 1122.41007
[15] Nakatsukasa, Y.; Sete, O.; Trefethen, L. N., The AAA algorithm for rational approximation, Preprint · Zbl 1390.41015
[16] Niederreiter, H., Random number generation and quasi-Monte Carlo methods, (1992), CBMS-NSF, SIAM Philadelphia · Zbl 0761.65002
[17] Safdari-Vaighani, A.; Heryudono, A.; Larsson, E., A radial basis function partition of unity collocation method for convection-diffusion equations arising in financial applications, J. Sci. Comput., 64, 2, 341-367, (Oct. 2014)
[18] Sarra, S. A., A Matlab radial basis function toolkit with symmetry exploitation, regularization, and extended precision
[19] Sarra, S. A., Regularized symmetric positive definite matrix factorizations for linear systems arising from RBF interpolation and differentiation, Eng. Anal. Bound. Elem., 44, 76-86, (2014) · Zbl 1297.65032
[20] Sarra, S. A., The Matlab radial basis function toolbox, J. Open Res. Softw., 5, (March 2017)
[21] Sarra, S. A.; Cogar, S., An examination of evaluation algorithms for the RBF method, Eng. Anal. Bound. Elem., 75, 36-45, (Feb. 2017)
[22] Sarra, S. A.; Kansa, E. J., Multiquadric radial basis function approximation methods for the numerical solution of partial differential equations, Adv. Comput. Mech., 2, (2009)
[23] Schaback, R., Error estimates and condition numbers for radial basis function interpolation, Adv. Comput. Math., 3, 251-264, (1995) · Zbl 0861.65007
[24] Shepard, D., A two-dimensional interpolation function for irregularly-spaced data, (Proceedings of the 1968 23rd ACM National Conference on, (1968), Association for Computing Machinery (ACM))
[25] Tadmor, E., Convergence of spectral methods for nonlinear conservation laws, SIAM J. Numer. Anal., 26, 30-44, (1989) · Zbl 0667.65079
[26] Wendland, H., Piecewise polynomial, positive definite and compactly supported radial basis functions of minimal degree, Adv. Comput. Math., 4, 389-396, (1995) · Zbl 0838.41014
[27] Wendland, H., Fast evaluation of radial basis functions: methods based on partition of unity, (Approximation Theory X: Wavelets, Splines, and Applications, (2002), Vanderbilt University Press), 473-483 · Zbl 1031.65022
[28] Wendland, H., Scattered data approximation, (2005), Cambridge University Press · Zbl 1075.65021
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.