GPU-based preconditioned conjugate gradient method for solving sparse linear systems
Shen Defei · Journal of Computer Applications · 2013
A GPU-accelerated preconditoned conjugate gradient method was studied to solve sparse linear equations.And the sparse matrix was stored in the Compressed Sparse Row(CSR) format.The programmes were coded on Compute Unified Device Architecture(CUDA) and tested on the device of nVidia GT430 GPU.According to the features of conjugate gradient method,strategies were investigated to optimize the sparse matrix vector multiplication and the data transfer between CPU and GPU.Compared with the implementation calling cusparseDcsrmv,the self-developed kernel code of sparse matrix vector multiplication can go to a speed-up of 2.1 in the best case.Equipped with this kernel,the preconditioned conjugate gradient code obtains a maximum speed-up of 7.4 against the CPU code,which is a bit advantageous over that using CUBLAS library and CUSPARSE library.