Hypergrammar Based Parallel Multi-frontal Solver for Grids with Point singularities

This paper describes the application of hypergraph grammars to drive a linear computational cost solver for grids with point singularities. Such graph gram- mar productions are the first mathematical formalisms used to describe solver algorithms, and each indicates the smallest atomic task that can be executed in parallel, which is very useful in the case of parallel execution. In particular, the partial order of execution of graph grammar productions can be found, and the sets of independent graph grammar productions can be localized. They can be scheduled set by set into a shared memory parallel machine. The graph- grammar-based solver has been implemented with NVIDIA CUDA for GPU. Graph grammar productions are accompanied by numerical results for a 2D case. We show that our graph-grammar-based solver with a GPU accelerator is, by order of magnitude, faster than the state-of-the-art MUMPS solver.

[1]  Annegret Habel,et al.  May we introduce to you: hyperedge replacement , 1986, Graph-Grammars and Their Application to Computer Science.

[2]  Grazyna Slusarczyk,et al.  Hypergraph Grammars in hp-Adaptive Finite Element Method , 2013, ICCS.

[3]  Mikaël Barboteu,et al.  An analytical and numerical approach to a bilateral contact problem with nonmonotone friction , 2013, Int. J. Appl. Math. Comput. Sci..

[4]  Ewa Grabska,et al.  A Graph Grammar Model of the hp Adaptive Three Dimensional Finite Element Method. Part II , 2012, Fundam. Informaticae.

[5]  Bruce M. Irons,et al.  A frontal solution program for finite element analysis , 1970 .

[6]  John K. Reid,et al.  The Multifrontal Solution of Indefinite Sparse Symmetric Linear , 1983, TOMS.

[7]  Robert Schaefer,et al.  Graph grammar‐driven parallel partial differential equation solver , 2010, Concurr. Comput. Pract. Exp..

[8]  L. Demkowicz.  One and two dimensional elliptic and Maxwell problems , 2006 .

[9]  Vipin Kumar,et al.  Highly Scalable Parallel Algorithms for Sparse Matrix Factorization , 1997, IEEE Trans. Parallel Distributed Syst..

[10]  Sivan Toledo,et al.  The design, implementation, and evaluation of a symmetric banded linear solver for distributed-memory parallel computers , 1998, TOMS.

[11]  Anna Paszynska,et al.  Graph Grammar Based Direct Solver for Hp-adaptive Finite Element Method with Point Singularities , 2013, ICCS.

[12]  Ewa Grabska,et al.  A Graph Grammar Model of the hp Adaptive Three Dimensional Finite Element Method. Part I , 2012, Fundam. Informaticae.

[13]  Ewa Grabska,et al.  Graph Transformations for Modeling hp-Adaptive Finite Element Method with Mixed Triangular and Rectangular Elements , 2009, ICCS.

[14]  Lexing Ying,et al.  A fast direct solver for elliptic problems on general meshes in 2D , 2012, J. Comput. Phys..

[15]  Victor M. Calo,et al.  The cost of continuity: A study of the performance of isogeometric finite elements using direct solvers , 2012 .

[16]  Maciej Paszynski.  On the Parallelization of Self-Adaptive hp-Finite Element Methods Part I. Composite Programmable Graph GrammarModel , 2009, Fundam. Informaticae.

[17]  Iain S. Duff,et al.  The Multifrontal Solution of Unsymmetric Sets of Linear Equations , 1984 .

[18]  Grzegorz Rozenberg,et al.  Handbook of Graph Grammars and Computing by Graph Transformations, Volume 1: Foundations , 1997 .

[19]  Patrick Hild,et al.  A sign preserving mixed finite element approximation for contact problems , 2011, Int. J. Appl. Math. Comput. Sci..

[20]  Annegret Habel,et al.  Some Structural Aspects of Hypergraph Languages Generated by Hyperedge Replacement , 1987, STACS.

[21]  Timothy A. Davis,et al.  Dynamic Supernodes in Sparse Cholesky Update/Downdate and Triangular Solves , 2009, TOMS.

[22]  Marcin Sieniek,et al.  Agent-oriented image processing with the hp-adaptive projection-based interpolation operator , 2011, ICCS.

[23]  Maciej Paszyński,et al.  Computing with hp-ADAPTIVE FINITE ELEMENTS: Volume II Frontiers: Three Dimensional Elliptic and Maxwell Problems with Applications , 2007 .

[24]  Victor M. Calo,et al.  A direct solver with reutilization of LU factorizations for hh-adaptive finite element grids with point singularities , 2013, Comput. Math. Appl..

[25]  T. Hughes,et al.  ISOGEOMETRIC ANALYSIS: APPROXIMATION, STABILITY AND ERROR ESTIMATES FOR h-REFINED MESHES , 2006 .

[26]  Thomas J. R. Hughes,et al.  Isogeometric Analysis: Toward Integration of CAD and FEA , 2009 .

[27]  Hartmut Ehrig,et al.  Handbook of graph grammars and computing by graph transformation: vol. 3: concurrency, parallelism, and distribution , 1999 .

[28]  Anshul Gupta,et al.  Recent advances in direct methods for solving unsymmetric sparse systems of linear equations , 2002, TOMS.

[29]  Otto von Estorff,et al.  BEM and FEM results of displacements in a poroelastic column , 2012, Int. J. Appl. Math. Comput. Sci..