IIT Home Page CNR Home Page

Adaptive solution of infinite linear systems by Krylov subspace methods

In this paper we consider the problem of approximating the solution of infinite linear systems, finitely expressed by a sparse coefficient matrix. We analyse an algorithm based on Krylov subspace methods embedded in an adaptive enlargement scheme. The management of the algorithm is not trivial, due to the irregular convergence behaviour frequently displayed by Krylov subspace methods for nonsymmetric systems. Numerical experiments, carried out on several test problems, indicate that the more robust methods, such as GMRES and QMR, embedded in the adaptive enlargement scheme, exhibit good performances.


JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS (05795J0), 2007

Authors: P. Favati, G. Lotti, O. Menchi, F. Romani
IIT authors:

Type: Article in ISI Journal
Field of reference: Mathematics
Da pagina 191 a pagina 199

Activity: Metodi numerici per problemi di grandi dimensioni