TY - JOUR
T1 - Using partial spectral information for block diagonal preconditioning of saddle-point systems
AU - Ramage, Alison
AU - Ruiz, Daniel
AU - Sartenaer, Annick
AU - Tannier, Charlotte
N1 - Funding Information:
We would like to very much thank all anonymous referees for their useful comments, which led in particular to a stronger result and shortened proof in Theorem . We would also like to thank Iain Duff and Andrew Wathen for their valuable remarks, questions and comments, that enabled us to substantially improve this paper. This work was partially granted by the ANR-BARESAFE project, ANR-11-MONU-004, Programme Modèles Numériques 2011, supported by the French National Agency for Research. This paper presents research results of the Belgian Network DYSCO (Dynamical Systems, Control, and Optimization), funded by the Interuniversity Attraction Poles Programme initiated by the Belgian Science Policy Office. The scientific responsibility rests with its authors.
Publisher Copyright:
© 2021, The Author(s), under exclusive licence to Springer Science+Business Media, LLC part of Springer Nature.
PY - 2021/3
Y1 - 2021/3
N2 - Considering saddle-point systems of the Karush–Kuhn–Tucker (KKT) form, we propose approximations of the “ideal” block diagonal preconditioner based on the exact Schur complement proposed by Murphy et al. (SIAM J Sci Comput 21(6):1969–1972, 2000). We focus on the case where the (1,1) block is symmetric and positive definite, but with a few very small eigenvalues that possibly affect the convergence of Krylov subspace methods like Minres. Assuming that these eigenvalues and their associated eigenvectors are available, we first propose a Schur complement preconditioner based on this knowledge and establish lower and upper bounds on the preconditioned Schur complement. We next analyse theoretically the spectral properties of the preconditioned KKT systems using this Schur complement approximation in two spectral preconditioners of block diagonal forms. In addition, we derive a condensed “two in one” formulation of the proposed preconditioners in combination with a preliminary level of preconditioning on the KKT system. Finally, we illustrate on a PDE test case how, in the context of a geometric multigrid framework, it is possible to construct practical block preconditioners that help to improve on the convergence of Minres.
AB - Considering saddle-point systems of the Karush–Kuhn–Tucker (KKT) form, we propose approximations of the “ideal” block diagonal preconditioner based on the exact Schur complement proposed by Murphy et al. (SIAM J Sci Comput 21(6):1969–1972, 2000). We focus on the case where the (1,1) block is symmetric and positive definite, but with a few very small eigenvalues that possibly affect the convergence of Krylov subspace methods like Minres. Assuming that these eigenvalues and their associated eigenvectors are available, we first propose a Schur complement preconditioner based on this knowledge and establish lower and upper bounds on the preconditioned Schur complement. We next analyse theoretically the spectral properties of the preconditioned KKT systems using this Schur complement approximation in two spectral preconditioners of block diagonal forms. In addition, we derive a condensed “two in one” formulation of the proposed preconditioners in combination with a preliminary level of preconditioning on the KKT system. Finally, we illustrate on a PDE test case how, in the context of a geometric multigrid framework, it is possible to construct practical block preconditioners that help to improve on the convergence of Minres.
KW - Block diagonal preconditioning
KW - Preconditioned Krylov methods
KW - Saddle-point linear systems
KW - Spectral preconditioning
UR - http://www.scopus.com/inward/record.url?scp=85098876862&partnerID=8YFLogxK
U2 - 10.1007/s10589-020-00246-3
DO - 10.1007/s10589-020-00246-3
M3 - Article
AN - SCOPUS:85098876862
SN - 0926-6003
VL - 78
SP - 353
EP - 375
JO - Computational Optimization and Applications
JF - Computational Optimization and Applications
IS - 2
ER -