OATAO - Open Archive Toulouse Archive Ouverte Open Access Week

Limited memory preconditioners for symmetric indefinite problems with application to structural mechanics

Gratton, Serge and Mercier, Sylvain and Tardieu, Nicolas and Vasseur, Xavier Limited memory preconditioners for symmetric indefinite problems with application to structural mechanics. (2016) Numerical Linear Algebra with Applications, 23 (5). 865-887. ISSN 1070-5325

(Document in English)

PDF (Author's version) - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader

Official URL: http://dx.doi.org/10.1002/nla.2058


This paper presents a class of limited memory preconditioners (LMP) for solving linear systems of equations with symmetric indefinite matrices and multiple right-hand sides. These preconditioners based on limited memory quasi-Newton formulas require a small number k of linearly independent vectors and may be used to improve an existing first-level preconditioner. The contributions of the paper are threefold. First, we derive a formula to characterize the spectrum of the preconditioned operator. A spectral analysis of the preconditioned matrix shows that the eigenvalues are all real and that the LMP class is able to cluster at least k eigenvalues at 1. Secondly, we show that the eigenvalues of the preconditioned matrix enjoy interlacing properties with respect to the eigenvalues of the original matrix provided that the k linearly independent vectors have been prior projected onto the invariant subspaces associated with the eigenvalues of the original matrix in the open right and left half-plane, respectively. Third, we focus on theoretical properties of the Ritz-LMP variant, where Ritz information is used to determine the k vectors. Finally, we illustrate the numerical behaviour of the Ritz limited memory preconditioners on realistic applications in structural mechanics that require the solution of sequences of large-scale symmetric saddle-point systems. Numerical experiments show the relevance of the proposed preconditioner leading to a significant decrease in terms of computational operations when solving such sequences of linear systems. A saving of up to 43% in terms of computational effort is obtained on one of these applications.

Item Type:Article
Audience (journal):International peer-reviewed journal
Uncontrolled Keywords:
Institution:French research institutions > Commissariat à l'Energie Atomique et aux énergies alternatives - CEA (FRANCE)
French research institutions > Centre National de la Recherche Scientifique - CNRS (FRANCE)
Other partners > EDF (FRANCE)
Other partners > Ecole Nationale Supérieure de Techniques Avancées - ENSTA (FRANCE)
Université de Toulouse > Institut National Polytechnique de Toulouse - Toulouse INP (FRANCE)
Université de Toulouse > Université Toulouse III - Paul Sabatier - UT3 (FRANCE)
Université de Toulouse > Université Toulouse - Jean Jaurès - UT2J (FRANCE)
Université de Toulouse > Université Toulouse 1 Capitole - UT1 (FRANCE)
Laboratory name:
Deposited On:04 Nov 2016 17:15

Repository Staff Only: item control page