Performant low-order matrix-free finite element kernels on GPU architectures

Randolph R. Settgast, Yohann Dudouit, Nicola Castelletto, William R. Tobin, Benjamin C. Corbett, Sergey Klevtsov
Atmospheric, Earth, and Energy Division, Lawrence Livermore National Laboratory, Livermore, CA, United States
arXiv:2308.09839 [math.NA], (22 Aug 2023)


   title={Performant low-order matrix-free finite element kernels on GPU architectures},

   author={Randolph R. Settgast and Yohann Dudouit and Nicola Castelletto and William R. Tobin and Benjamin C. Corbett and Sergey Klevtsov},






Download Download (PDF)   View View   Source Source   Source codes Source codes




Numerical methods such as the Finite Element Method (FEM) have been successfully adapted to utilize the computational power of GPU accelerators. However, much of the effort around applying FEM to GPU’s has been focused on high-order FEM due to higher arithmetic intensity and order of accuracy. For applications such as the simulation of subsurface processes, high levels of heterogeneity results in high-resolution grids characterized by highly discontinuous (cell-wise) material property fields. Moreover, due to the significant uncertainties in the characterization of the domain of interest, e.g. geologic reservoirs, the benefits of high order accuracy are reduced, and low-order methods are typically employed. In this study, we present a strategy for implementing highly performant low-order matrix-free FEM operator kernels in the context of the conjugate gradient (CG) method. Performance results of matrix-free Laplace and isotropic elasticity operator kernels are presented and are shown to compare favorably to matrix-based SpMV operators on V100, A100, and MI250X GPUs.
No votes yet.
Please wait...

Recent source codes

* * *

* * *

HGPU group © 2010-2024 hgpu.org

All rights belong to the respective authors

Contact us: