Assembly of finite element methods on graphics processors

Cris Cecka, Adrian J. Lew, E. Darve

Research output: Contribution to journalArticlepeer-review

173 Scopus citations


Recently, graphics processing units (GPUs) have had great success in accelerating many numerical computations. We present their application to computations on unstructured meshes such as those in finite element methods. Multiple approaches in assembling and solving sparse linear systems with NVIDIA GPUs and the Compute Unified Device Architecture (CUDA) are created and analyzed. Multiple strategies for efficient use of global, shared, and local memory, methods to achieve memory coalescing, and optimal choice of parameters are introduced. We find that with appropriate preprocessing and arrangement of support data, the GPU coprocessor using single-precision arithmetic achieves speedups of 30 or more in comparison to a well optimized double-precision single core implementation. We also find that the optimal assembly strategy depends on the order of polynomials used in the finite element discretization. © 2010 John Wiley & Sons, Ltd.
Original languageEnglish (US)
Pages (from-to)640-669
Number of pages30
JournalInternational Journal for Numerical Methods in Engineering
Issue number5
StatePublished - Aug 23 2010
Externally publishedYes


Dive into the research topics of 'Assembly of finite element methods on graphics processors'. Together they form a unique fingerprint.

Cite this