Enlarged GMRES for solving linear systems with one or multiple right-hand sides - Inria - Institut national de recherche en sciences et technologies du numérique Accéder directement au contenu
Article Dans Une Revue IMA Journal of Numerical Analysis Année : 2018

Enlarged GMRES for solving linear systems with one or multiple right-hand sides

Résumé

We propose a variant of the generalized minimal residual (GMRES) method for solving linear systems of equations with one or multiple right-hand sides. Our method is based on the idea of the enlarged Krylov subspace to reduce communication. It can be interpreted as a block GMRES method. Hence, we are interested in detecting inexact breakdowns. We introduce a strategy to perform the test of detection. Furthermore, we propose a technique for deflating eigenvalues that has two benefits. The first advantage is to avoid the plateau of convergence after the end of a cycle in the restarted version. The second is to have very fast convergence when solving the same system with different right-hand sides, each given at a different time (useful in the context of a constrained pressure residual preconditioner). We test our method with these deflation techniques on academic test matrices arising from solving linear elasticity and convection–diffusion problems as well as matrices arising from two real-life applications, seismic imaging and simulations of reservoirs. With the same memory cost we obtain a saving of up to 50% in the number of iterations required to reach convergence with respect to the original method.
Fichier non déposé

Dates et versions

hal-01963032 , version 1 (21-12-2018)

Identifiants

Citer

Hussam Al Daas, Laura Grigori, Pascal Hénon, Philippe Ricoux. Enlarged GMRES for solving linear systems with one or multiple right-hand sides. IMA Journal of Numerical Analysis, 2018, ⟨10.1093/imanum/dry054⟩. ⟨hal-01963032⟩
93 Consultations
0 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More