A new parallel sparse direct solver: Presentation and numerical experiments in large‐scale structural mechanics parallel computing

The main purpose of this work is to present a new parallel direct solver: Dissection solver. It is based on LU factorization of the sparse matrix of the linear system and allows to detect automatically and handle properly the zero-energy modes, which are important when dealing with DDM. A performance evaluation and comparisons with other direct solvers (MUMPS, DSCPACK) are also given for both sequential and parallel computations. Results of numerical experiments with a two-level parallelization of large-scale structural analysis problems are also presented: FETI is used for the global problem parallelization and Dissection for the local multithreading. In this framework, the largest problem we have solved is of an elastic solid composed of 400 subdomains running on 400 computation nodes (3200 cores) and containing about 165 millions dof. The computation of one single iteration consumes less than 20 min of CPU time. Several comparisons to MUMPS are given for the numerical computation of large-scale linear systems on a massively parallel cluster: performances and weaknesses of this new solver are highlighted. Copyright © 2011 John Wiley & Sons, Ltd.

[1]  C. Farhat,et al.  A method of finite element tearing and interconnecting and its parallel solution algorithm , 1991 .

[2]  Patrick R. Amestoy,et al.  Multifrontal parallel distributed symmetric and unsymmetric solvers , 2000 .

[3]  Ibrahima Gueye,et al.  Résolution de grands systèmes linéaires issus de la méthode des éléments finis sur des calculateurs massivement parallèles. (Solving large linear systems arising in finite element approximations on massively parallel computers) , 2009 .

[4]  J. Pasciak,et al.  Computer solution of large sparse positive definite systems , 1982 .

[5]  Frédéric Feyel,et al.  A parallel algorithm for direct solution of large sparse linear systems, well suitable to domain decomposition methods , 2009 .

[6]  Charbel Farhat,et al.  Implicit parallel processing in structural mechanics , 1994 .

[7]  Tarek P. Mathew,et al.  Domain Decomposition Methods for the Numerical Solution of Partial Differential Equations , 2008, Lecture Notes in Computational Science and Engineering.

[8]  J. Mandel.  Balancing domain decomposition , 1993 .

[9]  Jack J. Dongarra,et al.  A set of level 3 basic linear algebra subprograms , 1990, TOMS.

[10]  A. George.  Nested Dissection of a Regular Finite Element Mesh , 1973 .

[11]  T. Chan,et al.  Domain decomposition algorithms , 1994, Acta Numerica.

[12]  P. Tallec.  Domain decomposition methods in computational mechanics , 1994 .

[13]  C. Farhat,et al.  The second generation FETI methods and their application to the parallel solution of large-scale linear and geometrically non-linear structural analysis problems , 2000 .

[14]  P. Gosselet,et al.  Non-overlapping domain decomposition methods in structural mechanics , 2006, 1208.4209.

[15]  J. W. Walker,et al.  Direct solutions of sparse network equations by optimally ordered triangular factorization , 1967 .

[16]  Andrea Toselli,et al.  Domain decomposition methods : algorithms and theory , 2005 .