Parallel application of a novel domain decomposition preconditioner for the stable finite element solution of three-dimensional convection-dominated PDEs

We describe and analyze the parallel implementation of a novel domain decomposition preconditioner for the fast iterative solution of linear systems of algebraic equations arising from the discretization of elliptic partial differential equations (PDEs) in three dimensions. In previous theoretical work, [4], this preconditioner has been proved to be optimal for symmetric positive-definite (SPD) linear systems. In this paper we provide details of our 3-d parallel implementation and demonstrate that the technique may be generalized to the solution of non-symmetric algebraic systems, such as those arising when convection-diffusion problems are discretized using either Galerkin or stabilized finite element methods (FEMs), [11].