An adaptive finite element multigrid solver using GPU acceleration

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Liebchen, Manuel, Jendersie, Robert, Kaya, Utku, Lessig, Christian, Richter, Thomas
Format: Preprint
Publié: 2024
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866914211284647936
author Liebchen, Manuel
Jendersie, Robert
Kaya, Utku
Lessig, Christian
Richter, Thomas
author_facet Liebchen, Manuel
Jendersie, Robert
Kaya, Utku
Lessig, Christian
Richter, Thomas
contents Adaptive finite elements combined with geometric multigrid solvers are one of the most efficient numerical methods for problems such as the instationary Navier-Stokes equations. Yet despite their efficiency, computations remain expensive and the simulation of, for example, complex flow problems can take many hours or days. GPUs provide an interesting avenue to speed up the calculations due to their very large theoretical performance. However, the large degree of parallelism and non-standard API make the use of GPUs in scientific computing challenging. In this work, we develop a GPU acceleration for the adaptive finite element library Gascoigne and study its effectiveness for different systems of partial differential equations. Our goal is thereby to integrate the GPU acceleration into the existing code with minimal changes, even when this requires a penalty in the GPU acceleration. Through the systematic formulation of all computations as linear algebra operations, we can employ GPU-accelerated linear algebra libraries, which simplifies the implementation and ensures the maintainability of the code while achieving very efficient GPU utilizations. Our results for a transport-diffusion equation, linear elasticity, and the instationary Navier-Stokes equations show substantial speedups of up to 20X compared to multi-core CPU implementations.
format Preprint
id arxiv_https___arxiv_org_abs_2405_05047
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle An adaptive finite element multigrid solver using GPU acceleration
Liebchen, Manuel
Jendersie, Robert
Kaya, Utku
Lessig, Christian
Richter, Thomas
Numerical Analysis
Adaptive finite elements combined with geometric multigrid solvers are one of the most efficient numerical methods for problems such as the instationary Navier-Stokes equations. Yet despite their efficiency, computations remain expensive and the simulation of, for example, complex flow problems can take many hours or days. GPUs provide an interesting avenue to speed up the calculations due to their very large theoretical performance. However, the large degree of parallelism and non-standard API make the use of GPUs in scientific computing challenging. In this work, we develop a GPU acceleration for the adaptive finite element library Gascoigne and study its effectiveness for different systems of partial differential equations. Our goal is thereby to integrate the GPU acceleration into the existing code with minimal changes, even when this requires a penalty in the GPU acceleration. Through the systematic formulation of all computations as linear algebra operations, we can employ GPU-accelerated linear algebra libraries, which simplifies the implementation and ensures the maintainability of the code while achieving very efficient GPU utilizations. Our results for a transport-diffusion equation, linear elasticity, and the instationary Navier-Stokes equations show substantial speedups of up to 20X compared to multi-core CPU implementations.
title An adaptive finite element multigrid solver using GPU acceleration
topic Numerical Analysis
url https://arxiv.org/abs/2405.05047