<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Thu, Jul 12, 2018 at 10:00 AM Lawrence Mitchell <<a href="mailto:lawrence.mitchell@imperial.ac.uk">lawrence.mitchell@imperial.ac.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
> On 12 Jul 2018, at 12:27, Stefano Zampini <<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>> wrote:<br>
> <br>
> Lawrence,<br>
> <br>
> I wrote the periodic GMSH reader and I can answer for it.<br>
> <br>
> With DMPlex, when we have periodicity, we currently have two ways of encoding it<br>
> <br>
> 1) through the maxCell concept (or, as Lisandro's like to call it, the "maxCell crap")<br>
> 2) through the coordinatesSection<br>
> <br>
> Since periodicity with general meshed cannot be encoded with 1), I populate the cell part of the coordinates section with the "localized coordinates" (only for periodic cells), and leave all the other dofs associated with the vertices.<br>
> I believe this kind of localization is also what you get with the maxCell crap (check it out with the BoxMesh): even in this case in the coordinate section you have cell dofs only for cells touching the periodic boundary, and all the vertices have their own coordinates.<br>
> Look eg in DMPlexComputeLineGeometry_Internal.<br>
> <br>
> What you know for sure is that the first 6 entries of the closure are the ones you should look for (not the size of the closure)<br>
<br>
Thanks.  This seems somewhat odd to me.  At the periodic boundary, the coordinate field is not single-valued, so I don't understand how you can reasonably put coordinate dofs on the vertices.  I would have thought that on the periodic cells, the vertices in the closure have no dofs in the coordinate section, and it's just the cell that has dofs associated with it.<br>
<br>
Is there a reason why you leave things on the vertices at all?<br></blockquote><div><br></div><div>You can think of it being multi-values, or you can think of the vertex values as lying on a specific Riemann sheet.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Cheers,<br>
<br>
Lawrence</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div>