Past Event: Oden Institute Seminar

MFEM: Accelerating Efficient Solution of PDEs at Exascale

Tzanio Kolev, Computational Mathematician, Center for Applied Scientific Computing, Lawrence Livermore National Laboratory

Abstract

Upcoming exascale architectures require rethinking of the numerical algorithms used in large-scale PDE-based applications. These architectures favor algorithms, such as high-order finite elements, that expose fine-grain parallelism and maximize the ratio of floating point operations to energy intensive data movement. In this talk we present an overview of MFEM [1], a scalable library for high-order finite element discretization of PDEs on general unstructured grids. We also report on recent work in the Center for Efficient Exascale Discretizations [2], a co-design center in the US Exascale Computing Project focused on next-generation discretization software and algorithms. Our approach to efficient operator evaluation is based on a "matrix-free" representation of the finite element operator, that factors a bilinear form into a series of sparse and dense components corresponding to the parallelism, mesh topology, basis, geometry, and pointwise physics in the problem. The operator decomposition exposes several layers of parallelism, enables the use of batched dgemss and tensor contractions, and only requires quadrature point values to be assembled for computing the action. This "partial assembly" formulation is a natural fit for modern HPC hardware, because it results both in less (nearly optimal) computation and less (optimal) data movement compared to assembling a global sparse matrix, therefore increasing performance and reducing time to solution. In addition to discussing efficient operator evaluation, we will provide an overview of the MFEM capabilities and applications to compressible hydrodynamics and electromagnetics. We will also review our work on performance optimizations for GPU architectures, high-order benchmarks and miniapps, scalable unstructured adaptive mesh refinement, high-order mesh optimization and matrix-free preconditioning. [1] MFEM: Modular finite element library, http://mfem.org. [2] Center for Efficient Exascale Discretizations, http://ceed.exascaleproject.org. Bio Tzanio Kolev is a computational mathematician at the Center for Applied Scientific Computing (CASC) in Lawrence Livermore National Laboratory, where he works on scalable high-order finite element discretizations and solvers for problems in compressible shock hydrodynamics, multi-material ALE, computational electromagnetics, and more. He is the project leader of the MFEM project at LLNL and the director of the Center for Efficient Exascale Discretizations (CEED) in DOE’s Exascale Computing Project. Tzanio's research interests include the development and analysis of finite element methods, massively parallel preconditioners, discretization-enhanced algebraic multigrid algorithms, and the design and implementation of large-scale scientific software. **Note: Please join this Zoom seminar online with the "Audio Only" function (no video)**