We report on a two-scale approach for efficient matrix-free finite element simulations. The proposed method is based on surrogate element matrices constructed by low-order polynomial approximations. It is applied to a Stokes-type PDE system with variable viscosity as is a key component in mantle convection models. We set the ground for a rigorous performance analysis inspired by the concept of parallel textbook multigrid efficiency and study the weak scaling behavior on SuperMUC, a peta-scale supercomputer system. For a complex geodynamical model, we achieve a parallel efficiency of 95% on up to 47 250 compute cores. Our largest simulation uses a trillion (O1012)) degrees of freedom for a global mesh resolution of 1.7 km.
CITATION STYLE
Bauer, S., Huber, M., Mohr, M., Rüde, U., & Wohlmuth, B. (2018). A New Matrix-Free Approach for Large-Scale Geodynamic Simulations and its Performance. In Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics) (Vol. 10861 LNCS, pp. 17–30). Springer Verlag. https://doi.org/10.1007/978-3-319-93701-4_2
Mendeley helps you to discover research relevant for your work.