[petsc-users] Geometric multigrid with a cell-centered discretization

Buesing, Henrik hbuesing at eonerc.rwth-aachen.de
Tue Jan 9 07:59:47 CST 2018

Dear all,

I am using DMDACreate3d to generate my domain and a cell-centered two-point flux approximation as discretization. I use geometric Multigrid (-pc_type mg) with DMDASetInterpolationType(petsc_da, DMDA_Q0, petsc_ierr). With DM_BOUNDARY_GHOSTED as boundary type in DMDACreate3d, I get

[0]PETSC ERROR: Arguments are incompatible
[0]PETSC ERROR: Ratio between levels: (mx - 1)/(Mx - 1) must be integer: mx 1536 Mx 768

With DM_BOUNDARY_PERIODIC, I get invalid accesses when I assemble the matrix.

Can I somehow tell him, that I have a cell-centered grid, or do I have to resort to using N+1 in my dimensions?

Thank you!

Dipl.-Math. Henrik Büsing
Institute for Applied Geophysics and Geothermal Energy
E.ON Energy Research Center
RWTH Aachen University
Mathieustr. 10            |    Tel +49 (0)241 80 49907
52074 Aachen, Germany     |    Fax +49 (0)241 80 49889
hbuesing at eonerc.rwth-aachen.de

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180109/c3aeac58/attachment.html>

More information about the petsc-users mailing list