<div dir="ltr">Thanks, Mark and Jed. It is very helpful. <div>So, for present, Petsc doesn't support Q1 interperpolation for cell-centered multigrid, right?</div><div><div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Xiaodong Liu, PhD<br>X: Computational Physics Division<br>Los Alamos National Laboratory<br>P.O. Box 1663, <br>Los Alamos, NM 87544<br>505-709-0534<br></div></div></div></div><br></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Mar 24, 2020 at 8:22 AM Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>> writes:<br>
<br>
> Good question. It does look like there is Q1:<br>
><br>
> src/dm/impls/da/da.c:-  ctype - DMDA_Q1 and DMDA_Q0 are currently the only<br>
> supported forms<br>
><br>
> And in looking at a cell centered<br>
> example src/snes/examples/tutorials/ex20.c, it looks like only DMDA_Q1<br>
> works. I get an error when I set it to DMDA_Q0 (DMDA_Q1 is the default).<br>
> This is puzzling, Q0 is natural in cell centered.<br>
<br>
The comments in those examples are kinda wrong -- they never told the DM<br>
it was cell-centered so it uses a multigrid that isn't compatible with<br>
the boundary conditions.  The interpolation is Q1 on the dual grid, not<br>
conservative Q1 on cells.<br>
<br>
> I am not familiar with DMDA and I don't understand why, from ex20, that you<br>
> have an odd number of points on a cell centered grid and an even number for<br>
> vertex centered (eg, ex14). I would think that it should be the opposite.<br>
<br>
The example is bad.<br>
</blockquote></div>