<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Fri, Jan 4, 2019 at 10:04 AM Yann Jobic via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</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">Dear Petsc Users,<br>
<br>
I'm using DMPlexCreateFromCellList to create my DM. I would like to have <br>
an order 2 geometry.<br>
<br>
It's working fine in 2D for elements of type Q9.<br></blockquote><div><br></div><div>I do not see how that is possible. The topology does not work that way, and there</div><div>is no way a 9 vertex quad was interpolated, so something else is happening.</div><div><br></div><div>Moreover, anything you input still has affine geometry. Toby has been working on</div><div>non-affine geometry, and all the code is there, but right now you have to build and</div><div>populate the coordinate space by hand. What this means is that you create a</div><div>quadratic space for the coordinate DM (rather than the default linear), get the</div><div>global coordinate vector out, and stick in the Q9 coordinates. The topology is still</div><div>the same no matter what coordinates you read in.</div><div><br></div><div>We will get around to making an interface for this soon.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
I checked that it's working correctly by using <br>
DMPlexComputeCellGeometryFEM, and compute the value of the determinant <br>
(CheckMeshGeometry of dm/impls/plex/examples/tutorials/ex2.c)<br>
<br>
I can import H8 elements, it's working fine.<br>
<br>
But for H27 element, it's not working. I really don't know how to order <br>
my vertex in the cells array. So far, the determinant is zero...<br>
<br>
I don't know where to look in order to find this information. I tried <br>
the function DMPlexGetRawFaces_Internal of the file<br>
<br>
<a href="https://www.mcs.anl.gov/petsc/petsc-current/src/dm/impls/plex/plexinterpolate.c.html#DMPlexGetRawFaces_Internal" rel="noreferrer" target="_blank">https://www.mcs.anl.gov/petsc/petsc-current/src/dm/impls/plex/plexinterpolate.c.html#DMPlexGetRawFaces_Internal</a><br>
<br>
but it didn't help me.<br>
<br>
Could you please point me where to look ?<br>
<br>
Thanks, and happy new year !<br>
<br>
Yann<br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="m_-4020022596717474202gmail_signature"><div dir="ltr"><div><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.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>