Efficient implementation of Jacobi iterative method for large sparse linear systems on graphic processing units

Abstract : In this paper, an original Jacobi implementation is considered for the solution of sparse linear systems of equations. The proposed algorithm helps to optimize the parallel implementation on GPU. The performance analysis of GPU-based (using CUDA) algorithm of the implementation of this algorithm is compared to the corresponding serial CPU-based algorithm. Numerical experiments performed on a set of matrices arising from the finite element discretization of various equations (3D Laplace equation, 3D gravitational potential equation, 3D Heat equation) with different meshes, illustrate the performance, robustness and efficiency of our algorithm, with a speed up to 23\times in double-precision arithmetics.
Complete list of metadatas

https://hal-centralesupelec.archives-ouvertes.fr/hal-01708749
Contributor : Frédéric Magoulès <>
Submitted on : Wednesday, February 14, 2018 - 9:36:13 AM
Last modification on : Thursday, March 28, 2019 - 11:24:10 AM

Identifiers

Citation

Abal-Kassim Cheik Ahamed, Frédéric Magoulès. Efficient implementation of Jacobi iterative method for large sparse linear systems on graphic processing units. Journal of Supercomputing, Springer Verlag, 2017, 73 (8), pp.3411-3432. ⟨10.1007/s11227-016-1701-3⟩. ⟨hal-01708749⟩

Share

Metrics

Record views

126