
<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div class=""><br class=""></div>Hi Jerome,<div class=""><br class=""></div><div class="">Yes indeed - all our meshing algorithms currently generate 1st order geometrical elements, starting from 1st order geometrical elements on the boundary. The meshes are then "curved" later on. I will add the ability to re-use a given high-order discrete surface input when we curve to our TODO list. But it's not trivial.</div><div class=""><br class=""></div><div class="">Christophe<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 1 Jun 2018, at 14:56, Jerome Robert <<a href="mailto:jeromerobert@gmx.com" class="">jeromerobert@gmx.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">Hi all,<br class=""><br class="">I am trying to fill an existing order 2 3D triangular mesh with order 2 tetrahedron. In the result mesh, the order 2 vertices of my input mesh are all moved so they were linearly interpolated. So everything is as if my input mesh was used by gmsh as an order 1 surface.<br class=""><br class="">Here is my .geo:<br class=""><br class="">// Various options I tryed with no luck<br class="">//Mesh.SecondOrderExperimental=1;<br class="">//Mesh.ElementOrder = 2;<br class="">//Mesh.SecondOrderLinear = 0;<br class="">//Mesh.Algorithm3D = 1;<br class=""><br class="">// surf.unv contains only order 2 triangles<br class="">Merge "surf.unv";<br class="">Surface Loop(2) = {1};<br class="">Volume(3) = {2};<br class="">Characteristic Length {3} = 1000;<br class=""><br class="">Which I run with gmsh 3.0.6:<br class=""><br class="">$ gmsh -3 my.geo -o my.msh<br class=""><br class="">My understanding is that my input mesh is converted to order 1 by this SetOrder1 call:<br class=""><br class=""><a href="https://gitlab.onelab.info/gmsh/gmsh/blob/a4306e50a3916dd52f72edce56c05f36954f3f0a/Mesh/Generator.cpp#L979" class="">https://gitlab.onelab.info/gmsh/gmsh/blob/a4306e50a3916dd52f72edce56c05f36954f3f0a/Mesh/Generator.cpp#L979</a><br class=""><br class="">Interpolation on GEntity::DiscreteSurface or GEntity::DiscreteCurve is also always linear.<br class=""><br class="">So there seems to be no way to use an order 2 mesh as a gmsh GModel. <br class=""><br class="">Is that correct ? If not what would be the way to do it ?<br class=""><br class="">Regards,<br class=""><br class="">Jerome<br class=""><br class="">_______________________________________________<br class="">gmsh mailing list<br class="">gmsh@onelab.info<br class="">http://onelab.info/mailman/listinfo/gmsh<br class=""></div></div></blockquote></div><br class=""><div class="">

<div dir="auto" style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;">— </div><div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;">Prof. Christophe Geuzaine<br class="">University of Liege, Electrical Engineering and Computer Science <br class=""><a href="http://www.montefiore.ulg.ac.be/~geuzaine" class="">http://www.montefiore.ulg.ac.be/~geuzaine</a><br class=""><br class="">Free software: http://gmsh.info | http://getdp.info | http://onelab.info</div></div>

</div>

<br class=""></div></body></html>