[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