<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Jun 2, 2017 at 5:04 PM, Alberto Paganini <span dir="ltr"><<a href="mailto:Alberto.Paganini@maths.ox.ac.uk" target="_blank">Alberto.Paganini@maths.ox.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear PETSc developers,<br>
<br>
I'm Alberto and I'm a user of the finite element library Firedrake,<br>
which relies on DMPlex to import meshes.<br></blockquote><div><br></div><div>Great.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
In order to use higher-order FEs, it is desirable to import higher-order meshes.<br></blockquote><div><br></div><div>I really do not like that term. Let me try and convince you that it is wrong. The topology of the</div><div>mesh is unchanged. You are only talking about the order of the representation of the geometry</div><div>field. Thus, it is not the mesh that is "higher order", but the geometry.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I've been told that DMPlex does not offer this future (at present).<br></blockquote><div><br></div><div>Toby just merged this to master, so I think we can say that we have alpha support for this. How</div><div>does it work? We already have a coordinateDM and coordinates Vec, so you just choose a</div><div>higher order discretization for the DS inside the coordinateDM. Does that make sense?</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I would be happy to contribute and add higher-order coordinate fields.<br>
Would someone be interested in giving me some support?<br>
<br>
I'm not the greatest programmer and I have (a limited) experience in C(++),<br></blockquote><div><br></div><div>Gentlemen do not write C++.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
but I'm motivated and I have a certain knowledge of FEs (which might be relevant).<br></blockquote><div><br></div><div>The most valuable thing for us is writing a good test. Preferably one in which an exact</div><div>solution can only be obtained if the high order geometry is being used correctly. If you</div><div>write that, we can make it work. Again, preferably in the style of existing PETSc tests.</div><div><br></div><div>  Thanks,</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">
Best<br>
<span class="HOEnZb"><font color="#888888">Alberto<br>
<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><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">http://www.caam.rice.edu/~mk51/</a><br></div></div></div>
</div></div>