[Gmsh] A bug is found in meshing information!

Christophe Geuzaine cgeuzaine at ulg.ac.be
Wed Jan 30 08:43:44 CET 2008


ZHUANG X. wrote:
> 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

These stats are for the initial mesh only. Then the mesh gets refined; 
the final number of tets is given in the Tools->Statistics window.



> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh


-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine