This paper presents an edge-based parallel code for the data computation that arises when applying one of the most popular electromagnetic methods in geophysics, namely, the controlled-source electromagnetic method (CSEM). The computational implementation is based on the linear Edge Finite Element Method in 3D isotropic domains because it has the ability to eliminate spurious solutions and is claimed to yield accurate results. The framework structure is able to exploit the embarrassingly-parallel tasks and the advantages of the geometric flexibility as well as to work with three different orientations for the dipole, or excitation source, on unstructured tetrahedral meshes in order to represent complex geological bodies through a local refinement technique. We demonstrate the performance and accuracy of our tool on the Marenostrum supercomputer (Barcelona Supercomputing Center) through scaling tests and canonical tests, respectively.

This project received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 644202. Authors gratefully acknowledge the support from the Mexican National Council for Science and Technology (CONACyT). This work benefitted from helpful suggestions of Hélène Barucq and Julien Diaz from MAGIQUE-3D team at Institut National de Recherche en Informatique et en Automatique (INRIA) and two anonymous reviewers.

Peer Reviewed