×

Resolution of singularities by rational functions. (English) Zbl 1535.65049

Summary: Results on the rational approximation of functions containing singularities are presented. We build further on the “lightning method,” recently proposed by A. Gopal and L. N. Trefethen [SIAM J. Numer. Anal. 57, No. 5, 2074–2094 (2019; Zbl 1431.65223)], based on exponentially clustering poles close to the singularities. Our results are obtained by augmenting the lightning approximation set with either a low-degree polynomial basis or partial fractions with poles clustering toward infinity in order to obtain a robust approximation of the smooth behavior of the function. This leads to a significant increase in the achievable accuracy as well as the convergence rate of the numerical scheme. For the approximation of \(x^\alpha\) on \([0,1]\), the optimal convergence rate as shown by H. Stahl [Bull. Am. Math. Soc., New Ser. 28, No. 1, 116–122 (1993; Zbl 0768.41015)] is now achieved simply by least-squares fitting.

MSC:

65F20 Numerical solutions to overdetermined systems, pseudoinverses
41A20 Approximation by rational functions
65E05 General theory of numerical methods in complex analysis (potential theory, etc.)
65D05 Numerical interpolation

References:

[1] Adcock, B. and Huybrechs, D., Frames and numerical approximation, SIAM Rev., 61 (2019), pp. 443-473. · Zbl 1421.42015
[2] Adcock, B. and Huybrechs, D., Frames and numerical approximation II: Generalized sampling, J. Fourier Anal. Appl., 26 (2020), pp. 87-114. · Zbl 1462.42048
[3] Baddoo, P. J., Lightning solvers for potential flows, Fluids, 5 (2020), pp. 1-17.
[4] Brubeck, P. D., Nakatsukasa, Y., and Trefethen, L. N., Vandermonde with Arnoldi, SIAM Rev., 63 (2021), pp. 405-415. · Zbl 1484.65022
[5] Brubeck, P. D. and Trefethen, L. N., Lightning Stokes solver, SIAM J. Sci. Comput., 44 (2022), pp. A1205-A1226. · Zbl 1509.41008
[6] Gopal, A. and Trefethen, L. N., New Laplace and Helmholtz solvers, Proc. Natl. Acad. Sci. USA, 116 (2019), pp. 10223-10225. · Zbl 1431.65224
[7] Gopal, A. and Trefethen, L. N., Solving Laplace problems with corner singularities via rational functions, SIAM J. Numer. Anal., 57 (2019), pp. 2074-2094. · Zbl 1431.65223
[8] Newman, D. J., Rational approximation to \(|x|\), Mich. Math. J., 11 (1964), pp. 11-14. · Zbl 0138.04402
[9] Stahl, H., Best uniform rational approximation of \(x^\alpha\) on \([0,1]\), Bull. Amer. Math. Soc., 28 (1993), pp. 116-122. · Zbl 0768.41015
[10] Stahl, H. R., Poles and zeros of best rational approximants of \(|x|\), Constr. Approx., 10 (1994), pp. 469-522. · Zbl 0815.41015
[11] Stahl, H. R., Best uniform rational approximation of \(x^\alpha\) on \([0, 1]\), Acta Math., 190 (2003), pp. 241-306. · Zbl 1077.41012
[12] Stenger, F., Explicit nearly optimal linear rational approximation with preassigned poles, Math. Comput., 47 (1986), pp. 225-252. · Zbl 0592.41019
[13] Trefethen, L. N., Lightning Laplace Solver, https://people.maths.ox.ac.uk/trefethen/lightning.html.
[14] Trefethen, L. N., Approximation Theory and Approximation Practice, extended ed., SIAM, Philadelphia, 2019.
[15] Trefethen, L. N., Numerical conformal mapping with rational functions, Comput. Methods Funct. Theory, 20 (2020), pp. 369-387. · Zbl 1461.30031
[16] Trefethen, L. N., Nakatsukasa, Y., and Weideman, J. A. C., Exponential node clustering at singularities for rational approximation, quadrature, and PDEs, Numer. Math., 147 (2021), pp. 227-254. · Zbl 1467.65021
[17] Trefethen, L. N. and Weideman, J. A. C., The exponentially convergent trapezoidal rule, SIAM Rev., 56 (2014), pp. 385-458. · Zbl 1307.65031
[18] Vyacheslavov, N., Approximation of the function \(|x|\) by rational functions, Math. Notes, 16 (1974), pp. 680-685. · Zbl 0314.41007
[19] Wasow, W., Asymptotic development of the solution of Dirichlet’s problem at analytic corners, Duke Math. J., 24 (1957), pp. 47-56. · Zbl 0077.09802
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. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.