Open schnellerhase opened 1 week ago
pmg
already has a quite versatile implementation of the multigrid operations available. Especially pmg.hpp uses very much the required interface of a fully abstract MG implementation.
The transfer operators are the matrices produced by fem.interpolation_matrix
- this should allow for easy replication of a $p$-multigrid as an example.
The pmg-dolfinx, implementing $p$-multigrid, project could be very interesting for the geometric multigrid implementation.