skip to content

Date: December 7, 2022, 4–5:30 p.m. CET. Coffee and tea from 4–4:30 p.m.

Speaker: Prof. Dr. Martin Kronbichler (University of Augsburg), ORCID


     Weyertal 121, 50931 Cologne
     RRZK building (Google Maps, OpenStreetMap), Entrance
     Room 1.03

Title: Efficient Implementations of High-Order Finite-Element Methods

Abstract: The talk focuses on the efficient implementation of high-order methods in the context of continuous and discontinuous finite elements applied to linear and nonlinear partial differential equations on complex geometries. The core algorithmic ingredient is the matrix-free operator evaluation (matrix-vector product), for which fast quadrature schemes for cell and face integrals based on a technique called sum factorization are employed. The algorithm selection for achieving a high throughput is demonstrated, guided by performance analysis, as well as the scalability to large-scale parallel machines. The best algorithms come with an arithmetic intensity of one to five flop/byte, with memory transfer primarily due to the access into input and output vectors as well as some geometry or variable coefficient data. As a result of our optimizations, many downstream solvers, such as explicit time stepping, smoothers in multigrid methods, or conjugate gradient solvers, are now no longer dominated by the matrix-vector product but by vector operations instead. To overcome this limitation, we consider the whole solution chain and fuse the memory access between different stages of an algorithm. The effect of these optimizations on the performance on contemporary CPU and GPU architectures are illustrated and results are shown of fluid dynamics simulations implemented with computational kernels from the deal.II finite element library.