The saddle point problem and the Manteuffel algorithm

详细信息    查看全文

• 作者：Apostolos Hadjidimos
• 关键词：Singular saddle point problem ; Generalized SOR iterative methods ; Generalized MSSOR iterative methods ; Second ; order iterative method ; Manteuffel’s algorithm
• 刊名：BIT Numerical Mathematics
• 出版年：2016
• 出版时间：December 2016
• 年：2016
• 卷：56
• 期：4
• 页码：1281-1302
• 全文大小：
• 刊物类别：Mathematics and Statistics
• 刊物主题：Computational Mathematics and Numerical Analysis; Numeric Computing; Mathematics, general;
• 出版者：Springer Netherlands
• ISSN：1572-9125
• 卷排序：56

文摘

In the last two decades, the augmented linear systems and the saddle point problems have been solved by many researchers who have used the conjugate gradient method or the generalized SOR iterative method and variants of them. In the latter class of methods, when the block \(A \in {\mathrm{I\!R}}^{m\times m}\) of the matrix coefficient \(\mathcal {A} = \left[ \begin{array}{cc} A &{} B \\ -B^T &{} 0\end{array}\right] \in {\mathrm{I\!R}}^{(m+n) \times (m+n)}\), \(m \ge n,\) of the linear system to be solved, is symmetric positive definite and \({\mathrm{rank}}(B) =r \le n\), convergence regions and optimal values of the parameters involved have been determined. In this work, we consider the block A to be nonsymmetric positive definite, \({\mathrm{rank}}(B) =r < n \), and use a two-level stationary iterative method whose main step is the linear second-order stationary iterative method for the solution of this class of problems. This method leads to the singular Manteuffel algorithm and the determination of its optimal parameters. As a byproduct, the optimal parameters of the Generalized Modified SSOR method in a particular case are also determined. Numerical examples verify our theoretical findings.