On large-scale generalized inverses in solving two-by-two block linear systems (Q1940097): Difference between revisions
From MaRDI portal
Revision as of 05:40, 6 July 2024
scientific article
Language | Label | Description | Also known as |
---|---|---|---|
English | On large-scale generalized inverses in solving two-by-two block linear systems |
scientific article |
Statements
On large-scale generalized inverses in solving two-by-two block linear systems (English)
0 references
5 March 2013
0 references
In large scale nonsymmetric saddle-point problems, one has to solve a linear system whose matrix is of the form \[ \left[\begin{matrix} A&B_1^T\cr B_2&-C\end{matrix}\right] \] with \(A\) and \(C\) square. A projected Schur complement method (PSCM) is used (see [\textit{J.~Haslinger} et al., Numer.\ Linear Algebra Appl. 14, No.~9, 713--739 (2007; Zbl 1199.65102)]). This paper discusses the role of the generalized inverse in the Schur complement of \(A\) when \(A\) is singular. The key point is that a generalized inverse with prescribed range and null space (e.g.\ the Moore-Penrose inverse) can be obtained by choosing appropriate projectors. This is a generalization of \textit{R.~Kučera} et al. [``On the Moore-Penrose inverse in solving saddle-point systems with singular diagonal blocks'', ibid. 19, No.~4, 677--699 (2012)]. Since in practice not the generalized inverse but its action on a vector is needed, an algorithm is described to compute that. Here, also a regularization method is proposed in the case where the null space is known. It causes only little fill-in and avoids the need to check for zero pivots. Numerical examples illustrate the method proposed.
0 references
generalized inverses
0 references
projectors
0 references
block linear systems
0 references
Schur complement
0 references
null-space method
0 references
nonsymmetric saddle-point problem
0 references
Moore-Penrose inverse
0 references
algorithm
0 references
regularization method
0 references
numerical examples
0 references
0 references
0 references
0 references
0 references
0 references