The Convergence Rate of Block Preconditioned Systems Arising From LMF-Based ODE Codes

The solution of ordinary and partial differential equations using implicit linear multistep formulas (LMF) is considered. More precisely, boundary value methods (BVMs), a class of methods based on implicit formulas is taken into account in this paper. These methods require the solution of large and sparse linear systems Mˆx=b. Block-circulant preconditioners are proposed to solve these linear systems. By investigating the spectral condition number of Mˆ, we show that the conjugate gradient method, when applied to solving the normalized preconditioned system, converges in at most O(logs) steps, where the integration step size is O(1/s). Numerical results are given to illustrate the effectiveness of the analysis