Saddle Point Least Squares Preconditioning of Mixed Methods
This work provides a theoretically grounded method for efficiently solving mixed finite element problems, which is relevant to computational scientists working on PDEs.
The paper presents a new discretization and preconditioning approach for mixed variational formulations that avoids assembling a global saddle point system, with proven sharp approximation properties and convergence rate estimates.
We present a simple way to discretize and precondition mixed variational formulations. Our theory connects with, and takes advantage of, the classical theory of symmetric saddle point problems and the theory of preconditioning symmetric positive definite operators. Efficient iterative processes for solving the discrete mixed formulations are proposed and choices for discrete spaces that are always compatible are provided. For the proposed discrete spaces and solvers, a basis is needed only for the test spaces and assembly of a global saddle point system is avoided. We prove sharp approximation properties for the discretization and iteration errors and also provide a sharp estimate for the convergence rate of the proposed algorithm in terms of the condition number of the elliptic preconditioner and the discrete $\inf-\sup$ and $\sup-\sup$ constants of the pair of discrete spaces.