Premium
An efficient sparse matrix‐vector multiplication on CUDA‐enabled graphic processing units for finite element method simulations
Author(s) -
Altinkaynak Atakan
Publication year - 2016
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.5346
Subject(s) - kernel (algebra) , cuda , discretization , finite element method , sparse matrix , multiplication (music) , solver , computer science , parallel computing , node (physics) , computational science , matrix (chemical analysis) , floating point , algorithm , matrix multiplication , mathematics , discrete mathematics , mathematical analysis , physics , materials science , quantum mechanics , combinatorics , gaussian , composite material , quantum , programming language , thermodynamics
Summary Finite element method (FEM) is a well‐developed method to solve real‐world problems that can be modeled with differential equations. As the available computational power increases, complex and large‐size problems can be solved using FEM, which typically involves multiple degrees of freedom (DOF) per node, high order of elements, and an iterative solver requiring several sparse matrix‐vector multiplication operations. In this work, a new storage scheme is proposed for sparse matrices arising from FEM simulations with multiple DOF per node. A sparse matrix‐vector multiplication kernel and its variants using the proposed scheme are also given for CUDA‐enabled GPUs. The proposed scheme and the kernels rely on the mesh connectivity data from FEM discretization and the number of DOF per node. The proposed kernel performance was evaluated on seven test matrices for double‐precision floating point operations. The performance analysis showed that the proposed GPU kernel outperforms the ELLPACK (ELL) and CUSPARSE Hybrid (HYB) format GPU kernels by an average of 42% and 32%, respectively, on a Tesla K20c card. Copyright © 2016 John Wiley & Sons, Ltd.