A New Matrix-Free Approach for Large-Scale Geodynamic Simulations and its Performance

被引:9
|
作者
Bauer, Simon [1 ]
Huber, Markus [2 ]
Mohr, Marcus [1 ]
Ruede, Ulrich [3 ,4 ]
Wohlmuth, Barbara [2 ]
机构
[1] Ludwig Maximilians Univ Munchen, Dept Earth & Environm Sci, Munich, Germany
[2] Tech Univ Munich, Inst Numer Math M2, Munich, Germany
[3] FAU Erlangen Nurnberg, Dept Comp Sci 10, Erlangen, Germany
[4] CERFACS, Parallel Algorithms Project, Toulouse, France
来源
关键词
Two-scale PDE discretization; Massively parallel multigrid; Matrix-free on-the-fly assembly; Large scale geophysical application; HIERARCHICAL HYBRID GRIDS; MANTLE FLOW; PARALLEL; MODELS; SOLVERS;
D O I
10.1007/978-3-319-93701-4_2
中图分类号
TP301 [理论、方法];
学科分类号
081202 ;
摘要
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 (O(10(12))) degrees of freedom for a global mesh resolution of 1.7 km.
引用
收藏
页码:17 / 30
页数:14
相关论文
共 50 条
  • [1] Large-scale simulation of mantle convection based on a new matrix-free approach
    Bauer, S.
    Huber, M.
    Ghelichkhan, S.
    Mohr, M.
    Ruede, U.
    Wohlmuth, B.
    [J]. JOURNAL OF COMPUTATIONAL SCIENCE, 2019, 31 : 60 - 76
  • [2] Matrix-free large-scale Bayesian inference in cosmology
    Jasche, Jens
    Lavaux, Guilhem
    [J]. MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, 2015, 447 (02) : 1204 - 1212
  • [3] A new matrix-free algorithm for the large-scale trust-region subproblem
    Rojas, M
    Santos, SA
    Sorensen, DC
    [J]. SIAM JOURNAL ON OPTIMIZATION, 2001, 11 (03) : 611 - 646
  • [4] A matrix-free cone complementarity approach for solving large-scale, nonsmooth, rigid body dynamics
    Tasora, A.
    Anitescu, M.
    [J]. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2011, 200 (5-8) : 439 - 453
  • [5] A matrix-free smoothing algorithm for large-scale support vector machines
    Ni, Tie
    Zhai, Jun
    [J]. INFORMATION SCIENCES, 2016, 358 : 29 - 43
  • [6] A matrix-free approach to build band preconditioners for large-scale bound-constrained optimization
    De Simone, V.
    di Serafino, D.
    [J]. JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2014, 268 : 82 - 92
  • [7] ALPAQA: A matrix-free solver for nonlinear MPC and large-scale nonconvex optimization
    Pas, Pieter
    Schuurmans, Mathijs
    Patrinos, Panagiotis
    [J]. 2022 EUROPEAN CONTROL CONFERENCE (ECC), 2022, : 417 - 422
  • [8] smashGP: Large-Scale Spatial Modeling via Matrix-Free Gaussian Processes
    Erlandson, Lucas
    Gomez, Ana Maria Estrada
    Chow, Edmond
    Paynabar, Kamran
    [J]. JOURNAL OF COMPUTATIONAL AND GRAPHICAL STATISTICS, 2024,
  • [9] Matrix-free algorithm for the large-scale constrained trust-region subproblem
    Kearsley, AJ
    [J]. OPTIMIZATION METHODS & SOFTWARE, 2006, 21 (02): : 233 - 245
  • [10] A matrix-free augmented lagrangian algorithm with application to large-scale structural design optimization
    Arreckx, Sylvain
    Lambe, Andrew
    Martins, Joaquim R. R. A.
    Orban, Dominique
    [J]. OPTIMIZATION AND ENGINEERING, 2016, 17 (02) : 359 - 384