MPI+X: task-based parallelisation and dynamic load balance of finite element assembly. Issue 3 (16th March 2019)

• Record Type: Journal Article
• Authors: Garcia-Gasulla, Marta; Houzeaux, Guillaume; Ferrer, Roger; Artigues, Antoni; López, Victor; Labarta, Jesús; Vázquez, Mariano
• Is Part Of: International journal of computational fluid dynamics. Volume 33:Issue 3(2019)
• Page Start: 115
• Page End: 136
• Publication Date: 2019-03-16
• Subjects: CFD -- finite element -- MPI+X -- MPI -- OpenMP -- dynamic load balance -- shared-memory parallelism -- distributed memory parallelism -- hybrid parallelism; Fluid dynamics -- Data processing -- Periodicals; 532.05
• DOI: 10.1080/10618562.2019.1617856
• Languages: English
• ISSNs: 1061-8562

Abstract

ABSTRACT: The main computing phases of numerical methods for solving partial differential equations are the algebraic system assembly and the iterative solver. This work focuses on the first task, in the context of a hybrid MPI+X paradigm. The matrix assembly consists of a loop over the elements, faces, edges or nodes of the MPI partitions to compute element matrices and vectors and then of their assemblies. In a MPI+X hybrid parallelism context, X has consisted traditionally of loop parallelism using OpenMP, with different techniques to avoid the race condition, but presenting efficiency or implementation drawbacks. We propose an alternative, based on task parallelism using some extensions to the OpenMP programming model. In addition, dynamic load balance will be applied, especially efficient in the presence of hybrid meshes. This paper presents the proposed methodology, its implementation and its validation through the solution of large computational mechanics problems up to 16k cores.