A Static Parallel Multifrontal Solver for Finite Element Meshes [chapter]

Alberto Bertoldo, Mauro Bianco, Geppino Pucci
2006 Lecture Notes in Computer Science  
We present a static parallel implementation of the multifrontal method to solve unsymmetric sparse linear systems on distributed-memory architectures. We target Finite Element (FE) applications where numerical pivoting can be avoided, since an implicit minimum-degree ordering based on the FE mesh topology suffices to achieve numerical stability. Our strategy is static in the sense that work distribution and communication patterns are determined in a preprocessing phase preceding the actual
more » ... ical computation. To balance the load among the processors, we devise a simple model-driven partitioning strategy to precompute a high-quality balancing for a large family of structured meshes. The resulting approach is proved to be considerably more efficient than the strategies implemented by MUMPS and SuperLU DIST, two state-of-the-art parallel multifrontal solvers.
doi:10.1007/11946441_67 fatcat:w2hbrt2blfhhzasj3pimh26un4