Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

The use of iterative linear-equation solvers in codes for large systems of stiff IVPs for ODEs

Journal Article · · SIAM J. Sci. Stat. Comput.; (United States)
DOI:https://doi.org/10.1137/0907027· OSTI ID:5609802
Systems of linear algebraic equations must be solved at each integration step in all commonly used methods for the numerical solution of systems of stiff IVPs for ODEs. Frequently, a substantial portion of the total computational work and storage required to solve stiff IVPs is devoted to solving these linear algebraic systems, particularly if the systems are large. Over the past decade, several efficient iterative methods have been developed to solve large sparse (nonsymmetric) systems of linear algebraic equations. The authors study the use of a class of these iterative methods in codes for stiff IVPs. The theoretical estimates and preliminary numerical results show that the use of iterative linear-equation solvers in stiff-ODE codes improves the efficiency - in terms of both computational-work and storage - with which a significant class of stiff IVPs having large sparse Jacobians can be solved.
Research Organization:
Department of Computer Science, Yale Univ., New Haven, CT 06520
OSTI ID:
5609802
Journal Information:
SIAM J. Sci. Stat. Comput.; (United States), Journal Name: SIAM J. Sci. Stat. Comput.; (United States) Vol. 7:2; ISSN SIJCD
Country of Publication:
United States
Language:
English