La résolution de grands systèmes linéaires creux, fréquents en simulation numérique d'équation aux dérivées partielles, peut être efficacement abordée en combinant les méthodes de décomposition de domaine et le calcul en précision mixte. Les méthodes de décomposition de domaine permettent de diviser le problème en sous-domaines traitables en parallèle, améliorant ainsi la convergence. Néanmoins, ces méthodes sont limité par la quantitée d'espace mémoire disponible sur chaque noeud de calcul. L’utilisation de la précision mixte — faible précision pour les étapes peu sensibles et forte précision pour les calculs critiques — permet de réduire le coût en temps et en mémoire tout en maintenant la précision globale. Cette approche, bien adaptée aux architectures modernes (GPU, calcul parallèle), offre un bon compromis entre performance et fiabilité, à condition de maîtriser les effets du conditionnement et de la stabilité numérique. Le travail consiste donc à comprendre ces phénoménes pour les méthodes de décomposition de domaines.