GPU accelerated spectral finite elements on all-hex meshes
Abstract
This paper presents a spectral element finite element scheme that efficiently solves elliptic problems on unstructured hexahedral meshes. The discrete equations are solved using a matrix-free preconditioned conjugate gradient algorithm. An additive Schwartz two-scale preconditioner is employed that allows h-independence convergence. An extensible multi-threading programming API is used as a common kernel language that allows runtime selection of different computing devices (GPU and CPU) and different threading interfaces (CUDA, OpenCL and OpenMP). Performance tests demonstrate that problems with over 50 million degrees of freedom can be solved in a few seconds on an off-the-shelf GPU.
- Publication:
-
Journal of Computational Physics
- Pub Date:
- November 2016
- DOI:
- 10.1016/j.jcp.2016.08.005
- arXiv:
- arXiv:1506.05996
- Bibcode:
- 2016JCoPh.324..246R
- Keywords:
-
- Spectral finite elements;
- GPU computing;
- Hexahedral meshes;
- Computer Science - Computational Engineering;
- Finance;
- and Science;
- Computer Science - Distributed;
- Parallel;
- and Cluster Computing;
- Computer Science - Numerical Analysis;
- 65Y05;
- 65Y10;
- 65Y20
- E-Print:
- 23 pages, 7 figures