A hybrid parallel solver for finite element computations

Abstract

We present a parallel hybrid solver called dual Schur for finite element computations based on domain decomposition. Each subdomain is allocated to one-processor and the linear systems internal to the subdomains are solved with a direct method (LU or Cholesky). The interface problem is expressed in terms of dual variables and uses an iterative algorithm.

Typical results are illustrated for the solution of a Navier-Stokes flow on the Convex Exemplar SPP system. We also emphasize the influence of the decomposition quality on the parallel performances. A good decomposition must exhibit a high load balance with a low interface size and it should contain “compact” subdomains.

Presenting author. The Ph.D. work of Denis Vanderstraeten is supported by the Fonds National de la Recherche Scientifique of the Belgian State. This paper presents research results of the Belgian Programme on Interuniversity Poles of Attraction, initiated by the Belgian State, Prime Minister's Office for Science, Technology and Culture. The scientific responsibility rests with its authors.