Acceleration of finite element method for 3D DC resistivity modeling using multi-GPU

Abstract

In this paper finite element method for 3D DC resistivity modeling accelerated using multi-GPU (Graphics Processing Unit). Solution of the large system of linear equations is the most expensive computation in finite element method performed in GPUs to reduce the computational time. Conjugate gradient solver used to solve large system of linear equations. We developed kernel for conjugate gradient solver that exploit data vectorization and written in PTX assembly form. We perform test on GTX 750Ti GPU and Tesla C2050 GPU. Our kernel have better sparse matrix-vector performance than CUSPARSE library in the first GPU, but lower performance in the latter GPU. The performance comparisons to the library are about 1.4 times and 0.7 times respectively. Our multi-GPU implementation achieved about 1.9 times performance of single GPU by using 2 identical GPUs. In comparison to the serial CPU implementation, about 10 times speedup could be achieved by using 2 GPUs.