Numerical validation of solutions of saddle point matrix equations

From MaRDI portal





A method is presented for estimating the accuracy of approximate solutions to linear systems defined by saddle point matrices. A saddle point matrix comprises four blocks: \(A\), \(B\), \(B^ {\text{T}}\), and \(-C\), where \(A\) and \(-C\) lie on the main diagonal, and where, moreover, \(A\) is a symmetric positive (SP) definite square matrix and \(C\) is a SP semidefinite square matrix. Two variants of the method are derived: a direct approach, and a preconditioner-based technique. A numerical experiment (Stokes equations) is used to compare the performance of both variants of the proposed method with other three error estimation algorithms. The results indicate that the proposed method delivers less sharp error bounds than the three reference algorithms do, but it is significantly faster and can be applied to large systems, where the reference algorithms failed or ran out of computer memory.











This page was built for publication: Numerical validation of solutions of saddle point matrix equations

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5317905)