Interpolating discrete advection-diffusion propagators at Leja sequences

Journal of Computational and Applied Mathematics

Volumn 172, Issue 1, 2004, Pages 79-99

  Caliari, M ^a   Vianello, M ^b   Bergamaschi, L ^b  

^a UNIVERSITY OF VERONA   (Italy)

^b UNIVERSITY OF PADOVA   (Italy)

Author keywords

Advection diffusion problem; Exponential operator; Leja sequence; Polynomial interpolation; Sparse matrix

Indexed keywords

DIFFUSION; FINITE DIFFERENCE METHOD; INTERPOLATION; MATRIX ALGEBRA; SPECTRUM ANALYSIS;

ADVECTION; MATRIX INTERPOLATION;

COMPUTATIONAL METHODS;

ADVECTION-DIFFUSION EQUATION; MATHEMATICAL METHOD;

EID: 4444264388     PISSN: 03770427     EISSN: None     Source Type: Journal    
DOI: 10.1016/j.cam.2003.11.015     Document Type: Article

References (42)

1. Baglama J. Calvetti D. Reichel L. Fast Leja points Electron. Trans. Numer. Anal. 7 1998 124-140
2. Bergamaschi L. Caliari M. Vianello M. Efficient approximation of the exponential operator for discrete 2D advection-diffusion problems Numer. Linear Algebra Appl. 10 3 2003 271-289
3. Bergamaschi L. Vianello M. Efficient computation of the exponential operator for large, sparse, symmetric matrices Numer. Linear Algebra Appl. 7 1 2000 27-45
4. M. Caliari, Efficient implementation of exponential integrators for 2D and 3D advection-diffusion equations, Ph.D. Thesis in Computational Mathematics, University of Padova, advisors M. Vianello and L. Bergamaschi, 2003, URL http://www.math.unipd.it/~mcaliari/pdf/phd.pdf http://www.math.unipd.it/~mcaliari/pdf/phd.pdf
5. Crank J. Nicolson P. A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type Proc. Cambridge Philos. Soc. 43 1947 50-67
6. Druskin V.L. Knizhnerman L.A. Two polynomial methods for calculating functions of symmetric matrices Comput. Math. Math. Phys. 29 6 1989 112-121
7. Ellacott S.W. Computation of Faber series with application to numerical polynomial approximation in the complex plane Math. Comput. 40 162 1983 575-587
8. Fischer B. Reichel L. Newton interpolation in Fejér and Chebyshev points Math. Comput. 53 187 1989 265-278
9. Friesner R.A. Tuckerman L.S. Dornblaser B.C. Russo T.V. A method for exponential propagation of large system of stiff nonlinear differential equations J. Sci. Comput. 4 4 1989 327-354
10. Gallopoulos E. Saad Y. Efficient solution of parabolic equations by Krylov subspace methods SIAM J. Sci. Statist. Comput. 13 5 1992 1236-1264
11. Hochbruck M. Lubich C. On Krylov subspace approximations to the matrix exponential SIAM J. Numer. Anal. 34 5 1997 1911-1925
12. Hochbruck M. Lubich C. Selhofer H. Exponential integrators for large systems of differential equations SIAM J. Sci. Comput. 19 5 1998 1552-1574
13. W. Hundsdorfer, J.G. Verwer, Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations, Springer series in Computational Mathematics, Vol. 33, Springer, Berlin, 2003.
14. Kershaw D.S. The incomplete Cholesky-conjugate gradient method for the iterative solution of systems of linear equations J. Comput. Phys. 26 1 1978 43-65
15. Lanser D. Verwer J.G. Analysis of operator splitting for advection-diffusion-reaction problems for air-pollution modelling J. Comput. Appl. Math. 111 1-2 1999 201-216
16. Le Calvez C. Saad Y. Modified Krylov acceleration for parallel environments Appl. Numer. Math. 30 2-3 1999 191-212
17. R.B Lehoucq, D.C. Sorensen, C. Yang, ARPACK Users' Guide: Solution of Large Scale Eigenvalue Problems by Implicity Restarted Arnoldi Methods, 1997, URL http://www.caam.rice.edu/software/ARPACK/UG/ug.html http://www.caam.rice.edu/software/ARPACK/UG/ug.html
18. Leja F. Sur certaines suites liées aux ensembles plans et leur application à la représentation conforme Ann. Polon. Math. 4 1957 8-13
19. Manteuffel T.A. The Tchebychev iteration for nonsymmetric linear systems Numer. Math. 28 3 1977 307-327
20. Manteuffel T.A. Adaptive procedure for estimating parameters for the nonsymmetric Tchebychev iteration Numer. Math. 31 2 1978/1979 183-208
21. A.I. Markushevic, Theory of Functions of a Complex Variable, Vol. 3, Prentice-Hall, Inc., Englewood Cliffs, NJ, 1967.
22. Moler C. van Loan C. Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later SIAM Rev. 45 2003 3-49
23. Moret I. Novati P. The computation of functions of matrices by truncated Faber series Numer. Funct. Anal. Optim. 22 5-6 2001 697-719
24. Moret I. Novati P. An interpolatory approximation of the matrix exponential based on Faber polynomials J. Comput. Appl. Math. 131 1-2 2001 361-380
25. I. Moret, P. Novati, Rational approximations of the matrix exponential. To appear in BIT (2003).
26. P. Novati, Polynomial methods for the computation of functions of large unsymmetric matrices, Ph.D. Thesis in Computational Mathematics, University of Trieste, advisor I. Moret, 2000.
27. Novati P. A polynomial method based on Fejér points for the computation of functions of unsymmetric matrices Appl. Numer. Math. 44 1-2 2003 201-224
28. Novati P. Solving linear initial value problems by Faber polynomials Numer. Linear Algebra Appl. 10 3 2003 247-270
29. Pini G. Gambolati G. Arnoldi and Crank-Nicolson methods for integration in time of the transport equation Internat. J. Numer. Methods Fluids 35 2001 25-38
30. A. Quarteroni, A. Valli, Numerical Approximation of Partial Differential Equations, Springer Series in Computational Mathematics, Vol. 23, Springer, Berlin, 1994.
31. Reddy S.C. Trefethen L.N. Pseudospectra of the convection-diffusion operator SIAM J. Appl. Math. 54 6 1994 1634-1649
32. Reichel L. Newton interpolation at Leja points BIT 30 2 1990 332-346
33. Reichel L. Trefethen L.N. Eigenvalues and pseudo-eigenvalues of Toeplitz matrices Linear Algebra Appl. 162/164 1992 153-185
34. Saad Y. Analysis of some Krylov subspace approximations to the matrix exponential operator SIAM J. Numer. Anal. 29 1 1992 209-228
35. Schaefer M.J. A polynomial based iterative method for linear parabolic equations J. Comput. Appl. Math. 29 1 1990 35-50
36. Sidje R.B. Expokit, A software package for computing matrix exponentials ACM Trans. Math. Software 24 1 1998 130-156
37. Tal-Ezer H. Kosloff R. An accurate and efficient scheme for propagating the time dependend Schrödinger equation J. Chem. Phys. 81 9 1984 3967-3971
38. Tal-Ezer H. Kosloff R. Cerien C. Low-order polynomial approximation of propagators for the time-dependent Schrödinger equation J. Comput. Phys. 100 1 1992 179-187
39. L.N. Trefethen, Computation of pseudospectra, in: A. Iserles (Ed.), Acta Numerical, Vol. 8, Cambridge University Press, Cambridge, 1999, pp. 247-295.
40. van der Vorst H.A. Bi-CGSTAB: a fast and smoothly converging variant of Bi-CG for the solution of nonsymmetric linear systems SIAM J. Sci. Statist. Comput. 13 2 1992 631-644
41. C.B. Vreugdenhil, B. Koren (Eds.), Numerical methods for advection-diffusion problems, Notes on Numerical Fluid Mechanics, Vol. 45, Friedrivk, Vieweg and Sohn, Braunschweig, 1993.
42. J.L. Walsh, Interpolation and approximation by rational functions in the complex domain, American Mathematical Society Colloquium Publications, Vol. XX, AMS, Providence, RI, 1935.