A new investigation of the extended Krylov subspace method for matrix function evaluations

Numerical Linear Algebra with Applications

Volumn 17, Issue 4, 2010, Pages 615-638

  Knizhnerman, L ^a   Simoncini, V ^b,c  

^a Central Geophysical Expedition ^*  (Russian Federation)

^b UNIVERSITY OF BOLOGNA   (Italy)

^c CIRSA   (Italy)

Author keywords

Faber polynomials; Krylov subspaces; Large matrices; Matrix functions

Indexed keywords

EID: 77954843781     PISSN: 10705325     EISSN: 10991506     Source Type: Journal    
DOI: 10.1002/nla.652     Document Type: Article

References (55)

1. Higham NJ. Matrix Functions-Theory and Applications. SIAM: Philadelphia, U.S.A., 2008.
2. Friesner RA, Tuckerman LS, Dornblaser BC, Russo TV. A method for exponential propagation of large systems of stiff nonlinear differential equations. Journal of Scientific Computing 1989; 4:327-354.
3. Puzynin IV, Selin AV, Vinitsky SI. Magnus-factorized method for numerical solving the time-dependent Schr̈odinger equation. Computer Physics Communications 2000; 126(1-2):158-161.
4. Gallopoulos E, Saad Y. Efficient solution of parabolic equations by Krylov approximation methods. SIAM Journal on Scientific and Statistical Computing 1992; 13(5):1236-1264.
5. Gallopoulos E, Saad Y. A parallel block cyclic reduction algorithm for the fast solution of elliptic equations. Parallel Computing 1989; 10:143-159.
6. Hochbruck M, Ostermann A. Exponential Runge-Kutta methods for parabolic problems. Applied Numerical Mathematics 2005; 53:323-339.
7. Hochbruck M, Lubich C. Exponential integrators for quantum-classical molecular dynamics. BIT 1999; 39(1): 620-645.
8. Hochbruck M, Lubich C, Selhofer H. Exponential integrators for large systems of differential equations. SIAM Journal on Scientific Computing 1998; 19(5):1552-1574.
9. Borici A. Computational methods for UV suppresses fermions. Journal of Computational Physics 2003; 189: 454-462.
10. Nettesheim P, Scḧutte Ch. Numerical integrators for quantum-classical molecular dynamics. In Computational Molecular Dynamics: Challenges, Methods, Ideas, Deuflhard P, Hermans J, Leimkuhler B, Marks AE, Reich S, Skeel RD (eds). Lecture Notes in Computational Science and Engineering, vol. 4. Springer: Berlin, 1999; 396-411.
11. van den Eshof J, Frommer A, Lippert Th, Schilling K, van der Vorst HA. Numerical methods for the QCD overlap operator. I: sign-function and error bounds. Computer Physics Communications 2002; 146(2): 203-224.
12. Druskin VL, Knizhnerman LA, Lee P. New spectral Lanczos decomposition method for induction modeling in arbitrary 3-D geometry. Geophysics 1999; 64(3):701-706.
13. Druskin V, Knizhnerman L. Two polynomial methods of calculating functions of symmetric matrices. U.S.S.R. Computational Mathematics and Mathematical Physics 1989; 29:112-121.
14. Druskin V, Knizhnerman L. Krylov subspace approximation of eigenpairs and matrix functions in exact and computer arithmetic. Numerical Linear Algebra with Applications 1995; 2(3):205-217.
15. Knizhnerman L. Calculation of functions of unsymmetric matrices using Arnoldi's method. U.S.S.R. Computational Mathematics and Mathematical Physics 1991; 31:1-9.
16. Saad Y. Analysis of some Krylov subspace approximations to the matrix exponential operator. SIAM Journal on Numerical Analysis 1992; 29:209-228.
17. Bergamaschi L, Vianello M. Efficient computation of the exponential operator for large, sparse, symmetric matrices. Numerical Linear Algebra with Applications 2000; 7(1):27-45.
18. Druskin V, Knizhnerman L. Extended Krylov subspaces: approximation of the matrix square root and related functions. SIAM Journal on Matrix Analysis and Applications 1998; 19(3):755-771.
19. Eiermann M, Ernst O. A restarted Krylov subspace method for the evaluation of matrix functions. SIAM Journal on Numerical Analysis 2006; 44(6):2481-2504.
20. Frommer A, Simoncini V. Stopping criteria for rational matrix functions of Hermitian and symmetric matrices. SIAM Journal on Scientific Computing 2008; 30(3):1387-1412.
21. Ilic M, Simpson DP, Turner IW. A restarted Lanczos approximation to functions of a symmetric matrix. Technical Report 8011, Queensland University of Technology, Australia, 2007.
22. Moret I. Rational Lanczos approximations to the matrix square root and related functions. Technical Report, Dipartimento di Matematica e Informatica, Università degli Studi di Trieste, Trieste, Italy, 2005. To appear in Numerical Linear Algebra with Applications.
23. Moret I, Novati P. An interpolatory approximation of the matrix exponential based on Faber polynomials. Journal of Computational and Applied Mathematics 2001; 131:361-380.
24. Moret I, Novati P. RD-rational approximations of the matrix exponential. BIT 2004; 44(3):595-615.
25. Popolizio M, Simoncini V. Acceleration techniques for approximating the matrix exponential operator. SIAM Journal on Matrix Analysis and Applications 2008; 30(2):657-683.
26. Trefethen LN, Weideman JAC, Schmelzer T. Talbot quadratures and rational approximations. BIT 2006; 46(3): 653-670.
27. Simoncini V. A new iterative method for solving large-scale Lyapunov matrix equations. SIAM Journal on Scientific Computing 2007; 29(3):1268-1288.
28. Frommer A. Monotone convergence of the Lanczos approximations to matrix functions of Hermitian matrices. Technical Report, Fachbereich Universiẗat, Wuppertal-D, 2008.
29. Popolizio M. Acceleration techniques for approximating the matrix exponential. Ph.D. Thesis, Università degli Studi di Bari, March 2008.
30. Alpak FO, Torres-Verd ́i{dotless}n C, Sepehrnoori K, Fang S, Knizhnerman L. An extended Krylov subspace method to simulate single-phase fluid flow phenomena in axisymmetric and anisotropic porous media. Journal of Petroleum Science and Engineering 2003; 40:121-144.
31. Nauts A, Wyatt RE. New approach to many state quantum dynamics: the recursive residue generation method. Physical Review Letters 1983; 51:2238-2241.
32. Tal-Ezer H. Spectral methods in time for parabolic problems. SIAM Journal on Numerical Analysis 1989; 26: 1-11.
33. Horn RA, Johnson CR. Topics in Matrix Analysis. Cambridge University Press: Cambridge, 1991.
34. Druskin V, Greenbaum A, Knizhnerman L. Using nonorthogonal Lanczos vectors in the computation of matrix functions. SIAM Journal on Scientific Computing 1998; 19(1):38-54.
35. van den Eshof J, Hochbruck M. Preconditioning Lanczos approximations to the matrix exponential. SIAM Journal on Scientific Computing 2006; 27(4):1438-1457.
36. Druskin VL. On monotonicity of the Lanczos approximation to the matrix exponential. Linear Algebra and its Applications 2008; 429(7):1679-1683.
37. Druskin V, Knizhnerman L. Spectral approach to solving three-dimensional Maxwell's diffusion equations in the time and frequency domains. Radio Science 1994; 29(4):937-953.
38. Suetin PK. Series of Faber polynomials. Analytical Methods and Special Functions. Gordon and Breach Science: Amsterdam, 1998. Originally published in Russian in 1984 as Riady po mnogochlennam Fabera by Nauka, Moscow.
39. Crouzeix M. Numerical range and numerical calculus in Hilbert space. Journal of Functional Analysis 2007; 244:668-690.
40. Beckermann B. Image num ́erique, GMRES et polyn̂omes de Faber. Comptes Rendus Mathematique, Acad ́emie de Science Paris, Series I 2005; 340:855-860.
41. Beckermann B, Reichel L. Error estimation and evaluation of matrix functions via the Faber transform. Manuscript, December 2008.
42. Gustafson K, Rao DKM. Numerical Range: The Field of Values of Linear Operators and Matrices. Universitext. Springer: New York, NY, U.S.A., 1977.
43. Higham NJ. The Matrix Computation Toolbox. Available from: http://www.ma.man.ac.uk/higham/mctoolbox.
44. Driscoll T. Algorithm 756: a MATLAB toolbox for Schwarz-Christoffel mapping. ACM Transactions on Mathematical Software 1996; 22(2):168-186.
45. Faber V, Greenbaum A, Marshall DE. The polynomial numerical hulls of Jordan blocks and related matrices. Linear Algebra and its Applications 2003; 374:231-246.
46. Knizhnerman L. Adaptation of the Lanczos and Arnoldi methods to the spectrum, or why the two Krylov subspace methods are powerful. Chebyshev Digest (in Russian: Chebyshyovsky Sbornik) 2002; 3(2):141-164. Also available from: http://www.tspu.tula.ru/res/math/csbor/tom3/v2/knizhnerman.rar
47. The MathWorks, Inc. MATLAB 7, September 2004.
48. Diele F, Moret I, Ragni S. Error estimates for polynomial Krylov approximations to matrix functions. SIAM Journal on Matrix Analysis and Applications 2008; 30(4):1546-1565.
49. Petrushev PP, Popov VA. Rational Approximation of Real Functions. Cambridge University Press: Cambridge, 1987.
50. Matrix Market. A visual repository of test data for use in comparative studies of algorithms for numerical linear algebra, Mathematical and Computational Sciences Division, National Institute of Standards and Technology. Available Online at: http://math.nist.gov/MatrixMarket.
51. Borici A, Frommer A, Jo ́o B, Kennedy A, Pendleton B. QCD and numerical analysis III. Proceedings of the Third International Workshop on Numerical Analysis and Lattice QCD, Edinburgh, U.K., 30 June-4 July 2003. Lecture Notes in Computational Science and Engineering, vol. 47, xii. Springer: Berlin, 2005; 201.
52. Simoncini V, Szyld DB. Theory of inexact Krylov subspace methods and applications to scientific computing. SIAM Journal on Scientific Computing 2003; 25(2):454-477.
53. Lopez L, Simoncini V. Analysis of projection methods for rational function approximation to the matrix exponential. SIAM Journal on Numerical Analysis 2006; 44(2):613-635.
54. Sonneveld P, van Gijzen MB. IDR(s): a family of simple and fast algorithms for solving large nonsymmetric linear systems. SIAM Journal on Scientific Computing 2008; 31(2):1035-1062.
55. Boyle J, Mihajlovic MD, Scott JA. HSL MI20: an efficient AMG preconditioner. Technical Report RAL-TR-2007-021, Rutherford Appleton Laboratory, Chilton, England, 2007.