I don't think you want to use -pc_type gamg if you want to use geometric
multigrid. You can can use -pc_type mg and the DMDA. The only thing I think you
need to change from the default use of DMDA and -pc_type mg is to provide
custom code that computes the interpolation between levels to
Good morning,
Good afternoon,
I am testing the performance of some preconditioners on a problem I have
and I wanted to try geometric multigrid. I am using a DMDA object and
respective matrices (the mesh is structured, but not rectangular, it is
curvilinear)
I have used the default version of
