Abstract
To improve the performance of block triangular (BT) preconditioner, we develop a two-parameter BT (TPBT) preconditioner for a double saddle point problem arising from liquid crystal directors modeling. Theoretical analysis shows that all the eigenvalues of the TPBT preconditioned coefficient matrix are real and located in an interval (0, 2) no matter which value the spectral radius of matrix D− 1CA− 1CT is chosen. Moreover, an upper bound of the degree of the minimal polynomial of the TPBT preconditioned coefficient matrix is also obtained. Inasmuch as the efficiency of the TPBT preconditioner depends on the values of its parameters, we further derive a class of fast and effective formulas to compute the quasi-optimal values of the parameters involved in the TPBT preconditioner. Finally, numerical results are reported to illustrate the feasibility and the efficiency of the TPBT preconditioner.
Similar content being viewed by others
References
Stewart, I.W.: The static and dynamic continuum theory of liquid crystals: A mathematical introduction. Taylor and Francis, London (2004)
de Gennes, P.G., Prost, J.: The physics of liquid crystals, 2nd edn. Clarendon Press, Oxford (1993)
Ramage, A., Gartland Jr, E.C.: A preconditioned nullspace method for liquid crystal director modeling. SIAM J. Sci. Comput. 35(1), B226–B247 (2013)
Bramble, J.H., Pasciak, J.E., Vassilev, A.T.: Analysis of the inexact Uzawa algorithm for saddle point problems. SIAM J. Numer. Anal. 34 (3), 1072–1092 (1997)
Bai, Z.-Z., Wang, Z.-Q.: On parameterized inexact Uzawa methods for generalized saddle point problems. Linear Algebra Appl. 428(11-12), 2900–2932 (2008)
Golub, G.H., Wu, X., Yuan, J.-Y.: SOR-like methods for augmented systems. BIT Numer. Math. 41(1), 71–85 (2001)
Bai, Z.-Z., Parlett, B.N., Wang, Z.-Q.: On generalized successive overrelaxation methods for augmented linear systems. Numer. Math. 102(1), 1–38 (2005)
Guo, P., Li, C.-X., Wu, S.-L.: A modified SOR-like method for the augmented systems. J. Comput. Appl. Math. 274(1), 58–69 (2015)
Bai, Z.-Z., Golub, G.H., Ng, M.K.: Hermitian and skew-Hermitian splitting methods for non-Hermitian positive definite linear systems. SIAM J. Matrix Anal. Appl. 24, 603–626 (2003)
Benzi, M., Golub, G.H.: A preconditioner for generalized saddle point problems. SIAM J. Matrix Anal. Appl. 26(1), 20–41 (2005)
Bai, Z-Z, Golub, G.H.: Accelerated Hermitian and skew-Hermitian splitting iteration methods for saddle-point problems. IMA J. Numer. Anal. 27, 1–23 (2007)
Huang, T.-Z., Wu, S.-L., Li, C.-X.: The spectral properties of the Hermitian and skew-Hermitian splitting preconditioner for generalized saddle point problems. J. Comput. Appl. Math. 229(1), 37–46 (2009)
Jiang, M.-Q., Cao, Y.: On local Hermitian and skew-Hermitian splitting iteration methods for generalized saddle point problems. J. Comput. Appl. Math. 231(2), 973–982 (2009)
Bai, Z.-Z., Wang, Z.-Q.: Restrictive preconditioners for conjugate gradient methods for symmetric positive definite linear systems. J. Comput. Appl. Math. 187, 202–226 (2006)
Bai, Z.-Z.: Motivations and realizations of Krylov subspace methods for large sparse linear systems. J. Comput. Appl. Math. 283, 71–78 (2015)
Simoncini, V.: Krylov subspace methods for saddle point problems with indefinite preconditioning. SIAM J. Matrix Anal. Appl. 24, 368–391 (2002)
Saad, Y.: Iterative Methods for Sparse Linear Systems, 2nd edn. SIAM, Philadelphia, PA, USA (2003)
Dollar, H.S.: Constraint-style preconditioners for regularized saddle-point problems. SIAM J. Matrix Anal. Appl. 29, 672–684 (2006)
Bai, Z.-Z., Ng, M.K., Wang, Z.-Q.: Constraint preconditioners for symmetric indefinite matrices. SIAM J. Matrix Anal. Appl. 31, 410–433 (2009)
Cao, Z.-H.: A class of constraint preconditioners for nonsymmetric saddle-point matrices. Numer. Math. 103, 47–61 (2006)
Grigori, L., Niu, Q., Xu, Y.-X.: Stabilized dimensional factorization preconditioner for solving incompressible Navier-Stokes equations. Appl. Numer. Math. 146, 309–327 (2019)
Shen, S.-Q.: A note on PSS preconditioners for generalized saddle point problems. Appl. Math. Comput. 237, 723–729 (2014)
Shen, Q.-Q., Shi, Q.: Generalized shift-splitting preconditioners for nonsingular and singular generalized saddle point problems. Comput. Math. Appl. 72, 632–641 (2016)
Cao, Y., Li, S.: Block triangular preconditioners based on symmetric-triangular decomposition for generalized saddle point problems. Appl. Math. Comput. 358, 262–277 (2019)
Beik, F.P.A., Benzi, M.: Block preconditioners for saddle point systems arising from liquid crystal directors modeling. Calcolo 55(3), 29 (2018)
Beik, F.P.A., Benzi, M.: Iterative methods for double saddle point systems. SIAM J. Matrix Anal. Appl. 39(2), 902–921 (2018)
Liang, Z.-Z., Zhang, G.-F.: Alternating positive semidefinite splitting preconditioners for double saddle point problems. Calcolo 56(3), 26 (2019)
Benzi, M., Deparis, S., Grandperrin, G., Quarteroni, A.: Parameter estimates for the relaxed dimensional factorization preconditioner and application to hemodynamics. Comput. Methods Appl. Mech. Eng. 300, 129–145 (2016)
Horn, R.A., Johnson, C.R.: Matrix analysis. 2nd edn. Cambridge University Press, Cambridge (1994)
Zhu, J.-L., Yang, A.-L., Wu, Y.-J.: A parameterized deteriorated PSS preconditioner and its optimization for nonsymmetric saddle point problems. Comput. Math. Appl. 79, 1420–1434 (2020)
Acknowledgements
The authors are very much indebted to the referees for their constructive and valuable comments and suggestions which greatly improved the original manuscript of this paper.
Funding
This work is partially supported by the National Natural Science Foundation of China (Grant Nos. 11471150, 11401281) and the National Key Research and Development Program of China (Grant No. 2018YFC0406600).
Author information
Authors and Affiliations
Corresponding authors
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Zhu, JL., Wu, YJ. & Yang, AL. A two-parameter block triangular preconditioner for double saddle point problem arising from liquid crystal directors modeling. Numer Algor 89, 987–1006 (2022). https://doi.org/10.1007/s11075-021-01142-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11075-021-01142-5
Keywords
- Double saddle point problem
- TPBT preconditioner
- Eigenvalue analysis
- Preconditioned matrix
- Optimal parameter