Exploiting semantics of temporal multi‐scale methods to optimize multi‐level mesh partitioning

Multi-scale problems are often solved by decomposing the problem domain into multiple subdomains, solving them independently using different levels of spatial and temporal refinement, and coupling the subdomain solutions back to obtain the global solution. Most commonly, finite elements are used for spatial discretization, and finite difference time stepping is used for time integration. Given a finite element mesh for the global problem domain, the number of possible decompositions into subdomains and the possible choices for associated time steps is exponentially large, and the computational costs associated with different decompositions can vary by orders of magnitude. The problem of finding an optimal decomposition and the associated time discretization that minimizes computational costs while maintaining accuracy is nontrivial. Existing mesh partitioning tools, such as METIS, overlook the constraints posed by multi-scale methods and lead to suboptimal partitions with a high performance penalty. We present a multi-level mesh partitioning approach that exploits domain-specific knowledge of multi-scale methods to produce nearly optimal mesh partitions and associated time steps automatically. Results show that for multi-scale problems, our approach produces decompositions that outperform those produced by state-of-the-art partitioners like METIS and even those that are manually constructed by domain experts. Copyright © 2017 John Wiley & Sons, Ltd.

[1]  Bradley C. Kuszmaul,et al.  Cilk: an efficient multithreaded runtime system , 1995, PPOPP '95.

[2]  Vipin Kumar,et al.  A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs , 1998, SIAM J. Sci. Comput..

[3]  Charbel Farhat,et al.  Modeling and Simulation of Multiphysics Systems , 2005, J. Comput. Inf. Sci. Eng..

[4]  Arun Prakash,et al.  Computationally efficient multi-time-step method for partitioned time integration of highly nonlinear structural dynamics , 2014 .

[5]  Arun Prakash,et al.  Exploiting domain knowledge to optimize parallel computational mechanics codes , 2013, ICS '13.

[6]  Jun Zhang,et al.  BILUTM: A Domain-Based Multilevel Block ILUT Preconditioner for General Sparse Matrices , 1999, SIAM J. Matrix Anal. Appl..

[7]  Bruce Hendrickson,et al.  A Multi-Level Algorithm For Partitioning Graphs , 1995, Proceedings of the IEEE/ACM SC95 Conference.

[8]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[9]  Wolfgang A. Wall,et al.  Towards a taxonomy for multiscale methods in computational mechanics: building blocks of existing methods , 2007 .

[10]  Richard B. Lehoucq,et al.  An Automated Multilevel Substructuring Method for Eigenspace Computation in Linear Elastodynamics , 2004, SIAM J. Sci. Comput..

[11]  Chao-Fu Wang,et al.  A bubble‐inspired algorithm for finite element mesh partitioning , 2013 .

[12]  Jinchao Xu,et al.  Some Nonoverlapping Domain Decomposition Methods , 1998, SIAM Rev..

[13]  Anthony Gravouil,et al.  Heterogeneous asynchronous time integrators for computational structural dynamics , 2015 .

[14]  Jacob Fish,et al.  Bridging the scales in nano engineering and science , 2006 .

[15]  A. Prakash,et al.  A FETI‐based multi‐time‐step coupling method for Newmark schemes in structural dynamics , 2004 .

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

[17]  Masha Sosonkina,et al.  pARMS: a parallel version of the algebraic recursive multilevel solver , 2003, Numer. Linear Algebra Appl..

[18]  Wolfgang Hackbusch,et al.  On the Computation of Approximate Eigenvalues and Eigenfunctions of Elliptic Operators by Means of a Multi-Grid Method , 1979 .

[19]  U. Hetmaniuk,et al.  A comparison of eigensolvers for large‐scale 3D modal analysis using AMG‐preconditioned iterative methods , 2005 .

[20]  Manolis Papadrakakis,et al.  The mosaic of high performance domain Decomposition Methods for Structural Mechanics: Formulation, interrelation and numerical efficiency of primal and dual methods , 2003 .