Iterative Methods for Sparse Linear Systems on Graphics Processing Unit

被引:14
作者
Ahamed, Abal-Kassim Cheik [1 ]
Magoules, Frederic [1 ]
机构
[1] Ecole Cent Paris, Appl Math & Syst Lab, Chatenay Malabry, France
来源
2012 IEEE 14TH INTERNATIONAL CONFERENCE ON HIGH PERFORMANCE COMPUTING AND COMMUNICATIONS & 2012 IEEE 9TH INTERNATIONAL CONFERENCE ON EMBEDDED SOFTWARE AND SYSTEMS (HPCC-ICESS) | 2012年
关键词
Krylov methods; linear algebra; sparse matrix-vector multiplication; graphics processing unit; CUDA; CUBLAS; CUSPARSE; Cusp; DOMAIN DECOMPOSITION METHODS;
D O I
10.1109/HPCC.2012.118
中图分类号
TP301 [理论、方法];
学科分类号
081202 ;
摘要
Many engineering and science problems require a computational effort to solve large sparse linear systems. Krylov subspace based iterative solvers have been widely used in that direction. Iterative Krylov methods involve linear algebra operations such as summation of vectors, dot product, norm, and matrix-vector multiplication. Since these operations could be very costly in computation time on Central Processing Unit (CPU), we propose in this paper to focus on the design of iterative solvers to take advantage of massive parallelism of Graphics Processing Unit (GPU). We consider Stabilized BiConjugate Gradient (BiCGStab), Stabilized BiConjugate Gradient (L) (BiCGStabl), Generalized Conjugate Residual (P-GCR), Bi-Conjugate Gradient Conjugate Residual (P-BiCGCR), transpose-free Quasi Minimal Residual (P-tfQMR) for the solution of sparse linear systems with non symmetric matrices and Conjugate Gradient (CG) for symmetric positive definite matrices. We discuss data format and data structure for sparse matrices, and how to efficiently implement these solvers on the Nvidia's CUDA platform. The scalability and performance of the methods are tested on several engineering problems, together with numerous numerical experiments which clearly illustrate the robustness, competitiveness and efficiency of our own proper implementation compared to the existing libraries.
引用
收藏
页码:836 / 842
页数:7
相关论文
共 25 条
[1]  
Aliaga J.I., 2012, LECT NOTES COMPUTER, V7133, P162
[2]  
[Anonymous], BIT, DOI DOI 10.1007/BF01933580
[3]  
[Anonymous], 2008, NVIDIA Technical Report NVR-2008-004
[4]  
[Anonymous], NVID CUDA TOOLK REF
[5]   Analysis of a conjugated infinite element method for acoustic scattering [J].
Autrique, Jean-Christophe ;
Magoules, Frederic .
COMPUTERS & STRUCTURES, 2007, 85 (09) :518-525
[6]  
Bell N, 2009, STUDENTS GUIDE TO THE MA TESOL, P1
[7]  
Bell N., 2012, CUSP: Generic parallel algorithms for sparse matrix and graph computations
[8]  
BERRY M. W., 1996, LECT APPL MATH, V32, P99
[9]  
CIARLET P. G., 2002, Classics in Appl. Math., V40
[10]   The University of Florida Sparse Matrix Collection [J].
Davis, Timothy A. ;
Hu, Yifan .
ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 2011, 38 (01)