SIAM Digital Library
 
 
 

You are not logged in Logged Out Log In

SIAM J. Sci. Comput. 28, pp. 1651-1674 (24 pages)

Recycling Krylov Subspaces for Sequences of Linear Systems

Michael L. Parks, Eric de Sturler, Greg Mackey, Duane D. Johnson, and Spandan Maiti

Full Text: Download PDF | Buy PDF (US$25) | View Cart
Many problems in science and engineering require the solution of a long sequence of slowly changing linear systems. We propose and analyze two methods that significantly reduce the total number of matrix‐vector products required to solve all systems. We consider the general case where both the matrix and right‐hand side change, and we make no assumptions regarding the change in the right‐hand sides. Furthermore, we consider general nonsingular matrices, and we do not assume that all matrices are pairwise close or that the sequence of matrices converges to a particular matrix. Our methods work well under these general assumptions, and hence form a significant advancement with respect to related work in this area. We can reduce the cost of solving subsequent systems in the sequence by recycling selected subspaces generated for previous systems. We consider two approaches that allow for the continuous improvement of the recycled subspace at low cost. We consider both Hermitian and non‐Hermitian problems, and we analyze our algorithms both theoretically and numerically to illustrate the effects of subspace recycling. We also demonstrate the effectiveness of our algorithms for a range of applications from computational mechanics, materials science, and computational physics.

© 2006 Society for Industrial and Applied Mathematics

RELATED DATABASES

To view database links for this article, you need to log in.

PUBLICATION DATA

ISSN

1064-8275 (print)  
1095-7197 (online)

ARTICLE DATA

History
Received April 24, 2004
Accepted February 16, 2006
Published online October 06, 2006

For access to fully linked references, you need to log in.

For access to citing articles, you need to log in.


Close

close