I recently saw this NGSolve tutorial on FETI-DP domain decomposition methods
5.2 FETI-DP in NGSolve I: Working with Point-Constraints — NGS-Py 6.2.1903 documentation.

I am trying to implement in PETSc the method, hence I was wondering whether the code uses PETSc in the background to handle the parallelization.

Thanks a lot in advance.

Hi, this example is using pure MPI NGSolve, but it is also possible to couple to petsc