PCPATCH: software for the topological construction of multigrid relaxation methods

Abstract

Effective relaxation methods are necessary for good multigrid convergence. For many equations, standard Jacobi and Gauß-Seidel are inadequate, and more sophisticated space decompositions are required; examples include problems with semidefinite terms or saddle point structure. In this paper we present a unifying software abstraction, PCPATCH, for the topological construction of space decompositions for multigrid relaxation methods. Space decompositions are specified by collecting topological entities in a mesh (such as all vertices or faces) and applying a construction rule (such as taking all degrees of freedom in the cells around each entity). The software is implemented in PETSc and facilitates the elegant expression of a wide range of schemes merely by varying solver options at runtime. In turn, this allows for the very rapid development of fast solvers for difficult problems.

22 pages, minor fixes in bibliography

References in corpus (3)

• PCPATCH: software for the topological construction of multigrid relaxation methods
• Fully Parallel Mesh I/O using PETSc DMPlex with an Application to Waveform Modeling
• Flexible, Scalable Mesh and Data Management using PETSc DMPlex

Cited by in corpus (15)

• A Reynolds-robust preconditioner for the Scott-Vogelius discretization of the stationary incompressible Navier-Stokes equations
• PCPATCH: software for the topological construction of multigrid relaxation methods
• Low-order divergence-free approximations for the Stokes problem on Worsey-Farin and Powell-Sabin splits
• An augmented Lagrangian preconditioner for the magnetohydrodynamics equations at high Reynolds and coupling numbers
• Monolithic multigrid for implicit Runge-Kutta discretizations of incompressible fluid flow
• A scalable and robust vertex-star relaxation for high-order FEM
• Robust multigrid methods for nearly incompressible elasticity using macro elements
• Low-order preconditioning of the Stokes equations
• Smoothers with localized residual computations for geometric multigrid methods
• An augmented Lagrangian preconditioner for implicitly-constituted non-Newtonian incompressible flow
• Monolithic Multigrid for Magnetohydrodynamics
• asQ: parallel-in-time finite element simulations using ParaDiag for geoscientific models and beyond
• Augmented Lagrangian preconditioners for the Oseen-Frank model of nematic and cholesteric liquid crystals
• Code generation for productive portable scalable finite element simulation in Firedrake
• Finite element appoximation and augmented Lagrangian preconditioning for anisothermal implicitly-constituted non-Newtonian flow