<div dir="ltr"><div><div>Thanks for your quick answers. To be honest, I am not familiar at all with DMShells and DMPlexes. But since it is what I need, I am going to try it. <br></div>Thanks again  for your advices,<br></div>Natacha <br> </div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 21, 2017 at 2:27 PM, Lawrence Mitchell <span dir="ltr"><<a href="mailto:lawrence.mitchell@imperial.ac.uk" target="_blank">lawrence.mitchell@imperial.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class=""><br>
> On 21 Mar 2017, at 13:24, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
><br>
> I think the remedy is as easy as specifying a DMShell that has a PetscSection (DMSetDefaultSection) with your ordering, and<br>
> I think this is how Firedrake (<a href="http://www.firedrakeproject.org/" rel="noreferrer" target="_blank">http://www.firedrakeproject.<wbr>org/</a>) does it.<br>
<br>
</span>We actually don't use a section, but we do provide DMCreateFieldDecomposition_<wbr>Shell.<br>
<br>
If you have a section that describes all the fields, then I think if the DMShell knows about it, you effectively get the same behaviour as DMPlex (which does the decomposition in the same manner?).<br>
<span class=""><br>
> However, I usually use a DMPlex which knows about my<br>
> mesh, so I am not sure if this strategy has any holes.<br>
<br>
</span>I haven't noticed anything yet.<br>
<span class="HOEnZb"><font color="#888888"><br>
Lawrence</font></span></blockquote></div><br></div>