Linear Iterative Solvers for Implicit Ode Methods
Abstract
In this paper we consider the numerical solution of stiff initial value problems, which lead to the problem of solving large systems of mildly nonlinear equations. For many problems derived from engineering and science, a solution is possible only with methods derived from iterative linear equation solvers. A common approach to solving the nonlinear equations is to employ an approximate solution obtained from an explicit method. In this paper we shall examine the error to determine how it is distributed among the stiff and non- stiff components, which bears on the choice of an iterative method. Our conclusion is that error is (roughly) uniformly distributed, a fact that suggests the Chebyshev method (and the accompanying Manteuffel adaptive parameter algorithm). We describe this method, also commenting on Richardson's method and its advantages for large problems. We then apply Richardson's method and the Chebyshev method with the Manteuffel algorithm to the solution of the nonlinear equations by Newton's method. (Author) (KR)
Document Details
- Document Type
- Technical Report
- Publication Date
- Aug 01, 1990
- Accession Number
- ADA227189
Entities
People
- Robert D. Skeel
- Robert E. Saylor