Méthodes de décomposition de domaine. Application au calcul haute performance

English

Séminaire Modèles et Algorithmes Déterministes: EDP-MOISE-MGMI

4/09/2014 - 11:00 Mr Pierre Jolivet (Laboratoire Jacques-Louis Lions) Salle 1 - Tour IRMA

Cette présentation introduit une vision unifiée de plusieurs méthodes de décomposition de domaine : celles avec recouvrement, dites de Schwarz, et celles basées sur des compléments de Schur, dites de sous-structuration. Il est ainsi possible de changer de méthodes de manière abstraite et de construire différents préconditionneurs pour accélérer la résolution de grands systèmes linéaires creux par des méthodes itératives. On rencontre régulièrement ce type de systèmes dans des problèmes industriels ou scientifiques après discrétisation de modèles continus. Bien que de tels préconditionneurs exposent naturellement de bonnes propriétés de parallélisme sur les architectures distribuées, ils peuvent s'avérer être peu performants numériquement pour des décompositions complexes ou des problèmes physiques multi-échelles. On peut pallier ces défauts de robustesse en calculant de façon concurrente des problèmes locaux creux ou denses aux valeurs propres généralisées. D'aucuns peuvent alors identifier des modes qui perturbent la convergence des méthodes itératives sous-jacentes a priori. En utilisant ces modes, il est alors possible de définir des opérateurs de projection qui utilisent un problème dit grossier. L'utilisation de ces outils auxiliaires règle généralement les problèmes sus-cités, mais tend à diminuer les performances algorithmiques des préconditionneurs. Dans cette présentation, on montre en trois points que la nouvelle construction développée est performante : 1) grâce à des essais numériques à très grande échelle sur Curie-un supercalculateur européen, puis en le comparant à des solveurs de pointe 2) multi-grilles et 3) directs.