[Gmsh] A bug is found in meshing information!

ZHUANG X. xiao-ying.zhuang at durham.ac.uk
Thu Nov 29 17:08:27 CET 2007


GMSH working team
 
Hello!
 
I have  found a bug in meshing information.
I found the tetrahedra element counter was wrong.
There are actully 102 tetrahedra elements generted while the statistics given in meshing process gave us 101 (see below).
I call gmsh program in Matlab and use the element and node information for analysis. 
I found that the meshing information given by gmsh in counting. See the attached file test1.msh. There are 102 tetrahedra elements generated subtantiilly from 117 to 218.
218-117+1 = 102. But the statistics informatio gave us 101.
The .gmsh and .geo file is attached in this mail.

EFG3DAnalysis(gmsh_file_name, test_model)  //that is the matlab command line
eval(strcat('!','F:\gmsh-2.0.8-Windows\gmsh',gmsh_file_name,' -3')); //that is the matlab command line

Info    : 'F:\gmsh-2.0.8-Windows\gmsh F:\GMSH_TEST\test1.geo -3 ' started on Thu Nov 29 15:56:34 2007 
Info    : Reading 'F:\GMSH_TEST\test1.geo' 
Info    : Read 'F:\GMSH_TEST\test1.geo' 
Info    : Meshing 1D... 
Info    : Meshing curve 1 
Info    : Meshing curve 2 
Info    : Meshing curve 3 
Info    : Meshing curve 4 
Info    : Meshing curve 5 
Info    : Meshing curve 6 
Info    : Meshing curve 7 
Info    : Meshing curve 8 
Info    : Meshing curve 9 
Info    : Meshing curve 10 
Info    : Meshing curve 11 
Info    : Meshing curve 12 
Info    : Mesh 1D complete (0 s) 
Info    : Mesh 
Info    : Meshing 2D... 
Info    : Meshing surface 14 (Plane, MeshAdapt+Delaunay) 
Info    : Meshing surface 16 (Plane, MeshAdapt+Delaunay) 
Info    : Meshing surface 18 (Plane, MeshAdapt+Delaunay) 
Info    : Meshing surface 20 (Plane, MeshAdapt+Delaunay) 
Info    : Meshing surface 22 (Plane, MeshAdapt+Delaunay) 
Info    : Meshing surface 24 (Plane, MeshAdapt+Delaunay) 
Info    : Mesh 2D complete (0 s) 
Info    : Mesh 
Info    : Meshing 3D... 
Info    : Meshing volume 26 (Delaunay) 
Info    : 44 points 24 edges and 84 faces in the final mesh 
Info    : 44 points created -- Worst tet radius is 1.52745 
Info    : Mesh 3D complete (0.015625 s) 
Info    : Mesh 
Info    : 45 vertices 210 elements 
Info    : Writing 'F:\GMSH_TEST\test1.msh' 
Info    : Wrote 'F:\GMSH_TEST\test1.msh' 
Constructing Delaunay tetrahedralization. 
Delaunay seconds:  0.015 
Creating surface mesh. 
Perturbing vertices. 
Delaunizing segments. 
Constraining facets. 
Segment and facet seconds:  0 
Removing unwanted tetrahedra. 
Hole seconds:  0 
Repairing mesh. 
Repair seconds:  0 
 
Writing nodes. 
Writing elements. 
Writing faces. 
Writing edges. 
 
Output seconds:  0 
Total running seconds:  0.015 
 
Statistics: 
 
  Input points: 44 
  Input facets: 84 
  Input segments: 126 
  Input holes: 0 
  Input regions: 0 
 
  Mesh points: 44 
  Mesh tetrahedra: 101 
  Mesh triangles: 244 
  Mesh subfaces: 84 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20071129/09db37ea/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test1.geo
Type: application/octet-stream
Size: 850 bytes
Desc: test1.geo
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20071129/09db37ea/attachment.geo>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: test1.mesh
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20071129/09db37ea/attachment.mesh>