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.
