<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=utf-8">
<meta name="Generator" content="Microsoft Word 12 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:"Segoe UI";
        panose-1:2 11 5 2 4 2 4 2 2 3;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","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.StileMessaggioDiPostaElettronica17
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 2.0cm 2.0cm 2.0cm;}
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="IT" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hullo Ruth,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">but the same kind of operation works fine in other cases that I have. The surface loops that define the volume are not changed by the rotation!<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">The reason why I need to do that operation AFTER having defined the volume is that I have an automatic procedure which reads a geo file, perform
 a translation/rotation on one geometry component specified through its physical surfaces IDs, generates the mesh and carry out the simulation . This is done in an outer loop so that many geometry modifications are analysed in the same run. For this reason
 the procedure appends the commands at the end of the geo file.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">I guess the problem is related with the use of OpenCascade, because the same approach works fine in cases when the Built-In geometry kernel is
 used.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">A.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif"">Da:</span></b><span style="font-size:10.0pt;font-family:"Segoe UI","sans-serif""> Ruth Vazquez Sabariego [mailto:ruth.sabariego@kuleuven.be]
<br>
<b>Inviato:</b> venerdì 20 ottobre 2017 12:02<br>
<b>A:</b> Alessandro Vicini<br>
<b>Cc:</b> gmsh@onelab.info<br>
<b>Oggetto:</b> Re: [Gmsh] Problem in transformation of surfaces list<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Hi Alessandro, <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">The problem is the order of the operations. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">You do a rotation of surfaces that define volume 1.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">That’s why they are kept. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">The rotation is properly handled if you write it as follows:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">tmp[] = Physical Surface {102,103};<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Rotate {{1, 0, 0}, {-0.5, 0.65, -0.25}, Pi/6} {<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">   Surface {tmp[]};<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">}<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Surface Loop(5) = {28, 23, 25, 27, 26, 24};<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Surface Loop(6) = {9, 4, 6, 8, 7, 10, 1, 22, 18, 21, 20, 17, 19};<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Surface Loop(7) = {15, 11, 13, 12, 16, 14};<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Volume(1) = {5, 6, 7};<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Physical Volume(999) = {1}; <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Best regards, <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Ruth<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<div>
<p class="MsoNormal"><span style="color:black">—<br>
Prof. Ruth V. Sabariego<br>
KU Leuven  <br>
Dept. Electrical Engineering ESAT/Electa, EnergyVille<br>
<a href="https://urlsand.esvalabs.com/?u=http%3A%2F%2Fwww.esat.kuleuven.be%2Felecta&e=83f3a71b&h=2cc88284&f=y&p=y">http://www.esat.kuleuven.be/electa</a><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><a href="https://urlsand.esvalabs.com/?u=http%3A%2F%2Fwww.energyville.be&e=83f3a71b&h=5a4cb613&f=y&p=y">http://www.energyville.be</a><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="color:black">Free software: <a href="https://urlsand.esvalabs.com/?u=http%3A%2F%2Fgmsh.info&e=83f3a71b&h=b3002766&f=y&p=y">http://gmsh.info</a> | <a href="https://urlsand.esvalabs.com/?u=http%3A%2F%2Fgetdp.info&e=83f3a71b&h=993699b8&f=y&p=y">http://getdp.info</a> | <a href="https://urlsand.esvalabs.com/?u=http%3A%2F%2Fonelab.info&e=83f3a71b&h=017e41bc&f=y&p=y">http://onelab.info</a><br>
<br>
<br>
<br>
<br>
<br>
<o:p></o:p></span></p>
</div>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal">On 20 Oct 2017, at 10:46, Alessandro Vicini <<a href="mailto:alessandro.vicini@sitael.com">alessandro.vicini@sitael.com</a>> wrote:<o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal"><br>
Hullo everybody. In the last two lines of the attached geo file I retrieve the list of elementary surfaces ID from their physical IDs, and then apply a rotation to the list. However, only the edges are rotated, and the "old" surfaces and edges are not deleted
 after the operation. Why...? Thank you.<br>
<br>
Alessandro<br>
<br>
<br>
<Test.geo>_______________________________________________<br>
gmsh mailing list<br>
<a href="mailto:gmsh@onelab.info">gmsh@onelab.info</a><br>
<a href="http://onelab.info/mailman/listinfo/gmsh">http://onelab.info/mailman/listinfo/gmsh</a><o:p></o:p></p>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>