Abstract
AbstractNull-space methods have long been used to solve large sparse n × n symmetric saddle point systems of equations in which the (2, 2) block is zero. This paper focuses on the case where the (1, 1) block is ill conditioned or rank deficient and the k × k (2, 2) block is non zero and small (k ≪ n). Additionally, the (2, 1) block may be rank deficient. Such systems arise in a range of practical applications. A novel null-space approach is proposed that transforms the system matrix into a nicer symmetric saddle point matrix of order n that has a non zero (2, 2) block of order at most 2k and, importantly, the (1, 1) block is symmetric positive definite. Success of any null-space approach depends on constructing a suitable null-space basis. We propose methods for wide matrices having far fewer rows than columns with the aim of balancing stability of the transformed saddle point matrix with preserving sparsity in the (1, 1) block. Linear least squares problems that contain a small number of dense rows are an important motivation and are used to illustrate our ideas and to explore their potential for solving large-scale systems.
Publisher
Springer Science and Business Media LLC
Reference62 articles.
1. Benzi, M., Golub, G. H., Liesen, J.: Numerical solution of saddle point problems. Acta Numerica 14, 1–137 (2005)
2. Govaerts, W.: Solution of bordered singular systems in numerical continuation and bifurcation. In: Proceedings of the Fifth International Congress on Computational and Applied Mathematics (Leuven, 1992). https://doi.org/10.1016/0377-0427(94)90311-5, vol. 50, pp 339–347 (1994)
3. Howell, J. S.: Prestructuring sparse matrices with dense rows and columns via null space methods. Numer. Linear Algebra Appl. 25, 1–30 (2018). https://doi.org/10.1002/nla.2133
4. Benzi, M.: A direct row-projection method for sparse linear systems. PhD thesis, Department of Mathematics North Carolina State University (1993)
5. Benzi, M., Meyer, C. D.: A direct projection method for sparse linear systems. SIAM J. Sci. Comput. 16(5), 1159–1176 (1995)
Cited by
4 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献