论文信息 - BCCB preconditioners for systems of BVM‐based numerical integrators

BCCB preconditioners for systems of BVM‐based numerical integrators

Boundary value methods (BVMs) for ordinary di erential equations require the solution of nonsymmetric, large and sparse linear systems. In this paper, these systems are solved by using the generalized minimal residual (GMRES) method. A block-circulant preconditioner with circulant blocks (BCCB preconditioner) is proposed to speed up the convergence rate of the GMRES method. The BCCB preconditioner is shown to be invertible when the BVM is Ak1 ; k2 -stable. The spectrum of the preconditioned matrix is clustered and therefore, the preconditioned GMRES method converges fast. Moreover, the operation cost in each iteration of the preconditioned GMRES method by using our BCCB preconditioner is less than that required by using block-circulant preconditioners proposed earlier. In numerical experiments, we compare the number of iterations of various preconditioners. Copyright ? 2003 John Wiley & Sons, Ltd.

Siu-Long Lei | Xiao-Qing Jin | S. Lei | X. Jin

[1] Y. Saad,et al. GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[2] Raymond H. Chan,et al. Circulant preconditioners for second order hyperbolic equations , 1992 .

[3] G. Strang. A proposal for toeplitz matrix calculations , 1986 .

[4] Daniele Bertaccini,et al. A Circulant Preconditioner for the Systems of LMF-Based ODE Codes , 2000, SIAM J. Sci. Comput..

[5] T. Chan. An Optimal Circulant Preconditioner for Toeplitz Systems , 1988 .

[6] Raymond H. Chan,et al. Strang‐type preconditioners for systems of LMF‐based ODE codes , 2001 .

[7] Siu-Long Lei,et al. Sine transform based preconditioners for solving constant-coefficient first-order PDEs☆ , 2003 .

[8] Daniele Bertaccini. Reliable preconditioned iterative linear solvers for some numerical integrators , 2001 .

[9] Daniele Bertaccini,et al. The Spectrum of Circulant-Like Preconditioners for Some General Linear Multistep Formulas for Linear Boundary Value Problems , 2002, SIAM J. Numer. Anal..

[10] J. G. Verwer,et al. Boundary value techniques for initial value problems in ordinary differential equations , 1983 .

[11] Henk A. van der Vorst,et al. Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..