<div dir="ltr">Hi, I enjoy gmsh but I ran into a bug I think.<div><br></div><div>My test.stl file is closed and contains no self-intersections. In fact, gmsh will tetrahedralize it! But my FEA diverges, so I'd like to remesh the four surfaces in the "multisolid" STL to improve and simplify the 3D mesh. I just can't get the remeshing to work, except for STL "solid" number 3, which becomes surface(4) in gmsh when merged.<div><div><br></div><div><div>~/bin/gmsh-2.8.5-Linux/bin/gmsh -o test_remeshed.msh -2 remesh.geo </div></div><div><br></div><div>gives me </div><div><br></div><div>...</div><div>Error   : Unknown surface 1<br></div><div><div>Warning : Could not make generatrices list for compound surface 252</div><div>Error   : Unknown surface 2</div><div>Warning : Could not make generatrices list for compound surface 253</div><div>Error   : Unknown surface 3</div><div>Warning : Could not make generatrices list for compound surface 254</div></div><div>...</div><div><br></div><div>Note that surface 4 is *not* reported, in fact, if I do this surface alone, I can remesh it. The .geo file is shown below, it should obviously work (?). Also, full log is enclosed below. An additional complaint would be the infinite loop that gmsh dives into in this case (^C at the end)!<br></div></div><div><br></div><div>Did I do something overwhelmingly stupid? Or is this a bug?</div><div><br></div><div>/si</div><div><font face="monospace, monospace">------------------------------------------------------------------------------</font></div><div><font face="monospace, monospace"><br></font></div><div><div><font face="monospace, monospace">Mesh.CharacteristicLengthFactor=2;</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Mesh.RemeshParametrization=1;</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Mesh.RemeshAlgorithm=1; </font></div><div><font face="monospace, monospace">//--------------------------------------------</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">Merge "test.stl";</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">// Surface(i), for i {1 : number_surfaces}, results from merging stl!</font></div><div><font face="monospace, monospace">// How many are there?</font></div><div><font face="monospace, monospace">ss[] = Surface "*";</font></div><div><font face="monospace, monospace">number_surfaces = #ss[];</font></div><div><font face="monospace, monospace">Printf("#Surfaces: %g", number_surfaces);</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">CreateTopology;</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">new_cs = news;</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">For i In {1 : number_surfaces}</font></div><div><font face="monospace, monospace">  Compound Surface(new_cs + i) = {i};</font></div><div><font face="monospace, monospace">EndFor</font></div><div><font face="monospace, monospace"><br></font></div><div><font face="monospace, monospace">// Express each stl surface as a physical one</font></div><div><font face="monospace, monospace">For i In {1 : number_surfaces}</font></div><div><font face="monospace, monospace">  Physical Surface(i) = {new_cs + i};</font></div><div><font face="monospace, monospace">EndFor</font></div><div><br></div></div><div>-----------------------------------</div><div><br></div><div><br></div><div><br></div><div><div>Info    : Running '/home/simphead/bin/gmsh-2.8.5-Linux/bin/gmsh -o test_remeshed.msh -2 remesh.geo' [Gmsh 2.8.5, 1 node, max. 1 thread]</div><div>Info    : Started on Wed Jan  7 10:57:57 2015</div><div>Info    : Reading 'remesh.geo'...</div><div>Info    : Reading 'test.stl'...</div><div>Info    : 11244 facets in solid 0</div><div>Info    : 141845 facets in solid 1</div><div>Info    : 1013 facets in solid 2</div><div>Info    : 196 facets in solid 3</div><div>Info    : Done reading 'test.stl'</div><div>#Surfaces: 4</div><div>Info    : Creating topology from mesh...</div><div>Info    : Removing duplicate mesh vertices...</div><div>Info    : Found 0 duplicate vertices </div><div>Info    : No duplicate vertices found</div><div>Info    : Done creating topology from mesh (2.12813 s)</div><div>Error   : Unknown surface 1</div><div>Warning : Could not make generatrices list for compound surface 252</div><div>Error   : Unknown surface 2</div><div>Warning : Could not make generatrices list for compound surface 253</div><div>Error   : Unknown surface 3</div><div>Warning : Could not make generatrices list for compound surface 254</div><div>Info    : Done reading 'remesh.geo'</div><div>Info    : Meshing 1D...</div><div>Info    : Done meshing 1D (0 s)</div><div>Info    : Meshing 2D...</div><div>Info    : Meshing surface 252 (Compound surface, MeshAdapt)</div><div>Info    : Wrong topology: Genus=1, Nb boundaries=0, AR=-nan</div><div>Info    : -----------------------------------------------------------</div><div>Info    : --- Split surface 252 in 3 parts with Multilevel Mesh partitioner</div><div>Info    : Building graph...</div><div>Error   : No mesh elements were found</div><div>Info    : Partitioning graph...</div><div>Info    : Multiscale Partition SUCCESSFULLY PERFORMED : 0 parts (-1.86049e-16 s)</div><div>Info    : *** Starting parametrize compounds:</div><div>Info    : *** Parametrize compounds done (0.004 s)</div><div>Info    : *** Starting meshing 1D edges ...:</div><div>Info    : *** Meshing 1D edges done (-3.41307e-16s)</div><div>Info    : *** Starting Mesh of surface 252 ...</div><div>Info    : *** Mesh of surface 252 done by assembly 0 remeshed faces (-3.41307e-16 s)</div><div>Info    : -----------------------------------------------------------</div><div>Info    : Re-orient all 0 face normals coherently</div><div><br></div><div><br></div><div>^C</div></div><div><br></div></div></div>