<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="">Dear all, </div>
<div class=""><br class="">
</div>
<div class="">I am not familiar with FreeFEM++...</div>
<div class="">The mesh generation happens without any issue. </div>
<div class=""><br class="">
</div>
<div class="">If the problem is the normal orientation, maybe you could invert the normal at the end of the geo file, for whatever surface you want.</div>
<div class="">e.g. with the file you sent, for having outward normals in the mesh you should add</div>
<div class="">Reverse Surface{1,13,15,14}; </div>
<span class="">
<div class=""><span class=""><br class="">
</span></div>
<div class="">Alternative to changing the loops.</div>
<div class="">Best regards, </div>
<div class="">Ruth </div>
<br class="">
</span><span class=""><br class="">
</span>
<div class=""><br class="">
</div>
<br class="">
<div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
—<br class="">
Prof. Ruth V. Sabariego<br class="">
KU Leuven, Dept. Electrical Engineering ESAT/Electa, EnergyVille<br class="">
<a href="http://www.esat.kuleuven.be/electa" class="">http://www.esat.kuleuven.be/electa</a></div>
<div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<a href="http://www.energyville.be" class="">http://www.energyville.be</a></div>
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<br class="">
</div>
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
<br class="">
<br class="">
<br class="">
<br class="">
<br class="">
<br class="">
</div>
</div>
</div>
</div>
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 15 May 2020, at 08:55, Jose Juan Alonso del Rosario <<a href="mailto:josejuan.alonso@gm.uca.es" class="">josejuan.alonso@gm.uca.es</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" style="caret-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; text-decoration: none;" class="">
<div class="">Dear Ugis and colleages of GMSH,</div>
<div class=""><br class="">
</div>
<div class="">Thank you for your suggestions. I have changed the loop orientation of the upper and lower faces of the pentagon and the whole volumen was meshed (no hole at the end). In addition the same error still happens.</div>
<div class=""><br class="">
</div>
<div class="">"read mesh ok 0Mesh3, num Tetra:= 24011, num Vertice:= 5665 num boundary Triangles:= 6078<br class="">
Mesh3::meshS, num Triangles:= 6078, num Vertice:= 3039 num boundary Edges:= 354<br class="">
read mesh ok 0Mesh3, num Tetra:= 24011, num Vertice:= 5665 num boundary Triangles:= 6078<br class="">
Mesh3::meshS, num Triangles:= 6078, num Vertice:= 3039 num boundary Edges:= 354<br class="">
Bad orientation: The adj border element defined by [ 3 210 ] is oriented in the same direction in element 3519 and in the element 367 ****** bug in mesh construction? orientation parameter?<span class="Apple-converted-space"> </span><br class="">
S, num Triangles:= 6078, num Vertice:= 3039 num boundary Edges:= 354<br class="">
Bad orientation: The adj border element defined by [ 3 210 ] is oriented in the same direction in element 3519 and in the element 367 ****** bug in mesh construction? orientation parameter? "</div>
<div class=""><br class="">
</div>
<div class="">This is after several Optimize3D and Optimize 3d with NetGen and reading it with freefem++<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">I include the .geo code. It is a thick slice with a pentagonal hole. The gmsh is now 4.5.6 and the meshing options are MeshAdapt and Frontal. I have played with the signs of plane surfaces 14 and 15.</div>
<div class=""><br class="">
</div>
<div class="">Thanks in advance!!!</div>
<div class=""><br class="">
</div>
<div class="">Jose<br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
</div>
<br style="caret-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; text-decoration: none;" class="">
<div class="gmail_quote" style="caret-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; text-decoration: none;">
<div dir="ltr" class="gmail_attr">El jue., 14 may. 2020 a las 13:54, Ugis Lacis (<<a href="mailto:ugis.lacis@gmail.com" class="">ugis.lacis@gmail.com</a>>) escribió:<br class="">
</div>
<blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;">
<div class="">Dear Jose,<br class="">
<br class="">
I think I have had similar error, when importing GMSH mesh into FreeFEM++, few years ago (although it seems that the error message might have changed). It turned out that the direction of loop definition in GMSH was the issue for me. See below my exchange with
FreeFEM++ mailing list.<br class="">
<br class="">
Maybe it can be helpful for you.<br class="">
<br class="">
Best regards,<br class="">
Ugis Lacis<br class="">
Researcher<br class="">
KTH Mechanics<br class="">
<br class="">
------------------------------------------------------------------------------------------------------------------------------------------------------------------<br class="">
-------- Forwarded Message --------
<table cellspacing="0" cellpadding="0" border="0" class="">
<tbody class="">
<tr class="">
<th valign="BASELINE" nowrap="" align="RIGHT" class="">Subject:<span class="Apple-converted-space"> </span></th>
<td class="">Re: [Freefempp] Issue with loading mesh from GMSH.</td>
</tr>
<tr class="">
<th valign="BASELINE" nowrap="" align="RIGHT" class="">Date:<span class="Apple-converted-space"> </span></th>
<td class="">Fri, 3 Aug 2018 09:41:29 +0200</td>
</tr>
<tr class="">
<th valign="BASELINE" nowrap="" align="RIGHT" class="">From:<span class="Apple-converted-space"> </span></th>
<td class="">Uģis Lācis<span class="Apple-converted-space"> </span><a href="mailto:ugis.lacis@gmail.com" target="_blank" class=""><ugis.lacis@gmail.com></a></td>
</tr>
<tr class="">
<th valign="BASELINE" nowrap="" align="RIGHT" class="">To:<span class="Apple-converted-space"> </span></th>
<td class="">FreeFem mailing list<span class="Apple-converted-space"> </span><a href="mailto:freefempp@ljll.math.upmc.fr" target="_blank" class=""><freefempp@ljll.math.upmc.fr></a></td>
</tr>
</tbody>
</table>
<br class="">
<br class="">
Hi everyone,<br class="">
<br class="">
I managed to solve the issue myself. It turned out that the direction of GMSH loop was surviving till FreeFEM++ and while GMSH did not complain about anything, FreeFEM++ couldn't use the mesh. mesh2 could be fixed by changing signs in all lines within the line
loop, i.e.,<br class="">
change<br class="">
Line Loop(5) = {18, 19, 20, 12, 7, 8};<br class="">
to<br class="">
Line Loop(5) = {-18, -19, -20, -12, -7, -8};<br class="">
<br class="">
Maybe this helps to someone at some point.<br class="">
<br class="">
Best,<br class="">
Ugis<br class="">
<br class="">
On 02/08/18 15:30, Uģis Lācis wrote:<br class="">
Dear FreeFEM++ users and developers,<br class="">
<br class="">
I have run into an issue which I am unable to solve by myself. I have been using FreeFEM++ together with GMSH without any problems for quite some time. Now, however, trying to load a 2D mesh generated by GMSH, gives me following error:<br class="">
<br class="">
Assertion fail : (area>0)<br class="">
line :281, in file ./include/fem.hpp<br class="">
<br class="">
I attach the corresponding geo and msh files of two meshes (mesh1 - working fine, mesh2 - does not load), as well as edp file for reading and plotting the mesh.<br class="">
<br class="">
I have carefully checked that physical tags are applied correctly, I have tried exporting both in MSH and MESH (medit) formats, using GMSH versions 3.06 and 2.16.0, as well as using FreeFEM++ versions 3.560000 and 3.610002, all try has led me to the same behaviour.
I have identified that the issue comes from the surface mesh, if I export only line mesh (comment out "Physical Surface(11) = {3,4};" from mesh2.geo and run "gmsh -2 mesh2.geo"), then it loads fine. But checking the mesh within GMSH, everything seems to be
perfectly fine...<br class="">
<br class="">
Do you have any suggestions on how to solve this issue?<br class="">
<br class="">
Best regards,<br class="">
Ugis Lacis<br class="">
Researcher<br class="">
KTH Mechanics<br class="">
------------------------------------------------------------------------------------------------------------------------------------------------------------------<br class="">
<br class="">
<br class="">
<div class="">On 14/05/2020 12:03, Jose Juan Alonso del Rosario wrote:<br class="">
</div>
<blockquote type="cite" class="">
<div dir="ltr" class="">
<div class="">Dear colleages</div>
<div class=""><br class="">
</div>
<div class="">I am trying to build the mesh for a slice with a pentagonal hole.</div>
<div class="">GMSH runs properly with good quality tets but when trying to read it from freefem++ the error is:</div>
<div class=""><br class="">
</div>
<div class="">" Bad orientation: The adj border element defined by [ 16 246 ] is oriented in the same direction in element 3922 and in the element 446 ****** bug in mesh construction? orientation parameter?<span class="Apple-converted-space"> </span><br class="">
current line = 12 mpirank 1 / 2<br class="">
</div>
<div class="">"</div>
<div class="">Can you give any advice on this?</div>
<div class=""><br class="">
</div>
<div class="">Thanks in advance<br class="">
</div>
<div class=""><br class="">
</div>
<div class="">Jose<br class="">
</div>
<div class=""><br class="">
--<span class="Apple-converted-space"> </span><br class="">
<div dir="ltr" class="">
<div dir="ltr" class="">Prof. Dr. Jose Juan Alonso del Rosario<br class="">
Applied Physics Dept.,
<div class="">Naval and Oceanic Engineering. </div>
<div class="">University of Cadiz.<br class="">
Avda Rep Saharaui s/n. Puerto Real, 11510, Cadiz, Spain<br class="">
Tlf: +34 956 016054<br class="">
Fax: +34 956016079<br class="">
Secretary: +34 956016078<br class="">
<br class="">
</div>
<div class="">"All I wanna say is that they don't really care about us", Michael Jackson<br class="">
</div>
</div>
</div>
</div>
</div>
<br class="">
<fieldset class=""></fieldset>
<pre class="">_______________________________________________
gmsh mailing list
<a href="mailto:gmsh@onelab.info" target="_blank" class="">gmsh@onelab.info</a>
<a href="http://onelab.info/mailman/listinfo/gmsh" target="_blank" class="">http://onelab.info/mailman/listinfo/gmsh</a>
</pre>
</blockquote>
<br class="">
</div>
</blockquote>
</div>
<br clear="all" style="caret-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; text-decoration: none;" class="">
<br style="caret-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; text-decoration: none;" class="">
<span style="caret-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; text-decoration: none; float: none; display: inline !important;" class="">--<span class="Apple-converted-space"> </span></span><br style="caret-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; text-decoration: none;" class="">
<div dir="ltr" class="gmail_signature" style="caret-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; text-decoration: none;">
<div dir="ltr" class="">Prof. Dr. Jose Juan Alonso del Rosario<br class="">
Applied Physics Dept.,
<div class="">Naval and Oceanic Engineering. </div>
<div class="">University of Cadiz.<br class="">
Avda Rep Saharaui s/n. Puerto Real, 11510, Cadiz, Spain<br class="">
Tlf: +34 956 016054<br class="">
Fax: +34 956016079<br class="">
Secretary: +34 956016078<br class="">
<br class="">
</div>
<div class="">"All I wanna say is that they don't really care about us", Michael Jackson<br class="">
</div>
</div>
</div>
<span id="cid:f_ka7ua43g0"><pentagono.geo></span><span style="caret-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; text-decoration: none; float: none; display: inline !important;" class="">_______________________________________________</span><br style="caret-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; text-decoration: none;" class="">
<span style="caret-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; text-decoration: none; float: none; display: inline !important;" class="">gmsh
mailing list</span><br style="caret-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; text-decoration: none;" class="">
<a href="mailto:gmsh@onelab.info" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">gmsh@onelab.info</a><br style="caret-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; text-decoration: none;" class="">
<a href="http://onelab.info/mailman/listinfo/gmsh" style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">http://onelab.info/mailman/listinfo/gmsh</a></div>
</blockquote>
</div>
<br class="">
</body>
</html>