Analysis of projection methods for rational function approximation to the matrix exponential

SIAM Journal on Numerical Analysis

Volumn 44, Issue 2, 2006, Pages 613-635

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

^a UNIVERSITY OF BARI   (Italy)

^b UNIVERSITY OF BOLOGNA   (Italy)

^c CIRSA   (Italy)

Author keywords

Exponential; Iterative methods; Krylov subspaces; Matrix; Rational functions

Indexed keywords

CONVERGENCE OF NUMERICAL METHODS; ERROR ANALYSIS; ITERATIVE METHODS; MATRIX ALGEBRA; POLYNOMIALS;

ERROR BOUNDS; EXPONENTIALS; KRYLOV SUBSPACES; RATIONAL FUNCTIONS;

APPROXIMATION THEORY;

EID: 33845347760     PISSN: 00361429     EISSN: None     Source Type: Journal    
DOI: 10.1137/05062590     Document Type: Article

References (53)

1. G. A. BAKER AND P. GRAVES-MORRIS, Padé Approximants, Encyclopedia of Mathematics and its Applications, Cambridge University Press, Cambridge, 1996.
2. C. BALDWIN, R. FREUND, AND E. GALLOPOULOS, A parallel iterative method for exponential propagation, in Proceedings of the Seventh SIAM Conference on Parallel Processing for Scientific Computing, D. Bailey et al., eds., SIAM, Philadelphia, 1995, pp. 534-539.
3. L. BERGAMASCHI AND M. VIANELLO, Efficient computation of the exponential operator for large, sparse, symmetric matrices, Numer. Linear Algebra Appl., 7 (2000), pp. 27-45.
4. A. BORICI, Computational methods for UV suppressed fermions, J. Comput. Phys., 189 (2003), pp. 454-462.
5. D. CALVETTI, E. GALLOPOULOS, AND L. REICHEL, Incomplete partial fractions for parallel evaluation of rational matrix functions, J. Comput. Appl. Math., 59 (1995), pp. 349-380.
6. P. CASTILLO AND Y. SAAD, Preconditioning the matrix exponential operator with applications, J. Sci. Comput., 13 (1999), pp. 225-302.
7. E. CELLEDONI AND I. MORET, A Krylov projection method for systems of ODEs, Appl. Numer. Math., 24 (1997), pp. 365-378.
8. -x in [0, +∞) and applications to heat-conduction problems, J. Approx. Theory, 2 (March 1969), pp. 50-65.
9. N. DEL BUONO, L. LOPEZ, AND R. PELUSO, Computation of exponentials of large sparse skew symmetric matrices, SIAM J. Sci. Comput., 27 (2005), pp. 278-293.
10. V. DRUSKIN, A. GREENBAUM, AND L. KNIZHNERMAN, Using nonorthogonal Lanczos vectors in the computation of matrix functions, SIAM J. Sci. Comput., 19 (1998), pp. 38-54.
11. V. DRUSKIN AND L. KNIZHNERMAN, Two polynomial methods of calculating functions of symmetric matrices, USSR Comput. Math. Math. Phys., 29 (1989), pp. 112-121.
12. V. DRUSKIN AND L. KNIZHNERMAN, Krylov subspace approximation of eigenpairs and matrix functions in exact and computer arithmetic, Numer. Linear Algebra Appl., 2 (1995), pp. 205-217.
13. V. DRUSKIN AND L. KNIZHNERMAN, Extended Krylov subspaces: Approximation of the matrix square root and related functions, SIAM J. Matrix Anal. Appl., 19 (1998), pp. 755-771.
14. M. EIERMANN AND O. ERNST, A Restarted Krylov Subspace Method for the Evaluation of Matrix Functions, technical report, Institut für Numerische Mathematik und Optimierung, Technische Universität Bergakademie Freiberg, Germany, 2005.
15. R. W. FREUND, On conjugate gradient type methods and polynomial preconditioners for a class of complex non-Hermitian matrices, Numer. Math., 57 (1990), pp. 285-312.
16. R. W. FREUND, Conjugate gradient-type methods for linear systems with complex symmetric coefficient matrices, SIAM J. Sci. Stat. Comput., 13 (1992), pp. 425-448.
17. R. FRIESNER, L. TUCKERMAN, B. DORNBLASER, AND T. RUSSO, A method for exponential propagation of large systems of stiff nonlinear differential equations, J. Sci. Comput., 4 (1989), pp. 327-354.
18. E. GALLOPOULOS AND Y. SAAD, A parallel block cyclic reduction algorithm for the fast solution of elliptic equations, Parallel Comput., 10 (1989), pp. 143-159.
19. E. GALLOPOULOS AND Y. SAAD, Efficient solution of parabolic equations by polynomial approximation methods, SIAM J. Sci. Stat. Comput., 13 (1992), pp. 1236-1264.
20. W. GAUTSCHI, Numerical Analysis. An Introduction, Birkhäuser, Boston, 1997.
21. G. GOLUB AND C. F. VAN LOAN, Matrix Computations, The Johns Hopkins University Press, 3rd ed., Baltimore, MD, 1996.
22. A. GREENBAUM, Using the Cauchy Integral Formula and the Partial fractions Decomposition of the Resolvent to Estimate ||f(A)||, technical report, Mathematics Department, University of Washington, Seattle, WA, 2000.
23. N. J. HIGHAM, The scaling and squaring method for the matrix exponential revisited, SIAM J. Matrix Anal. Appl., 26 (2005), pp. 1179-1193.
24. M. HOCHBRUCK AND C. LUBICH, On Krylov subspace approximations to the matrix exponential operator, SIAM J. Numer. Anal., 34 (1997), pp. 1911-1925.
25. M. HOCHBRUCK AND C. LUBICH, Exponential integrators for quantum-classical molecular dynamics, BIT, 39 (1999), pp. 620-645.
26. M. HOCHBRUCK, C. LUBICH, AND H. SELHOFER, Exponential integrators for large systems of differential equations, SIAM J. Sci. Comput., 19 (1998), pp. 1552-1574.
27. M. HOCHBRUCK AND J. NIEHOFF, private communication, 2005.
28. M. HUHTANEN, A Hermitian Lanczos method for normal matrices, SIAM J. Matrix Anal. Appl., 23 (2002), pp. 1092-1108.
29. L. KNIZHNERMAN, On Adaptation of the Lanczos Method to the Spectrum, Technical report #EMG-001-95-12, Schlumberger-Doll Research, Ridgefield, CT, 1995.
30. L. LOPEZ AND V. SIMONCINI, Preserving geometric properties of the exponential matrix by block Krylov subspace methods, Technical report, Dipartimento di Matematica, Università di Bologna, December, 2005.
31. Y. Y. LU, Computing a matrix function for exponential integrators, J. Comput. Appl. Math., 161 (2003), pp. 203-216.
32. THE MATH WORKS, INC., MATLAB 7, September 2004.
33. C. MOLER AND C. VAN LOAN, Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later, SIAM Rev., 45 (2003), pp. 3-49.
34. I. MORET AND P. NOVATI, RD-rational approximations of the matrix exponential, BIT, 44 (2004), pp. 595-615.
35. I. MORET AND P. NOVATI, Interpolating functions of matrices on zeros of quasi-kernel polynomials, Numer. Linear Algebra Appl., 11 (2005), pp. 337-353.
36. A. NAUTS AND R. WYATT, New approach to many state quantum dynamics: The recurisive residue generation method, Phys. Rev. Lett., 51 (1983), pp. 2238-2241.
37. S. P. NøRSETT, Restricted Fade approximations to the exponential function, SIAM J. Numer. Anal., 15 (Oct. 1978), pp. 1008-1029.
38. S. P. NøRSETT AND A. WOLFBRANDT, Attainable order of rational approximation to the exponential function with only real poles, BIT, 17 (1977), pp. 200-208.
39. T. PARK AND J. LIGHT, Unitary quantum time evolution by iterative Lanczos reduction, J. Chem. Phys., 85 (1986), pp. 5870-5876.
40. M. POPOLIZIO AND V. SIMONCINI, Acceleration techniques for the approximation of the matrix exponential with Krylov subspace methods, in preparation.
41. Y. SAAD, Analysis of some Krylov subspace approximations to the matrix exponential operator, SIAM J. Numer. Anal., 29 (1992), pp. 209-228.
42. Y. SAAD, Iterative Methods for Sparse Linear Systems, 2nd ed., SIAM, 2003.
43. R. B. SIDJE, Expokit: A software package for computing matrix exponentials, ACM Trans. Math. Software, 24 (1998), pp. 130-156.
44. V. SIMONCINI, On the convergence of restarted Krylov subspace methods, SIAM J. Matrix Anal. Appl., 22 (2000), pp. 430-452.
45. V. SIMONCINI, Restarted full orthogonalization method for shifted linear systems, BIT, 43 (2003), pp. 459-466.
46. V. SIMONCINI AND D. B. SZYLD, Theory of inexact Krylov subspace methods and applications to scientific computing, SIAM J. Sci. Comput., 25 (2003), pp. 454-477.
47. D. E. STEWART AND T. S. LEYK, Error estimates for Krylov subspace approximations of matrix exponentials, J. Comput. Appl. Math., 72 (1996), pp. 359-369.
48. H. TAL-EZER, Spectral methods in time for parabolic problems, SIAM J. Numer. Anal., 26 (1989), pp. 1-11.
49. J. VAN DEN ESHOF AND M. HOCHBRUCK, Preconditioning Lanczos Approximations to the Matrix Exponential, Technical report, Mathematical Institute, Heinrich Heine Universität, Düsseldorf, Germany, 2004. To appear in SIAM J. Sci. Comput.
50. H. VAN DER VORST, Iterative Krylov Methods for Large Linear Systems, Cambridge University Press, Cambridge, 2003.
51. H. A. VAN DER VORST, An iterative solution method for solving f(A)x=b, using Krylov subspace information obtained for the symmetric positive definite matrix A, J. Comput. Appl. Math., 18 (1987), pp. 249-263.
52. R. S. VARGA, On higher order stable implicit methods for solving parabolic partial differential equations, J. Math. Phys., XL (1961), pp. 220-231.
53. G. WALZ, Computing the matrix exponential and other matrix functions, J. Comput. Appl. Math., 21 (1988), pp. 119-123.