HPDDM and PETSc 3.25 for CPU-GPU Architectures
Robust Schwarz and GenEO-type preconditioners available in CPU-GPU PETSc workflows
HPDDM domain-decomposition preconditioners are now available in modern PETSc workflows, including CPU-GPU contexts, through developments integrated in PETSc 3.25.
1. Context
Domain decomposition methods are central to solving the very large linear systems that arise in scientific and industrial simulations. HPDDM provides advanced solvers and preconditioners, including multilevel Schwarz methods and GenEO-type coarse spaces.
The 2026 work focused on making these robust methods directly usable inside the modern PETSc ecosystem, with particular attention to heterogeneous CPU-GPU architectures.
2. Key Result
HPDDM internal structures were revised to rely more directly on PETSc-native objects and operations. The result is more than a technical GPU port: advanced domain-decomposition preconditioners can now be used in production PETSc workflows while benefiting from PETSc accelerator abstractions.
3. Impact
This removes a major software barrier for NumPEx and Exa-MA by bringing advanced preconditioning methods closer to the heterogeneous architectures targeted by exascale computing.
Feel++ has already integrated PETSc 3.25 and can now exploit HPDDM strategies through PETSc in real application cases. Other Exa-MA frameworks using PETSc, such as FreeFEM or TRUST/CEA, should be able to benefit from these developments with reduced integration cost.
4. Next Steps
-
Validate the methods at larger scale.
-
Optimize performance on NVIDIA and AMD GPUs.
-
Extend support toward Kokkos.
-
Evaluate mixed and multiphysics problems representative of Exa-MA applications.