
<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"><base href="x-msg://1807/"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>Le 6 avr. 2013 à 21:57, John Pease Moore <<a href="mailto:jpmoore4@MIT.EDU">jpmoore4@MIT.EDU</a>> a écrit :</div><br class="Apple-interchange-newline"><blockquote type="cite"><div ocsi="0" fpstyle="1" style="font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="direction: ltr; font-family: Tahoma; font-size: 10pt; ">Hi All,<br><br>I'm trying to create a very simple high order curved mesh of a sphere. When I generate the mesh in gmsh and visualize the surface, everything looks fine. However, when I take the mesh, read it into MATLAB and plot it myself, it looks terrible--it is no longer high order. The polynomial order is 5.<br><br></div></div></blockquote><div><br></div><div>Are you sure you plot it correctly in matlab : in gmsh, we subdivide a lot high order elements for visualization….</div><div>There is indeed an option in Gmsh for controlling the subdivision : gmsh -numsubedges 10 refines every mesh edge </div><div>10 times (which represent about 100 subtriangles for one high order triangle).</div><br><blockquote type="cite"><div ocsi="0" fpstyle="1" style="font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="direction: ltr; font-family: Tahoma; font-size: 10pt; ">I have included three images: two of the sphere generated with gmsh, and another visualized within MATLAB.<br><br>I can't figure out what is causing the discrepancy. I have tried using both monomial and orthonormal (koorwinder) nodal basis functions, but that doesn't seem to help.<span class="Apple-converted-space"> </span><br><br>What type of basis functions does gmsh use? Is it producing its visualization using Bezier curves or nurbs? That's the only explanation I can think of.<br><br></div></div></blockquote><div><br></div><div><br></div>we use equidistant lagrange basis (see the documentation for node ordering)<br><div><br></div><div><br></div><div>JF</div><br><blockquote type="cite"><div ocsi="0" fpstyle="1" style="font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="direction: ltr; font-family: Tahoma; font-size: 10pt; ">I am using the following command to generate the mesh:<br>gmsh -2 -order 5 -hoOptimize sphere_generator.geo<br><br>Does anyone know what might be going on here? Any help would be greatly appreciated!<br><br>Thanks for your time!<br><br>John<br></div><span><mysphereplot.png></span><span><spherewedges.png></span><span><spherewoedges.png></span>_______________________________________________<br>gmsh mailing list<br><a href="mailto:gmsh@geuz.org">gmsh@geuz.org</a><br><a href="http://www.geuz.org/mailman/listinfo/gmsh">http://www.geuz.org/mailman/listinfo/gmsh</a><br></div></blockquote></div><br><div apple-content-edited="true">

<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div><span class="Apple-style-span" style="font-size: medium; "><div>------------------------------------------------------------------</div><div><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Prof. Jean-Francois Remacle</span></font></div><div><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Universite catholique de Louvain (UCL)</span></font></div><div><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Ecole Polytechnique de Louvain (EPL) - Louvain School of Engineering</span></font></div><div><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Institute of Mechanics, Materials and Civil Engineering (iMMC)</span></font></div><div><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Center for Systems Engineering and Applied Mechanics (CESAME)</span></font></div><div><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Tel : +32-10-472352 -- Mobile : +32-473-909930 </span></font></div><div><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "><br></span></font></div><div><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "> </span></font></div><div><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "><br></span></font></div><div><br></div></span></div></div></div></span></div><br class="Apple-interchange-newline"></span><br class="Apple-interchange-newline">

</div>

<br></body></html>