I'm about to implement the multigrid preconditioner for my FEM code for the case of regulare 2D triangular meshes. My meshes have a union jack pattern, thus the nodes are equidistributed in both directions but the stencil varies between a 9-point and a 5-point stencil. Is it possible to fit this into the DMDA framework?
Thomas
