Abstract
When the Newton method or the fixed-point method is employed to solve the systems of nonlinear equations arising in the sinc-Galerkin discretization of certain time-dependent partial differential equations, in each iteration step we need to solve a structured subsystem of linear equations iteratively by, for example, a Krylov subspace method such as the preconditioned GMRES. In this paper, based on the tensor and the Toeplitz structures of the linear subsystems we construct structured preconditioners for their coefficient matrices and estimate the eigenvalue bounds of the preconditioned matrices under certain assumptions. Numerical examples are given to illustrate the effectiveness of the proposed preconditioning methods. It has been shown that a combination of the Newton/fixed-point iteration with the preconditioned GMRES method is efficient and robust for solving the systems of nonlinear equations arising from the sinc-Galerkin discretization of the time-dependent partial differential equations.
Original language | English |
---|---|
Pages (from-to) | 1019-1037 |
Number of pages | 19 |
Journal | SIAM Journal on Numerical Analysis |
Volume | 47 |
Issue number | 2 |
DOIs | |
Publication status | Published - 13 Feb 2009 |
Scopus Subject Areas
- Numerical Analysis
- Computational Mathematics
- Applied Mathematics
User-Defined Keywords
- Eigenvalue bound
- GMRES method
- Preconditioning
- Sinc-Galerkin discretization
- Time-dependent partial differential equation
- Toeplitzlike matrix