Strathprints Home | Open Access | Browse | Search | User area | Copyright | Help | Library Home | SUPrimo

Analysis of the Enright-Kamel partitioning method for stiff ODEs

Higham, D.J. (1989) Analysis of the Enright-Kamel partitioning method for stiff ODEs. IMA Journal of Numerical Analysis, 9 (1). pp. 1-14. ISSN 0272-4979

Full text not available in this repository. (Request a copy from the Strathclyde author)

Abstract

The use of implicit formulae in the solution of stiff ODEs gives rise to systems of nonlinear equations which are usually solved iteratively by a modified Newton scheme. The linear algebra costs associated with such schemes may form a substantial part of the overall cost of the solution. The work of W. H. Enright and M. S. Kamel attempts to reduce the cost of the iteration by automatically transforming and partitioning the system. We provide new theoretical justification for this method in the case where the stiff eigenvalues of the Jacobian matrix used in the modified Newton iteration are small in number and well separated from the other eigenvalues. The theory of Y. Saad is introduced and adapted to show that the method uses the projection of the Jacobian onto a Krylov subspace which virtually contains the dominant subspace. This is shown to have favourable consequences. Numerical evidence is provided to support the theory.

Item type: Article
ID code: 214
Keywords: nonlinear equations, numerical mathematics, Newton scheme, Krylov subspace, Mathematics, Computational Mathematics, Applied Mathematics, Mathematics(all)
Subjects: Science > Mathematics
Department: Faculty of Science > Mathematics and Statistics
Related URLs:
    Depositing user: Ms Sarah Scott
    Date Deposited: 08 Mar 2006
    Last modified: 04 Sep 2014 09:52
    URI: http://strathprints.strath.ac.uk/id/eprint/214

    Actions (login required)

    View Item