
<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">

<meta name="Generator" content="Microsoft Word 14 (filtered medium)">

<style><!--

/* Font Definitions */

@font-face

        {font-family:Calibri;

        panose-1:2 15 5 2 2 2 4 3 2 4;}

/* Style Definitions */

p.MsoNormal, li.MsoNormal, div.MsoNormal

        {margin:0in;

        margin-bottom:.0001pt;

        font-size:11.0pt;

        font-family:"Calibri","sans-serif";}

a:link, span.MsoHyperlink

        {mso-style-priority:99;

        color:blue;

        text-decoration:underline;}

a:visited, span.MsoHyperlinkFollowed

        {mso-style-priority:99;

        color:purple;

        text-decoration:underline;}

span.EmailStyle17

        {mso-style-type:personal-compose;

        font-family:"Calibri","sans-serif";

        color:windowtext;}

.MsoChpDefault

        {mso-style-type:export-only;

        font-family:"Calibri","sans-serif";}

@page WordSection1

        {size:8.5in 11.0in;

        margin:1.0in 1.0in 1.0in 1.0in;}

div.WordSection1

        {page:WordSection1;}

--></style><!--[if gte mso 9]><xml>

<o:shapedefaults v:ext="edit" spidmax="1026" />

</xml><![endif]--><!--[if gte mso 9]><xml>

<o:shapelayout v:ext="edit">

<o:idmap v:ext="edit" data="1" />

</o:shapelayout></xml><![endif]-->

</head>

<body lang="EN-US" link="blue" vlink="purple">

<div class="WordSection1">

<p class="MsoNormal">Hello,<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">I use gmsh version 2.9.0 and came across some unexpected behavior:<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">Attached are two example .geo files for a simple geometry with two parts: a square on the left and a square on the right hand side, sharing one edge. I want the boundary of both squares to be oriented counter clockwise. Thus, the shared

 edge has to be reversed for one square (the one on the right hand side in the example). Based on the manual, this is done adding a minus sign in the  “Physical Line” command. This works fine for linear and quadratic meshes. But when I look at 3<sup>rd</sup>

 or 4<sup>th</sup> order meshes, I have the problem, that the interior edge nodes are not reversed. For the example.msh and the example_reversed.msh I get the following node ordering of line element 3:<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">3<sup>rd</sup> order: 2,3,9,10 gets reversed to 3,2,9,10<o:p></o:p></p>

<p class="MsoNormal">4<sup>th</sup> order: 2,3,10,11,12 gets reversed to 3,2,10,11,12<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">Instead, I expected the inner edge nodes to be reversed as well, such that

<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">3<sup>rd</sup> order: 2,3,9,10  would get reversed to 3,2,10,9<o:p></o:p></p>

<p class="MsoNormal">4<sup>th</sup> order: 2,3,10,11,12 would get reversed to 3,2,12,11,10                  

<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">Is this a bug? If not, is there a way to obtain the behavior I expected?<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

<p class="MsoNormal">Thank you for your help!<o:p></o:p></p>

<p class="MsoNormal">Kind regards,<o:p></o:p></p>

<p class="MsoNormal">Verena<o:p></o:p></p>

<p class="MsoNormal"><o:p> </o:p></p>

</div>

</body>

</html>