<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jan 9, 2018 at 10:12 AM, Buesing, Henrik <span dir="ltr"><<a href="mailto:hbuesing@eonerc.rwth-aachen.de" target="_blank">hbuesing@eonerc.rwth-aachen.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Sorry (again), for posting only part of the error message...<br>
<br>
So, here it comes. I will first [1] show the code snippets I am using, then second [2] I will show the PETSc options, and finally [3] the error message I get.<br>
The grid is 1536 x 1 x 256 in x,y,z-direction. Thus, we see the first level of coarsening in x-direction here.<br>
<br>
Thank you!<br>
Henrik<br>
<br>
[1]<br>
CALL DMDACreate3d(PETSC_COMM_WORLD,<wbr>DM_BOUNDARY_GHOSTED,DM_<wbr>BOUNDARY_GHOSTED,DM_BOUNDARY_<wbr>GHOSTED,&<br>
DMDA_STENCIL_STAR, I0g,J0g,K0g, PETSC_DECIDE, PETSC_DECIDE,PETSC_DECIDE, petsc_two, petsc_one, &<br>
PETSC_NULL_INTEGER,PETSC_NULL_<wbr>INTEGER,PETSC_NULL_INTEGER, petsc_da, petsc_ierr)<br>
CALL DMSetUp(petsc_da,petsc_ierr)<br>
...<br>
CALL DMSetFromOptions(petsc_da, petsc_ierr)<br>
CALL DMDASetInterpolationType(<wbr>petsc_da, DMDA_Q0, petsc_ierr)<br>
CALL SNESSetDM(petsc_snes, petsc_da, petsc_ierr)<br>
CALL SNESSetFromOptions(petsc_snes,<wbr>petsc_ierr)<br>
<br>
[2] -snes_monitor -snes_linesearch_monitor -snes_type newtonls -snes_linesearch_type basic -ksp_type fgmres -ksp_max_it 1000 -snes_stol 1e-06 -snes_rtol 1e-06 -snes_atol 1e-6 -ksp_gmres_restart 100 -ksp_monitor -pc_type mg -da_refine_y 1 -ksp_converged_reason -snes_converged_reason -ksp_monitor -ksp_view -pc_mg_levels 2 -options_left<br>
<br>
[3]<br>
[0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--<br>
<span class="">[0]PETSC ERROR: Arguments are incompatible<br>
[0]PETSC ERROR: Ratio between levels: (mx - 1)/(Mx - 1) must be integer: mx 1536 Mx 768<br></span></blockquote><div><br></div><div>It looks like you give it the number of vertices, not cells, so +1 to your grid size.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">
</span>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" rel="noreferrer" target="_blank">http://www.mcs.anl.gov/petsc/<wbr>documentation/faq.html</a> for trouble shooting.<br>
[0]PETSC ERROR: Petsc Development GIT revision: v3.8.2-48-g851ec02 GIT Date: 2017-12-05 09:52:17 -0600<br>
[0]PETSC ERROR: /rwthfs/rz/cluster/work/<wbr>hb111949/EoCoE/aixcelerate17/<wbr>gge_gnu/bench_run/000023/<wbr>000000_copy_sub_launch/work/<wbr>shem_fw64gnu_const.x on a gnu_openmpi named <a href="http://linuxbmc0008.rz.RWTH-Aachen.DE" rel="noreferrer" target="_blank">linuxbmc0008.rz.RWTH-Aachen.DE</a> by hb111949 Tue Jan 9 16:04:39 2018<br>
[0]PETSC ERROR: Configure options --download-fblaslapack --with-cc=mpicc -with-fc=mpif90 --with-cxx=mpicxx --download-hypre --download-superlu_dist --download-suitesparse --download-scalapack --download-blacs --download-hdf5 --download-parmetis --download-metis --with-debugging=0 --download-mumps<br>
[0]PETSC ERROR: #1 DMCreateInjection_DA_3D() line 1195 in /rwthfs/rz/cluster/work/<wbr>hb111949/Code/petsc/src/dm/<wbr>impls/da/dainterp.c<br>
[0]PETSC ERROR: #2 DMCreateInjection_DA() line 1283 in /rwthfs/rz/cluster/work/<wbr>hb111949/Code/petsc/src/dm/<wbr>impls/da/dainterp.c<br>
[0]PETSC ERROR: #3 DMCreateInjection() line 1082 in /rwthfs/rz/cluster/work/<wbr>hb111949/Code/petsc/src/dm/<wbr>interface/dm.c<br>
[0]PETSC ERROR: #4 DMCoarsen_DA() line 1219 in /rwthfs/rz/cluster/work/<wbr>hb111949/Code/petsc/src/dm/<wbr>impls/da/da.c<br>
[0]PETSC ERROR: #5 DMCoarsen() line 2427 in /rwthfs/rz/cluster/work/<wbr>hb111949/Code/petsc/src/dm/<wbr>interface/dm.c<br>
[0]PETSC ERROR: #6 PCSetUp_MG() line 618 in /rwthfs/rz/cluster/work/<wbr>hb111949/Code/petsc/src/ksp/<wbr>pc/impls/mg/mg.c<br>
[0]PETSC ERROR: #7 PCSetUp() line 924 in /rwthfs/rz/cluster/work/<wbr>hb111949/Code/petsc/src/ksp/<wbr>pc/interface/precon.c<br>
[0]PETSC ERROR: #8 KSPSetUp() line 381 in /rwthfs/rz/cluster/work/<wbr>hb111949/Code/petsc/src/ksp/<wbr>ksp/interface/itfunc.c<br>
[0]PETSC ERROR: #9 KSPSolve() line 612 in /rwthfs/rz/cluster/work/<wbr>hb111949/Code/petsc/src/ksp/<wbr>ksp/interface/itfunc.c<br>
[0]PETSC ERROR: [1]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------------------------------<wbr>--<br>
<span class=""><br>
<br>
--<br>
Dipl.-Math. Henrik Büsing<br>
Institute for Applied Geophysics and Geothermal Energy<br>
E.ON Energy Research Center<br>
RWTH Aachen University<br>
------------------------------<wbr>------------------------<br>
Mathieustr. 10 | Tel <a href="tel:%2B49%20%280%29241%2080%2049907" value="+492418049907">+49 (0)241 80 49907</a><br>
52074 Aachen, Germany | Fax <a href="tel:%2B49%20%280%29241%2080%2049889" value="+492418049889">+49 (0)241 80 49889</a><br>
------------------------------<wbr>------------------------<br>
<a href="http://www.eonerc.rwth-aachen.de/GGE" rel="noreferrer" target="_blank">http://www.eonerc.rwth-aachen.<wbr>de/GGE</a><br>
<a href="mailto:hbuesing@eonerc.rwth-aachen.de">hbuesing@eonerc.rwth-aachen.de</a><br>
------------------------------<wbr>------------------------<br>
<br>
</span>> -----Ursprüngliche Nachricht-----<br>
> Von: Smith, Barry F. [mailto:<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>]<br>
> Gesendet: 09 January 2018 15:10<br>
> An: Buesing, Henrik <<a href="mailto:hbuesing@eonerc.rwth-aachen.de">hbuesing@eonerc.rwth-aachen.<wbr>de</a>><br>
> Cc: petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>><br>
> Betreff: Re: [petsc-users] Geometric multigrid with a cell-centered<br>
> discretization<br>
<div class="HOEnZb"><div class="h5">><br>
><br>
> With cell centered elements each element is sliced in half in each of<br>
> the three directions for refinement (or you can refine only in 1 or 2 of<br>
> the dimensions if you like). This means that the ratio of mx to MX needs<br>
> to be 2 (or 4).So your numbers are fine.<br>
><br>
> > On Jan 9, 2018, at 7:59 AM, Buesing, Henrik <hbuesing@eonerc.rwth-<br>
> <a href="http://aachen.de" rel="noreferrer" target="_blank">aachen.de</a>> wrote:<br>
> ><br>
> > Dear all,<br>
> ><br>
> > I am using DMDACreate3d to generate my domain and a cell-centered<br>
> > two-point flux approximation as discretization. I use geometric<br>
> > Multigrid (-pc_type mg) with DMDASetInterpolationType(<wbr>petsc_da,<br>
> > DMDA_Q0, petsc_ierr). With DM_BOUNDARY_GHOSTED as boundary type in<br>
> > DMDACreate3d, I get<br>
> ><br>
> > [0]PETSC ERROR: Arguments are incompatible [0]PETSC ERROR: Ratio<br>
> > between levels: (mx - 1)/(Mx - 1) must be integer: mx 1536 Mx 768<br>
><br>
> We need ALL the output from the error (always). It should never be<br>
> doing this check in your case so something is wrong with the order of<br>
> your code perhaps.<br>
><br>
> Barry<br>
><br>
> ><br>
> > With DM_BOUNDARY_PERIODIC, I get invalid accesses when I assemble the<br>
> matrix.<br>
> ><br>
> > Can I somehow tell him, that I have a cell-centered grid, or do I have<br>
> to resort to using N+1 in my dimensions?<br>
> ><br>
> > Thank you!<br>
> > Henrik<br>
> ><br>
> ><br>
> > --<br>
> > Dipl.-Math. Henrik Büsing<br>
> > Institute for Applied Geophysics and Geothermal Energy E.ON Energy<br>
> > Research Center RWTH Aachen University<br>
> > ------------------------------<wbr>------------------------<br>
> > Mathieustr. 10 | Tel +49 (0)241 80 49907<br>
> > 52074 Aachen, Germany | Fax +49 (0)241 80 49889<br>
> > ------------------------------<wbr>------------------------<br>
> > <a href="http://www.eonerc.rwth-aachen.de/GGE" rel="noreferrer" target="_blank">http://www.eonerc.rwth-aachen.<wbr>de/GGE</a><br>
> > <a href="mailto:hbuesing@eonerc.rwth-aachen.de">hbuesing@eonerc.rwth-aachen.de</a><br>
> > ------------------------------<wbr>------------------------<br>
<br>
</div></div></blockquote></div><br></div></div>