Strang-type preconditioners for systems of LMF-based ODE codes

Raymond H. CHAN*, Michael K. NG, Xiao Qing JIN

*Corresponding author for this work

Research output: Journal PublicationsJournal Article (refereed)peer-review

38 Citations (Scopus)

Abstract

We consider the solution of ordinary differential equations (ODEs) using boundary value methods. These methods require the solution of one or more unsymmetric, large and sparse linear systems. The GMRES method with the Strang-type block-circulant preconditioner is proposed for solving these linear systems. We show that if an Ak1,k2-stable boundary value method is used for an m-by-m system of ODEs, then our preconditioners are invertible and all the eigenvalues of the preconditioned systems are 1 except for at most 2m(k1 + k2) outliers. It follows that when the GMRES method is applied to solving the preconditioned systems, the method will converge in at most 2m(k1 + k2) + 1 iterations. Numerical results are given to illustrate the effectiveness of our methods.

Original languageEnglish
Pages (from-to)451-462
Number of pages12
JournalIMA Journal of Numerical Analysis
Volume21
Issue number2
DOIs
Publication statusPublished - 1 Apr 2001
Externally publishedYes

Fingerprint

Dive into the research topics of 'Strang-type preconditioners for systems of LMF-based ODE codes'. Together they form a unique fingerprint.

Cite this