Les méthodes hiérarchiques constituent aujourd'hui un outil incontournable pour accélérer les calculs sur de grandes matrices issues de simulations scientifiques. Elles permettent de réduire significativement les coûts en mémoire et en calcul en exploitant la structure de rang faible des interactions éloignées. Parmi les principales approches figurent les H-matrices (algébriques) et les méthodes analytiques de type Fast Multipole Method (FMM).
Bien que développées par des communautés différentes, ces approches reposent sur des principes communs : partitionnement hiérarchique, critères d'admissibilité, approximations de rang faible et organisation des bases. Les mettre en évidence ouvre la voie à des implémentations plus génériques, modulaires et réutilisables.
Dans ce contexte, un récent travail de recherche a proposé un cadre unifié et composable pour le produit matrice-vecteur hiérarchique :
Emmanuel Agullo, Bérenger Bramas, Olivier Coulaud, Antoine Gicquel. Unified and Composable Rank-Structured Hierarchical Matrix–Vector Multiplication. RR-9611, Inria, 2026.
Ce travail introduit une abstraction commune permettant d'exprimer aussi bien les H-matrices que les méthodes FMM au sein d'un même algorithme générique. Cette approche a été implémentée dans la bibliothèque Composyx, développée à Inria.
Les développements seront réalisés dans la bibliothèque Composyx, une bibliothèque C++ développée à Inria permettant l'expression modulaire d'algorithmes numériques.
