A Parallel Graph Algorithm for Detecting Mesh Singularities in Distributed Memory Ice Sheet Simulations

We present a new, distributed-memory parallel algorithm for detection of degenerate mesh features that can cause singularities in ice sheet mesh simulations. Identifying and removing mesh features such as disconnected components (icebergs) or hinge vertices (peninsulas of ice detached from the land) can significantly improve the convergence of iterative solvers. Because the ice sheet evolves during the course of a simulation, it is important that the detection algorithm can run in situ with the simulation --- running in parallel and taking a negligible amount of computation time --- so that degenerate features (e.g., calving icebergs) can be detected as they develop. We present a distributed memory, BFS-based label-propagation approach to degenerate feature detection that is efficient enough to be called at each step of an ice sheet simulation, while correctly identifying all degenerate features of an ice sheet mesh. Our method finds all degenerate features in a mesh with 13 million vertices in 0.0561 seconds on 1536 cores in the MPAS Albany Land Ice (MALI) model. Compared to the previously used serial pre-processing approach, we observe a 46,000x speedup for our algorithm, and provide additional capability to do dynamic detection of degenerate features in the simulation.

[1]  Kishore Kothapalli,et al.  Efficient Multicore Algorithms For Identifying Biconnected Components , 2016, Int. J. Netw. Comput..

[2]  Kesheng Wu,et al.  Optimizing two-pass connected-component labeling algorithms , 2009, Pattern Analysis and Applications.

[3]  R. Tuminaro,et al.  MPAS-Albany Land Ice (MALI): a variable-resolution ice sheet model for Earth system modeling using Voronoi grids , 2018, Geoscientific Model Development.

[4]  Robert E. Tarjan,et al.  An Efficient Parallel Biconnectivity Algorithm , 2011, SIAM J. Comput..

[5]  Daniel F. Martin,et al.  Adaptive mesh, finite volume modeling of marine ice sheets , 2013, J. Comput. Phys..

[6]  B. Nour-Omid,et al.  A study of the factorization fill‐in for a parallel implementation of the finite element method , 1994 .

[7]  Andrew J. Plater,et al.  Book reviewSea-level change: Roger Revelle; Studies in Geophysics, National Research Council, National Academy Press, Washington, DC, 1990; xii + 246 pp.; USD 29.95, GBP 25.75; ISBN 0-309-04039 , 1992 .

[8]  J. Hopcroft,et al.  Efficient algorithms for graph manipulation , 1971 .

[9]  Kamesh Madduri,et al.  Simple parallel biconnectivity algorithms for multicore platforms , 2014, 2014 21st International Conference on High Performance Computing (HiPC).

[10]  Houjun Tang,et al.  Parallel In Situ Detection of Connected Components in Adaptive Mesh Refinement Data , 2015, 2015 15th IEEE/ACM International Symposium on Cluster, Cloud and Grid Computing.

[11]  Andrew G. Salinger,et al.  Albany/FELIX : a parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis , 2014 .

[12]  Michael A. Heroux,et al.  A new overview of the Trilinos project , 2012, Sci. Program..

[13]  Celso C. Ribeiro,et al.  Power optimization in ad hoc wireless network topology control with biconnectivity requirements , 2013, Comput. Oper. Res..

[14]  J. Hopcroft,et al.  Algorithm 447: efficient algorithms for graph manipulation , 1973, CACM.

[15]  Eric Rignot,et al.  Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM) , 2012 .

[16]  Garth M. Reese,et al.  Mechanism free domain decomposition , 2003 .

[17]  Kelly P. Gaither,et al.  A Distributed-Memory Algorithm for Connected Components Labeling of Simulation Data , 2015, Topological and Statistical Methods for Complex Data, Tackling Large-Scale, High-Dimensional, and Multivariate Data Spaces.

[18]  Mehmet Deveci,et al.  Zoltan2: Next-Generation Combinatorial Toolkit. , 2012 .

[19]  Michael A. Heroux,et al.  Tpetra, and the use of generic programming in scientific computing , 2012, Sci. Program..

[20]  Andrew G. Salinger,et al.  A Matrix Dependent/Algebraic Multigrid Approach for Extruded Meshes with Applications to Ice Sheet Modeling , 2016, SIAM J. Sci. Comput..