EN FR
EN FR


Section: New Results

High performance numerical computing

High order HDG schemes and domain decomposition solvers for frequency-domain electromagnetics

Participants : Emmanuel Agullo [HIEPACS project-team, Inria Bordeaux - Sud-Ouest] , Cristobal Samaniego Alvarado [HIEPACS project-team, Inria Bordeaux - Sud-Ouest] , Mathieu Faverge [HIEPACS project-team, Inria Bordeaux - Sud-Ouest] , Luc Giraud [HIEPACS project-team, Inria Bordeaux - Sud-Ouest] , Matthieu Kuhn [HIEPACS project-team, Inria Bordeaux - Sud-Ouest] , Stéphane Lanteri, Grégoire Pichon [HIEPACS project-team, Inria Bordeaux - Sud-Ouest] , Pierre Ramet [HIEPACS project-team, Inria Bordeaux - Sud-Ouest] .

This work is undertaken in the context of PRACE 5IP (http://www.prace-ri.eu/prace-5ip/) project and aims at the development of scalable frequency-domain electromagnetic wave propagation solvers, in the framework of the HORSE simulation software. HORSE is based on a high order HDG scheme formulated on an unstructured tetrahedral grid for the discretization of the system of three-dimensional Maxwell equations in heterogeneous media, leading to the formulation of large sparse undefinite linear system for the hybrid variable unknowns. This system is solved with domain decomposition strategies that can be either a purely algebraic algorithm working at the matrix operator level (i.e. a black-box solver), or a tailored algorithm designed at the continuous PDE level (i.e. a PDE-based solver). In the former case, we collaborate with the HIEPACS project-team at Inria Bordeaux - Sud-Ouest in view of adapting and exploiting the MaPHyS (Massively Parallel Hybrid Solver - https://gitlab.inria.fr/solverstack/maphys) algebraic hybrid iterative-direct domain decomposition solver. More precisely, this collaboration is concerned with two topics: one one hand, the improvement of the iterative convergence of MaPHyS for the HDG hybrid variable linear system and, on the other hand, the leveraging of low rank compression techniques for reducing the memory footprint of the factorization of subdomain problems using the PaStiX (Parallel Sparse matriX package - http://pastix.gforge.inria.fr/) package.

Figure 8. Scattering of a plane wave by a squadron Lockheed F-104 Starfighter. Contour lines of the amplitude of 𝐄 field. Simulations are performed with a HDG scheme based on a cubic interpolation of the electric and magnetic field unknowns, combined with a PDE-based domain decomposition solver.
IMG/F104Squad_P3_View1.png

High order HDG schemes and domain decomposition solvers for frequency-domain electromagnetics

Participants : Stéphane Lanteri, Laércio Lima Pilla [CORSE project-team, Inria Grenoble - Rhône Alpes] , Jean-François Méhaut [CORSE project-team, Inria Grenoble - Rhône Alpes] .

This work is undertaken in the context of PRACE 5IP (http://www.prace-ri.eu/prace-5ip/) project and aims at the development of a hybrid MPI/OpenMP parallellization of the DGTD solver of the DIOGENeS software suite. In practice, we concentrated our efforts on identifying and evaluating the best approaches for implementing fine grain parallism of the main DG numerical kernels, based on OpenMP features for loop-based parallelism on one hand, and task-based parallelism on the other hand.