Abstract
We establish an equivalence between two classes of methods for solving fractional diffusion problems, namely, Reduced Basis Methods (RBM) and Rational Krylov Methods (RKM). In particular, we demonstrate that several recently proposed RBMs for fractional diffusion can be interpreted as RKMs. This changed point of view allows us to give convergence proofs for some methods where none were previously available.
We also propose a new RKM for fractional diffusion problems with poles chosen using the best rational approximation of the function z−s with z ranging over the spectral interval of the spatial discretization matrix. We prove convergence rates for this method and demonstrate numerically that it is competitive with or superior to many methods from the reduced basis, rational Krylov, and direct rational approximation classes. We provide numerical tests for some elliptic fractional diffusion model problems.
Funding statement: The first author has been funded by the Austrian Science Fund (FWF) through grant numbers F 65 and W1245. The second author has been partially supported by the Austrian Science Fund (FWF) grant P 33956-NBL.
Acknowledgment
The constructive comments by the anonymous referees which have led to a marked improvement of the paper are gratefully acknowledged.
References
[1] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Vol. 55, National Bureau of Standards Applied Mathematics Series, 1964.10.1115/1.3625776Search in Google Scholar
[2] L. Aceto, D. Bertaccini, F. Durastante, and P. Novati, Rational Krylov methods for functions of matrices with applications to fractional partial differential equations, J. Comput. Phys. 396 (2019), 470–482.10.1016/j.jcp.2019.07.009Search in Google Scholar
[3] L. Aceto and P. Novati, Rational approximations to fractional powers of self-adjoint positive operators, Numerische Mathematik 143 (2019), No. 1, 1–16.10.1007/s00211-019-01048-4Search in Google Scholar
[4] N. I. Achieser, Theory of Approximation, Dover Books on Advanced Mathematics, Dover Publications, 1992.Search in Google Scholar
[5] H. Antil, Y. Chen, and A. C. Narayan, Reduced basis methods for fractional Laplace equations via extension, SIAM J. Sci. Comput. 41 (2018), A3552–A3575.10.1137/18M1204802Search in Google Scholar
[6] A. Ben-Israel and T. Greville, Generalized Inverses: Theory and Applications, Vol. 15, Springer, New York, 2003.Search in Google Scholar
[7] A. Bonito, D. Guignard, and A. R. Zhang, Reduced basis approximations of the solutions to spectral fractional diffusion problems, J. Numer. Math. 28 (2020), No. 3, 147–160.10.1515/jnma-2019-0053Search in Google Scholar
[8] A. Bonito, W. Lei, and J. E. Pasciak, On sinc quadrature approximations of fractional powers of regularly accretive operators, J. Numer. Math. 27 (2019), No. 2, 57–68.10.1515/jnma-2017-0116Search in Google Scholar
[9] A. Bonito and J. E. Pasciak, Numerical approximation of fractional powers of elliptic operators, Math. Comput. 84 (2015), No. 295, 2083–2110.10.1090/S0025-5718-2015-02937-8Search in Google Scholar
[10] O. Burkovska and M. Gunzburger, Aflne approximation of parametrized kernels and model order reduction for nonlocal and fractional Laplace models, SIAM J. Numer. Anal. 58 (2020), No. 3, 1469–1494.10.1137/19M124321XSearch in Google Scholar
[11] L. Caffarelli and L. Silvestre, An extension problem related to the fractional Laplacian, Commun. Part. Differ. Equ. 32 (2007), No. 8, 1245–1260.10.1080/03605300600987306Search in Google Scholar
[12] T. Chiu, T. Hsieh, C. Huang, and T. Huang, Note on the Zolotarev optimal rational approximation for the overlap Dirac operator, Phys. Review D 66 (2002), No. 11, 114502.10.1103/PhysRevD.66.114502Search in Google Scholar
[13] M. Crouzeix, Numerical range and functional calculus in Hilbert space, J. Funct. Anal. 244 (2007), No. 2, 668–690.10.1016/j.jfa.2006.10.013Search in Google Scholar
[14] T. Danczul, C. Hofreither, and J. Schöberl, A unified rational Krylov method for elliptic and parabolic fractional diffusion problems, arXiv Preprint, 2021, arXiv:2103.13068.Search in Google Scholar
[15] T. Danczul and J. Schöberl, A reduced basis method for fractional diffusion operators, I, Numer. Math. (2022), DOI: 10.1007/s00211-022-01287-y10.1007/s00211-022-01287-ySearch in Google Scholar
[16] T. Danczul and J. Schöberl, A reduced basis method for fractional diffusion operators, II, J. Numer. Math., 29 (2021), No. 4, 269–287.10.1515/jnma-2020-0042Search in Google Scholar
[17] R. DeVore, G. Petrova, and P. Wojtaszczyk, Greedy algorithms for reduced bases in Banach spaces, Construct. Approx. 37 (2013), No. 3, 455–466.10.1007/s00365-013-9186-2Search in Google Scholar
[18] H. Dinh, H. Antil, Y. Chen, E. Cherkaev, and A. Narayan, Model reduction for fractional elliptic problems using Kato’s formula, Mathematical Control & Related Fields, 12 (2022), No. 1, 115–146.10.3934/mcrf.2021004Search in Google Scholar
[19] T. A. Driscoll, N. Hale, and L. N. Trefethen, Chebfun Guide, Pafnuty Publications, 2014.Search in Google Scholar
[20] V. Druskin, L. Knizhnerman, and M. Zaslavsky, Solution of large scale evolutionary problems using rational Krylov subspaces with optimized shifts, SIAM J. Sci. Comput. 31 (2009), No. 5, 3760–3780.10.1137/080742403Search in Google Scholar
[21] V. Druskin, C. Lieberman, and M. Zaslavsky, On adaptive choice of shifts in rational Krylov subspace reduction of evolutionary problems, SIAM J. Sci. Comput. 32 (2010), No. 5, 2485–2496.10.1137/090774082Search in Google Scholar
[22] V. Druskin and V. Simoncini, Adaptive rational Krylov subspaces for large-scale dynamical systems, Systems & Control Letters 60 (2011), No. 8, 546–560.10.1016/j.sysconle.2011.04.013Search in Google Scholar
[23] S.-I. Filip, Y. Nakatsukasa, L. N. Trefethen, and B. Beckermann, Rational minimax approximation via adaptive barycentric representations, SIAM J. Sci. Comput. 40 (2018), No. 4, A2427–A2455.10.1137/17M1132409Search in Google Scholar
[24] A. A. Gonchar, Zolotarëv problems connected with rational functions, Mathematics of the USSR-Sbornik 78 (1969), No. 120, 640–654.Search in Google Scholar
[25] S. Güttel, Rational Krylov Methods for Operator Functions, Ph.D. thesis, Technische Universität Bergakademie Freiberg, Germany, 2010.Search in Google Scholar
[26] S. Güttel, Rational Krylov approximation of matrix functions: Numerical methods and optimal pole selection, GAMM-Mitteilungen 36 (2013), No. 1, 8–31.10.1002/gamm.201310002Search in Google Scholar
[27] S. Güttel and L. Knizhnerman, A black-box rational Arnoldi variant for Cauchy–Stieltjes matrix functions, BIT Numer. Math. 53 (2013), No. 3, 595–616.10.1007/s10543-013-0420-xSearch in Google Scholar
[28] S. Harizanov, R. Lazarov, S.Margenov, and P. Marinov, The best uniform rational approximation: Applications to solving equations involving fractional powers of elliptic operators, arXiv Preprint, 2019, arXiv:1910.13865.Search in Google Scholar
[29] S. Harizanov, R. Lazarov, S.Margenov, P.Marinov, and J. Pasciak, Analysis of numerical methods for spectral fractional elliptic equations based on the best uniform rational approximation, J. Comput. Phys. 408 (2020), 109285.10.1016/j.jcp.2020.109285Search in Google Scholar
[30] S. Harizanov, R. Lazarov, S.Margenov, P.Marinov, and Y. Vutov, Optimal solvers for linear systems with fractional powers of sparse SPD matrices, Numer. Linear Algebra Appl. 25 (2018), No. 5, e2167.10.1002/nla.2167Search in Google Scholar
[31] S. Harizanov, R. Lazarov, P.Marinov, S.Margenov, and J. Pasciak, Comparison analysis on two numerical methods for fractional diffusion problems based on rational approximations of tγ, 0 ⩽ t ⩽ 1, arXiv Preprint, 2018, arXiv:1805.00711.10.1007/978-3-030-14244-5_9Search in Google Scholar
[32] J. S. Hesthaven, G. Rozza, and B. Stamm, Certified Reduced Basis Methods for Parametrized Partial Differential Equations, Springer, Switzerland, 2015.10.1007/978-3-319-22470-1Search in Google Scholar
[33] C. Hofreither, A unified view of some numerical methods for fractional diffusion, Comput. & Math. Appl. 80 (2020), No. 2, 332–350.10.1016/j.camwa.2019.07.025Search in Google Scholar
[34] C. Hofreither, An algorithm for best rational approximation based on barycentric rational interpolation, Numer. Algorithms 88 (2021), No. 1, 365–388.10.1007/s11075-020-01042-0Search in Google Scholar
[35] A. C. Ionită, Lagrange Rational Interpolation and Its Applications to Approximation of Large-Scale ,Dynamical Systems, Ph.D. thesis, Rice University, Houston, TY, 2013.Search in Google Scholar
[36] L. Knizhnerman, V. Druskin, and M. Zaslavsky, On optimal convergence rate of the rational Krylov subspace reduction for electromagnetic problems in unbounded domains, SIAM J. Numer. Anal. 47 (2009), 953–971.10.1137/080715159Search in Google Scholar
[37] A. Lischke, G. Pang, M. Gulian, F. Song, C. Glusa, X. Zheng, Z.Mao, W. Cai, M. M. Meerschaert, M. Ainsworth, and G. E. Karniadakis, What is the fractional Laplacian? A comparative review with new results, J. Comput. Phys. 404 (2020), 109009.10.1016/j.jcp.2019.109009Search in Google Scholar
[38] I. Moret and P. Novati, Krylov subspace methods for functions of fractional differential operators, Math. Comput. 88 (2018), No. 315, 293–312.10.1090/mcom/3332Search in Google Scholar
[39] R. H. Nochetto, E. Otárola, and A. J. Salgado, A PDE approach to fractional diffusion in general domains: A priori error analysis, Found. Comput. Math. 15 (2015), No. 3, 733–791.10.1007/s10208-014-9208-xSearch in Google Scholar
[40] I. V. Oseledets, Lower bounds for separable approximations of the Hilbert kernel, Sbornik: Mathematics, 198 (2007), No. 3, 425–432.10.1070/SM2007v198n03ABEH003842Search in Google Scholar
[41] A. Quarteroni, A.Manzoni, and F. Negri, Reduced Basis Methods for Partial Differential Equations, Springer International Publishing, 2016.10.1007/978-3-319-15431-2Search in Google Scholar
[42] A. Ruhe, Rational Krylov sequence methods for eigenvalue computation, Linear Algebra Appl. 58 (1984), 391–405.10.1016/0024-3795(84)90221-0Search in Google Scholar
[43] J. Schöberl, NETGEN – An advancing front 2D/3D-mesh generator based on abstract rules, Comput. Visual. Sci. 1 (1997), 41–52.10.1007/s007910050004Search in Google Scholar
[44] J. Schöberl, C++11 Implementation of Finite Elements in NGSolve, Institute for Analysis and Scientific Computing, Vienna University of Technology, 2014.Search in Google Scholar
[45] H. R. Stahl, Best uniform rational approximation of xα on [0, 1], Acta Mathematica 190 (2003), No. 2, 241–306.10.1007/BF02392691Search in Google Scholar
[46] P. N. Vabishchevich, Numerically solving an equation for fractional powers of elliptic operators, J. Comput. Phys. 282 (2015), 289–302.10.1016/j.jcp.2014.11.022Search in Google Scholar
[47] P. N. Vabishchevich, Numerical solving unsteady space-fractional problems with the square root of an elliptic operator, Math. Modelling Anal. 21 (2016), No. 2, 220–238.10.3846/13926292.2016.1147000Search in Google Scholar
[48] R. S. Varga and A. J. Carpenter, Some numerical results on best uniform rational approximation of xα on [0, 1], Numer. Algorithms 2 (1992), No. 2, 171–185.10.1007/BF02145384Search in Google Scholar
[49] D. R. Witman, M. Gunzburger, and J. Peterson, Reduced-order modeling for nonlocal diffusion problems, Int. J. Numer. Methods Fluids 83 (2016), No. 3, 307–327.10.1002/fld.4269Search in Google Scholar
[50] E. I. Zolotarëv, Collected works, St. Petersburg Academy of Sciences, St.-Petersburg, 1877 (in Russian).Search in Google Scholar
© 2022 Walter de Gruyter GmbH, Berlin/Boston
