<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:"MS Mincho";
panose-1:2 2 6 9 4 2 5 8 3 4;}
@font-face
{font-family:"MS Mincho";
panose-1:2 2 6 9 4 2 5 8 3 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:"\@MS Mincho";
panose-1:2 2 6 9 4 2 5 8 3 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
@font-face
{font-family:Verdana;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* 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;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:Consolas;}
span.EmailStyle19
{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:72.0pt 72.0pt 72.0pt 72.0pt;}
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-GB link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Hello Kate,<o:p></o:p></span></p><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 style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>The issue might be related to some recent bug in 2.9.3.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>I am saying this because if I try to remesh with that I get a crash, while remeshing works quite well with an older Gmsh that I have (2.8.6).<o:p></o:p></span></p><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 style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Also, I think you get a better mesh if you follow Tutorial 13 (<a href="http://geuz.org/gmsh/doc/texinfo/gmsh.html#t13_002egeo">http://geuz.org/gmsh/doc/texinfo/gmsh.html#t13_002egeo</a>) and reconstruct the edges as well.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>If not, only the interiors of the faces will remeshed.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>I’ve added an updated test.geo with does this.<o:p></o:p></span></p><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 style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Best regards,<br>Serban<o:p></o:p></span></p><p class=MsoNormal><span 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 lang=EN-US style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>From:</span></b><span lang=EN-US style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'> gmsh [mailto:gmsh-bounces@ace20.montefiore.ulg.ac.be] <b>On Behalf Of </b>look<br><b>Sent:</b> 25 June 2015 11:57<br><b>To:</b> gmsh@geuz.org<br><b>Subject:</b> [Gmsh] Remesh STL<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Verdana","sans-serif";background:#F7F7F7'>Hello everybody,</span><span style='font-size:10.0pt;font-family:"Verdana","sans-serif"'><br><br><span style='background:#F7F7F7'>I am trying to remesh a very easy trivial *.stl geometry which I previously generated in Blender. Therefore I followed the instructions from a gmsh mailing list post: </span></span><a href="http://geuz.org/pipermail/gmsh/2011/006887.html" target="_blank"><span style='font-size:10.0pt;font-family:"Verdana","sans-serif";color:#990066;text-decoration:none'>http://geuz.org/pipermail/gmsh/2011/006887.html</span></a><span style='font-size:10.0pt;font-family:"Verdana","sans-serif"'><br><br><span style='background:#F7F7F7'>1) Open inout.stl in gmsh (with gui)</span><br><span style='background:#F7F7F7'>2) Click on Mesh->Reclassify 2D</span><br><span style='background:#F7F7F7'>3) Click on "All" in step 1 and step 3, then click on "Reclassify"</span><br><span style='background:#F7F7F7'>4) Save the reclassified mesh as a new msh file -> reclassified_inout.msh</span><br><span style='background:#F7F7F7'>5) I can see several distinct surfaces when I open this new file and go to Tools/Visibility.</span><br><span style='background:#F7F7F7'>6) Write a new file test.geo:</span></span><o:p></o:p></p><div><div style='mso-element:para-border-div;border:inset 1.0pt;padding:5.0pt 5.0pt 5.0pt 5.0pt;background:#F7F7FF'><pre style='background:#F7F7FF;border:none;padding:0cm;overflow:auto'>Mesh.RemeshAlgorithm=1;<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Mesh.CharacteristicLengthFactor=0.05;<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'><o:p> </o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Merge "reclassified_inout.msh";<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'><o:p> </o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>CreateTopology;<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'><o:p> </o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(23)={3};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(24)={4};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(25)={5};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(26)={6};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(27)={7};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(28)={8};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(29)={9};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(30)={10};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(31)={11};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(32)={12};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(33)={13};<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Compound Surface(34)={14};<o:p></o:p></pre></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Verdana","sans-serif";background:#F7F7F7'>7) Run </span><o:p></o:p></p><div style='margin-left:15.0pt;margin-top:3.75pt;margin-right:15.0pt;margin-bottom:15.0pt'><div style='margin-bottom:1.5pt'><p class=MsoNormal><span style='font-size:8.5pt;font-family:"Verdana","sans-serif"'>Code:<o:p></o:p></span></p></div><div style='mso-element:para-border-div;border:inset 1.0pt;padding:5.0pt 5.0pt 5.0pt 5.0pt;background:#F7F7FF'><pre style='background:#F7F7FF;border:none;padding:0cm;overflow:auto'>gmsh -2 test.geo<o:p></o:p></pre></div></div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Verdana","sans-serif";background:#F7F7F7'>Unfortunately I get the following error message:</span><o:p></o:p></p><div style='margin-left:15.0pt;margin-top:3.75pt;margin-right:15.0pt;margin-bottom:15.0pt'><div style='margin-bottom:1.5pt'><p class=MsoNormal><span style='font-size:8.5pt;font-family:"Verdana","sans-serif"'>Code:<o:p></o:p></span></p></div><div style='mso-element:para-border-div;border:inset 1.0pt;padding:5.0pt 5.0pt 5.0pt 5.0pt;background:#F7F7FF'><pre style='background:#F7F7FF;border:none;padding:0cm;overflow:auto'>kate@kate-VirtualBox:/media/sf_documents/gmsh/Daten/inout_cubes_vonBlender$ gmsh -2 test.geo<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Running 'gmsh -2 test.geo' [Gmsh 2.9.3, 1 node, max. 1 thread]<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Started on Thu Jun 25 11:34:00 2015<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Reading 'test.geo'...<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Reading 'reclassified_inout.msh'...<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : 16 vertices<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Vertex numbering is dense<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : 36 elements<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Done reading 'reclassified_inout.msh'<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Creating topology from mesh...<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Removing duplicate mesh vertices...<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Found 0 duplicate vertices <o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : No duplicate vertices found<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Done creating topology from mesh (0.000554 s)<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Done reading 'test.geo'<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Meshing 1D...<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Done meshing 1D (1.7e-05 s)<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Meshing 2D...<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Meshing surface 23 (Compound surface, MeshAdapt)<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Info : Parametrizing surface 23 with 'convex map'<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: Object is in wrong state<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: Not for unassembled matrix<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: Petsc Release Version 3.5.2, Sep, 08, 2014 <o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: gmsh on a linux_complex_mumps_seq named kate-VirtualBox by kate Thu Jun 25 11:34:00 2015<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: Configure options --with-pic --with-blas-lapack-lib="[/usr/local/lib/liblapack.a,/usr/local/lib/libblas.a,-lpthread]" --with-clanguage=cxx --with-debugging=0 --with-mpi=0 --with-x=0 --with-ssl=0 --with-mpiuni-fortran-binding=0 --with-mumps-include=/home/geuzaine/src/mumps_seq/include --with-mumps-lib="[/home/geuzaine/src/mumps_seq/lib/libcmumps.a,/home/geuzaine/src/mumps_seq/lib/libsmumps.a,/home/geuzaine/src/mumps_seq/lib/libdmumps.a,/home/geuzaine/src/mumps_seq/lib/libzmumps.a,/home/geuzaine/src/mumps_seq/lib/libmumps_common.a,/home/geuzaine/src/mumps_seq/lib/libmetis.a,/home/geuzaine/src/mumps_seq/lib/libmpiseq.a,/home/geuzaine/src/mumps_seq/lib/libpord.a]" --with-mumps-serial --with-shared-libraries=0 --with-scalar-type=complex<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: #1 MatGetOrdering() line 189 in /home/geuzaine/src/petsc-3.5.2/src/mat/order/sorder.c<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: #2 PCSetUp_ILU() line 196 in /home/geuzaine/src/petsc-3.5.2/src/ksp/pc/impls/factor/ilu/ilu.c<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: #3 PCSetUp() line 902 in /home/geuzaine/src/petsc-3.5.2/src/ksp/pc/interface/precon.c<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: #4 KSPSetUp() line 305 in /home/geuzaine/src/petsc-3.5.2/src/ksp/ksp/interface/itfunc.c<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: #5 KSPSolve() line 417 in /home/geuzaine/src/petsc-3.5.2/src/ksp/ksp/interface/itfunc.c<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>[0]PETSC ERROR: #6 _try() line 9 in /home/geuzaine/src/gmsh/Solver/linearSystemPETSc.hpp<o:p></o:p></pre><pre style='background:#F7F7FF;border:none;padding:0cm'>Aborted (core dumped)<o:p></o:p></pre></div></div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Verdana","sans-serif";background:#F7F7F7'>I have attached the necessary files to reproduce the problem. I tried with gmsh 2.9.3 on OS X and gmsh 2.9.3 on Ubuntu.</span><span style='font-size:10.0pt;font-family:"Verdana","sans-serif"'><br><br><span style='background:#F7F7F7'>I would appreciate any help/hint/comment!</span><br><br><span style='background:#F7F7F7'>Best regards,</span><br><br><span style='background:#F7F7F7'>Kate</span></span><o:p></o:p></p></div></div></div></body></html>
<table><tr><td bgcolor=#ffffff><font color=#000000><pre>______________________________________________________________________
Fujitsu Laboratories of Europe Limited
Hayes Park Central, Hayes End Road, Hayes, Middlesex, UB4 8FE
Registered No. 4153469
This e-mail and any attachments are for the sole use of addressee(s) and
may contain information which is privileged and confidential. Unauthorised
use or copying for disclosure is strictly prohibited. The fact that this
e-mail has been scanned by Trendmicro Interscan does not guarantee that
it has not been intercepted or amended nor that it is virus-free.</pre></font></td></tr></table>