Extended and rational Hessenberg methods for the evaluation of matrix functions

被引:2
作者
Ramezani, Z. [1 ]
Toutounian, F. [1 ,2 ]
机构
[1] Ferdowsi Univ Mashhad, Sch Math Sci, Dept Appl Math, Mashhad, Razavi Khorasan, Iran
[2] Ferdowsi Univ Mashhad, Ctr Excellence Modeling & Control Syst, Mashhad, Razavi Khorasan, Iran
关键词
Krylov subspace methods; Extended Krylov subspace; Rational Krylov subspace; Hessenberg process; Matrix function; Shifted linear system; 65F10; NONSYMMETRIC LINEAR-SYSTEMS; KRYLOV SUBSPACE APPROXIMATIONS; BLOCK CMRH METHOD; RESTARTED GMRES; SQUARE-ROOT; ALGORITHM; IMPLEMENTATION; EQUATIONS;
D O I
10.1007/s10543-018-0742-9
中图分类号
TP31 [计算机软件];
学科分类号
081202 ; 0835 ;
摘要
Some Krylov subspace methods for approximating the action of matrix functions are presented in this paper. The main idea of these techniques is to project the approximation problem onto a subspace of much smaller dimension. Then the matrix function operation is performed with a much smaller matrix. These methods are projection methods that use the Hessenberg process to generate bases of the approximation spaces. We also use the introduced methods to solve shifted linear systems. Some numerical experiments are presented in order to show the efficiency of the proposed methods.
引用
收藏
页码:523 / 545
页数:23
相关论文
共 73 条
[1]  
Afanasjew M, 2007, ELECTRON T NUMER ANA, V28, P206
[2]   Numerical approximation of the product of the square root of a matrix with a vector [J].
Allen, EJ ;
Baglama, J ;
Boyd, SK .
LINEAR ALGEBRA AND ITS APPLICATIONS, 2000, 310 (1-3) :167-181
[3]   Weighted and flexible versions of block CMRH method for solving nonsymmetric linear systems with multiple right-hand sides [J].
Amini, S. ;
Toutounian, F. .
COMPUTERS & MATHEMATICS WITH APPLICATIONS, 2018, 76 (08) :2011-2021
[4]   The block CMRH method for solving nonsymmetric linear systems with multiple right-hand sides [J].
Amini, S. ;
Toutounian, F. ;
Gachpazan, M. .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2018, 337 :166-174
[5]  
[Anonymous], 2008, Functions of matrices
[6]  
[Anonymous], MATR MARK
[7]   Superlinear convergence of the rational Arnoldi method for the approximation of matrix functions [J].
Beckermann, Bernhard ;
Guettel, Stefan .
NUMERISCHE MATHEMATIK, 2012, 121 (02) :205-236
[8]   ERROR ESTIMATES AND EVALUATION OF MATRIX FUNCTIONS VIA THE FABER TRANSFORM [J].
Beckermann, Bernhard ;
Reichel, Lothar .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2009, 47 (05) :3849-3883
[9]  
Benzi M, 2007, ELECTRON T NUMER ANA, V28, P16
[10]   Lanczos-based exponential filtering for discrete ill-posed problems [J].
Calvetti, D ;
Reichel, L .
NUMERICAL ALGORITHMS, 2002, 29 (1-3) :45-65