Applications of Mathematics, Vol. 62, No. 6, pp. 537-559, 2017


Preconditioning of two-by-two block matrix systems with square matrix blocks, with applications

Owe Axelsson

Received August 16, 2017.   First published December 4, 2017.

Abstract:  Two-by-two block matrices of special form with square matrix blocks arise in important applications, such as in optimal control of partial differential equations and in high order time integration methods. Two solution methods involving very efficient preconditioned matrices, one based on a Schur complement reduction of the given system and one based on a transformation matrix with a perturbation of one of the given matrix blocks are presented. The first method involves an additional inner solution with the pivot matrix block but gives a very tight condition number bound when applied for a time integration method. The second method does not involve this matrix block but only inner solutions with a linear combination of the pivot block and the off-diagonal matrix blocks. Both the methods give small condition number bounds that hold uniformly in all parameters involved in the problem, i.e. are fully robust. The paper presents shorter proofs, extended and new results compared to earlier publications.
Keywords:  preconditioning; Schur complement; transformation; optimal control; implicit time integration
Classification MSC:  65F08


References:
[1] O. Axelsson: On the efficiency of a class of $A$-stable methods. BIT, Nord. Tidskr. Inf.-behandl. 14 (1974), 279-287. DOI 10.1007/BF01933227 | MR 0391518 | Zbl 0289.65028
[2] O. Axelsson, R. Blaheta, R. Kohut: Preconditioning methods for high-order strongly stable time integration methods with an application for a DAE problem. Numer. Linear Algebra Appl. 22 (2015), 930-949. DOI 10.1002/nla.2015 | MR 3426322 | Zbl 06604516
[3] O. Axelsson, S. Farouq, M. Neytcheva: Comparison of preconditioned Krylov subspace iteration methods for PDE-constrained optimization problems. Poisson and convection-diffusion control. Numer. Algorithms 73 (2016), 631-663. DOI 10.1007/s11075-016-0111-1 | MR 3564863 | Zbl 1353.65059
[4] O. Axelsson, S. Farouq, M. Neytcheva: Comparison of preconditioned Krylov subspace iteration methods for PDE-constrained optimization problems. Stokes control. Numer. Algorithms 74 (2017), 19-37. DOI 10.1007/s11075-016-0136-5 | MR 3590387 | Zbl 1365.65167
[5] O. Axelsson, A. Kucherov: Real valued iterative methods for solving complex symmetric linear systems. Numer. Linear Algebra Appl. 7 (2000), 197-218. DOI 10.1002/1099-1506(200005)7:4<197::AID-NLA194>3.0.CO;2-S | MR 1762967 | Zbl 1051.65025
[6] O. Axelsson, M. Neytcheva, B. Ahmad: A comparison of iterative methods to solve complex valued linear algebraic systems. Numer. Algorithms 66 (2014), 811-841. DOI 10.1007/s11075-013-9764-1 | MR 3240302 | Zbl 1307.65034
[7] O. Axelsson, P. S. Vassilevski: A black box generalized conjugate gradient solver with inner iterations and variable-step preconditioning. SIAM J. Matrix Anal. Appl. 12 (1991), 625-644. DOI 10.1137/0612048 | MR 1121697 | Zbl 0748.65028
[8] Z.-Z. Bai: On preconditioned iteration methods for complex linear systems. J. Eng. Math. 93 (2015), 41-60. DOI 10.1007/s10665-013-9670-5 | MR 3386091 | Zbl 1360.65089
[9] Z.-Z. Bai, M. Benzi, F. Chen, Z.-Q. Wang: Preconditioned MHSS iteration methods for a class of block two-by-two linear systems with applications to distributed control problems. IMA J. Numer. Anal. 33 (2013), 343-369. DOI 10.1093/imanum/drs001 | MR 3020961 | Zbl 1271.65100
[10] Z.-Z. Bai, F. Chen, Z.-Q. Wang: Additive block diagonal preconditioning for block two-by-two linear systems of skew-Hamiltonian coefficient matrices. Numer. Algorithms 62 (2013), 655-675. DOI 10.1007/s11075-013-9696-9 | MR 3034831 | Zbl 1267.65034
[11] F. Brezzi, M. Fortin: Mixed and Hybrid Finite Element Methods. Springer Series in Computational Mathematics 15, Springer, New York (1991). DOI 10.1007/978-1-4612-3172-1 | MR 1115205 | Zbl 0788.73002
[12] J. C. Butcher: Numerical Methods for Ordinary Differential Equations. John Wiley & Sons, Chichester (2008). DOI 10.1002/9780470753767 | MR 2401398 | Zbl 1167.65041
[13] J. Cahouet, J.-P. Chabard: Some fast 3D finite element solvers for the generalized Stokes problem. Int. J. Numer. Methods Fluids 8 (1988), 869-895. DOI 10.1002/fld.1650080802 | MR 0953141 | Zbl 0665.76038
[14] A. Greenbaum, V. Pták, Z. Strakoš: Any nonincreasing convergence curve is possible for GMRES. SIAM J. Matrix Anal. Appl. 17 (1996), 465-469. DOI 10.1137/S0895479894275030 | MR 1397238 | Zbl 0857.65029
[15] J. L. Lions: Some Aspects of the Optimal Control of Distributed Parameter Systems. CBMS-NSF Regional Conference Series in Applied Mathematics 6, Society for Industrial and Applied Mathematics, Philadelphia (1972). DOI 10.1137/1.9781611970616 | MR 0479526 | Zbl 0275.49001
[16] C. C. Paige, M. A. Saunders: Solution of sparse indefinite systems of linear equations. SIAM J. Numer. Anal. 12 (1975), 617-629. DOI 10.1137/0712047 | MR 0383715 | Zbl 0319.65025
[17] Y. Saad: A flexible inner-outer preconditioned GMRES-algorithm. SIAM J. Sci. Comp. 14 (1993), 461-469. DOI 10.1137/0914028 | MR 1204241 | Zbl 0780.65022

Affiliations:   Owe Axelsson, Institute of Geonics of the Czech Academy of Sciences, Studentská 1768, 708 00 Ostrava, Czech Republic, e-mail: owe.axelsson@it.uu.se


 
PDF available at: