The global extended-rational Arnoldi method for matrix function approximation
Abstract
The numerical computation of matrix functions such as , where is an large and sparse square matrix, is an block with and is a nonlinear matrix function, arises in various applications such as network analysis ( or , machine learning , theory of quantum chromodynamics , electronic structure computation, and others. In this work, we propose the use of global extended-rational Arnoldi method for computing approximations of such expressions. The derived method projects the initial problem onto an global extended-rational Krylov subspace of a low dimension. An adaptive procedure for the selection of shift parameters is given. The proposed method is also applied to solve parameter dependent systems. Numerical examples are presented to show the performance of the global extended-rational Arnoldi for these problems.
Cite
@article{arxiv.2004.00059,
title = {The global extended-rational Arnoldi method for matrix function approximation},
author = {A. H. Bentbib and M. El Ghomari and K. Jbilou},
journal= {arXiv preprint arXiv:2004.00059},
year = {2020}
}