2-D electromagnetic modelling by finite element method on GPU

2-D electromagnetic model with a line source uses Preconditioned Conjugate Gradient to solve the stiffness matrix and Herrmann's pseudo-delta function to simulate the neighbor-hood of the singular source. We employ this Finite Element Method (FEM) on Fermi-GPU, which is the third generation GPU for the general computing purpose. Our algorithm takes advantage of the GPU's particular characteristics of Single Instruction Multiple Thread (SIMT) architecture and runs entirely on the GPU to minimize the transfer time between CPU and GPU. Meanwhile, we design a reasonable pre-conditioner (Diagonal Pre-conditioner) on the GPU and inject a mixed scheme into our algorithm for obtaining the good performance. In the numerical experiments, the correctness of the algorithm is acknowledged, and it achieves high performance by outperforming CPU solvers over hundreds of times. Indeed, it allows finer mesh generation so as to obtain more accurate resolutions and advance the ability of data interpretation of geophysical prospecting. (C) 2016 Elsevier GmbH. All rights reserved.