Finite element numerical integration for first order approximations on multi- and many-core architectures

The paper presents investigations on the performance of the finite element numerical integration algorithm for first order approximations and three processor architectures, popular in scientific computing, classical x86_64 CPU, Intel Xeon Phi and NVIDIA Kepler GPU. We base the discussion on theoreti...

Full description

Saved in:
Bibliographic Details
Published in:Computer methods in applied mechanics and engineering Vol. 305; pp. 827 - 848
Main Authors: Banasa, Krzysztof, Kruzelb, Filip, Bielanskia, Jan
Format: Journal Article
Language:English
Published: Elsevier B.V 15.06.2016
Subjects:
Algorithms
Approximation
Architecture
Architecture (computers)
Finite element method
First order approximation
Graphics processors
Mathematical analysis
Mathematical models
Multi-core processors
Multi-threaded programming
Numerical integration
Finite element method
First order approximation
Numerical integration
Multi-threaded programming
Graphics processors
Multi-core processors
ISSN:0045-7825, 1879-2138
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The paper presents investigations on the performance of the finite element numerical integration algorithm for first order approximations and three processor architectures, popular in scientific computing, classical x86_64 CPU, Intel Xeon Phi and NVIDIA Kepler GPU. We base the discussion on theoretical performance models and our own implementations for which we perform a range of computational experiments. For the latter, we consider a unifying programming model and portable OpenCL implementation for all architectures. Variations of the algorithm due to different problems solved and different element types are investigated and several optimizations aimed at proper optimization and mapping of the algorithm to computer architectures are demonstrated. The experimental results show the varying levels of performance for different architectures, but indicate that the algorithm can be effectively ported to all of them. The conclusions indicate the factors that limit the performance for different problems and types of approximation and the performance ranges that can be expected for FEM numerical integration on different processor architectures.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 23
ISSN:0045-7825
1879-2138
DOI:10.1016/j.cma.2016.03.038