Hi! The "Possibilities for extensions" section of step-75 gives you a few hints on how to achieve p-multigrid preconditioning with matrix-based methods: https://www.dealii.org/current/doxygen/deal.II/step_75.html#Solvewithmatrixbasedmethods
If you look for an example implementation, the "hpbox" tool demonstrates both matrix-based and matrix-free implementations of the p-multigrid method. https://zenodo.org/record/6425947 https://github.com/marcfehling/hpbox/ Best, Marc On Monday, October 3, 2022 at 3:39:51 AM UTC-6 yy.wayne wrote: > Hello, > > I plan to test the converge of polynomial multigrid preconditioning on a > problem, with matrix-based framework. However, step-75 is a bit confusing > since I cannot clearly separate classses contributed by matrix-free and > p-multigrid. > > Specifically, in step-16 (with MeshWorker) and step-56 (assemble > cell_matrix by hand) there are assemble_system() and assemble_multigrid(), > while in step-75 there are solve_system() which calls solve_with_gmg and > mg_solve. It is due to matrix-free or p-multigrid? By the way, is there a > code gallery for matrix-based polynomial-multigrid ? > > Thanks > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/19802e08-cd43-49a9-a060-2a62ff68c7e6n%40googlegroups.com.