An Optimal Iterative Solver for Symmetric Indefinite Systems stemming from Mixed Approximation

Silvester, David J. and Simoncini, Valeria (2011) An Optimal Iterative Solver for Symmetric Indefinite Systems stemming from Mixed Approximation. ACM Transactions on Mathematical Software, 37 (4). ISSN 1749-9097

This is the latest version of this item.

[thumbnail of a42-silvester.pdf] PDF

Download (458kB)


We discuss the design and implementation of a suite of functions for solving symmetric indefinite linear systems associated with mixed approximation of systems of PDEs. The novel feature of our iterative solver is the incorporation of error control in the natural "energy" norm in combination with an a posteriori estimator for the PDE approximation error. This leads to a robust and optimally efficient stopping criterion: the iteration is terminated as soon as the algebraic error is insignificant compared to the approximation error. We describe a "proof of concept" MATLAB implementation of this algorithm and we illustrate its effectiveness when integrated into the Incompressible Flow Iterative Solution Software (IFISS) package (cf. ACM Transactions on Mathematical Software 33, Article 14, 2007).

Item Type: Article
Additional Information: This is the published version of this manuscript.
Uncontrolled Keywords: Finite elements, incompressible flow, iterative solvers, stopping criteria, EST_MINRES, MATLAB
Subjects: MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis
Depositing User: professor david silvester
Date Deposited: 23 Apr 2011
Last Modified: 20 Oct 2017 14:12

Available Versions of this Item

Actions (login required)

View Item View Item