[petsc-users] Setting intepolation and restriction operators in DMMG

[domenico.borzacchiello at univ-st-etienne.fr]

Hi,

I'm still in the process of coding my Stokes solver (FV+MAC
Discretisation) with DMMG.
So far the I've been testing it with direct solvers to check that
everything was right for the functional and jacobian creation and it runs
fine both in parallel and sequential mode. I now need to implement an
iterative solver  that will most likely be multigrid, and since I am using
fully staggered arrangement I need to define my own grid transfer
operators because the interpolation/restriction will be different for u v
w p. How can I do this?

Would you suggest another solution strategy instead?

thank you
Domenico.