Abstract
AbstractIn this paper we propose a method to compute the solution to the fractional diffusion equation on directed networks, which can be expressed in terms of the graph Laplacian L as a product $$f(L^T) \varvec{b}$$
f
(
L
T
)
b
, where f is a non-analytic function involving fractional powers and $$\varvec{b}$$
b
is a given vector. The graph Laplacian is a singular matrix, causing Krylov methods for $$f(L^T) \varvec{b}$$
f
(
L
T
)
b
to converge more slowly. In order to overcome this difficulty and achieve faster convergence, we use rational Krylov methods applied to a desingularized version of the graph Laplacian, obtained with either a rank-one shift or a projection on a subspace.
Publisher
Springer Science and Business Media LLC
Subject
Applied Mathematics,Computational Mathematics,Computer Networks and Communications,Software
Reference33 articles.
1. Beckermann, B., Reichel, L.: Error estimates and evaluation of matrix functions via the Faber transform. SIAM J. Numer. Anal. 47(5), 3849–3883 (2009). https://doi.org/10.1137/080741744
2. Benzi, M., Bertaccini, D., Durastante, F., Simunec, I.: Non-local network dynamics via fractional graph Laplacians. J. Complex Netw. 8(3), cnaa017 (2020). https://doi.org/10.1093/comnet/cnaa017
3. Benzi, M., Fika, P., Mitrouli, M.: Graphs with absorption: numerical methods for the absorption inverse and the computation of centrality measures. Linear Algebra Appl. 574, 123–152 (2019). https://doi.org/10.1016/j.laa.2019.03.026
4. Berljafa, M., Elsworth, S., Güttel, S.: A rational Krylov toolbox for MATLAB. MIMS EPrint 2014.56, Manchester Institute for Mathematical Sciences, University of Manchester, Manchester, UK (2014). http://rktoolbox.org
5. Berljafa, M., Güttel, S.: Generalized rational Krylov decompositions with an application to rational approximation. SIAM J. Matrix Anal. Appl. 36(2), 894–916 (2015). https://doi.org/10.1137/140998081
Cited by
7 articles.
订阅此论文施引文献
订阅此论文施引文献，注册后可以免费订阅5篇论文的施引文献，订阅后可以查看论文全部施引文献