Modeling of Three-Dimensional Magnetotelluric Fields by Vector Finite Element Method Using GPU

Abstract

SUMMARYThree-dimensional modeling of magnetotelluric fields using the vector finite element method is carried out. We use the vector differential equation for the single vector-potential and the system of vector-scalar equations for the coupled vector and scalar potentials. We use the edge-based vector basis functions for the vector potential and nodal basis functions for the scalar potential. We use the COCR solver for the solution of the large sparse system of linear equations with Jacobi and folded preconditioner. The CPU and GPU implementations of the developed methods are compared. We test our methods on complex realistic problem of three-dimensional magnetotelluric sounding providing in Russian Far East. The conducted experiments show the advantage of GPU using.