In this paper we present the performance of our parallel multi-frontal direct solver when applied to solve linear systems of equations resulting from discretizations of a Finite Element Method ( -FEM). The -FEM generates a sequence of computational meshes delivering exponential convergence of the numerical error with respect to the mesh size (number of degrees of freedom). A sequence of meshes is obtained by performing several refinements starting from an arbitrary initial mesh. The solver constructs initial elimination tree for an arbitrary initial mesh, and expands the elimination tree each time the mesh is refined. The solver has been tested on 3D Direct Current (DC) borehole resistivity measurement simulations problems. We compare the solver with two versions of the MUMPS parallel solver: with (1) distributed entries executed over the entire problem, and (2) the direct sub-structuring method with parallel MUMPS solver utilized to solve the interface problem. We show that by providing to the solver the knowledge about the structure of the -FEM, the order of elimination is obtained straightforward, and leads to a better performance than by submitting the entire matrix to the solver and executing a connectivity graph based ordering algorithm.
[1]
Luc Giraud,et al.
Iterative versus direct parallel substructuring methods in semiconductor device modelling
,
2005,
Numer. Linear Algebra Appl..
[2]
Victor M. Calo,et al.
Fourier series expansion in a non-orthogonal system of coordinates for the simulation of 3D-DC borehole resistivity measurements
,
2008
.
[3]
J. T. Oden,et al.
A parallel multifrontal algorithm and its implementation
,
1997
.
[4]
Bruce M. Irons,et al.
A frontal solution program for finite element analysis
,
1970
.
[5]
Leszek Demkowicz,et al.
Parallel, fully automatic hp-adaptive 2d finite element package
,
2006
.
[6]
Vipin Kumar,et al.
A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs
,
1998,
SIAM J. Sci. Comput..
[7]
Maciej Paszynski.
Performance of Multi Level Parallel Direct Solver for hpFinite Element Method
,
2007,
PPAM.