<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Jun 2, 2017 at 7:55 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
<br>
> On Fri, Jun 2, 2017 at 5:04 PM, Alberto Paganini <<br>
> <a href="mailto:Alberto.Paganini@maths.ox.ac.uk">Alberto.Paganini@maths.ox.ac.<wbr>uk</a>> wrote:<br>
><br>
>> 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>
>><br>
><br>
> Great.<br>
><br>
><br>
>> In order to use higher-order FEs, it is desirable to import higher-order<br>
>> meshes.<br>
>><br>
><br>
> I really do not like that term. Let me try and convince you that it is<br>
> wrong. The topology of the<br>
> mesh is unchanged. You are only talking about the order of the<br>
> representation of the geometry<br>
> field. Thus, it is not the mesh that is "higher order", but the geometry.<br>
><br>
><br>
>> I've been told that DMPlex does not offer this future (at present).<br>
>><br>
><br>
> Toby just merged this to master, so I think we can say that we have alpha<br>
> support for this. How<br>
> does it work? We already have a coordinateDM and coordinates Vec, so you<br>
> just choose a<br>
> higher order discretization for the DS inside the coordinateDM. Does that<br>
> make sense?<br>
<br>
</span>Can it load quadratic geometry from a file (ExodusII or otherwise)?<br></blockquote><div><br></div><div>If someone requests a given file format, we can do it. That's how we always proceed.</div><div> </div><div>   Matt</div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
That's the interesting part -- the quadratic geometry is rarely<br>
interesting except in elements that touch the boundary and you can't<br>
infer the true curved boundary without access to a geometric model (a<br>
much more cumbersome software stack).<br>
</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>