Improving edge finite element assembly for geophysical electromagnetic modelling on shared-memory architectures
Document typeConference report
Rights accessOpen Access
European Commission's projectGEAGAM - Geophysical Exploration using Advanced GAlerkin Methods (EC-H2020-644202)
HPC4E - HPC for Energy (EC-H2020-689772)
This work presents a set of node-level optimizations to perform the assembly of edge finite element matrices that arise in 3D geophysical electromagnetic modelling on shared-memory architectures. Firstly, we describe the traditional and sequential assembly approach. Secondly, we depict our vectorized and shared-memory strategy which does not require any low level instructions because it is based on an interpreted programming language, namely, Python. As a result, we obtained a simple parallel-vectorized algorithm whose runtime performance is considerably better than sequential version. The set of optimizations have been included to the work-flow of the Parallel Edge-based Tool for Geophysical Electromagnetic Modelling (PETGEM) which is developed as open-source at the Barcelona Supercomputing Center. Finally, we present numerical results for a set of tests in order to illustrate the performance of our strategy.
CitationCastillo-Reyes, Octavio; de la Puente, Josep; Cela, Jose M. Improving edge finite element assembly for geophysical electromagnetic modelling on shared-memory architectures. A: IEEE Annual Conference in Ubiquitous Computing, Electronics & Mobile Communication Conference (UEMCON), 20-22 Oct. 2016. "Ubiquitous Computing, Electronics & Mobile Communication Conference (UEMCON), IEEE Annual". IEE, 2016, p. 1-6.