[Gmsh] high order p=10 tetrahedron : a internal node numbering problem ?

Pierre Saramito Pierre.Saramito at imag.fr
Tue Oct 18 16:18:15 CEST 2011


Hi Everybody,

I perhaps founded a bug while using high order tetrahedra meshes.
I just generated by hand the small "tetra_p1.msh" file (attached).

Then :
  gmsh -3 -order 10 tetra_p1.msh -o tetra_p10.msh
       => the file is well generated
	  but it could be a problem with the most nested
	  internal node numbers: something like a bad
	  node permutation.
  gmsh tetra_p10.msh
       => core dump.
  gmsh -check tetra_p10.msh
	=> ok
  gmsh -version
	2.5.0 

Note that for p=9 all is ok:
  gmsh -3 -order 9 tetra_p1.msh -o tetra_p9.msh
  gmsh tetra_p9.msh

Please could you take an eye on the node numbering when p=10 ?

I take this opportunity to congratulate all the developers
for the wonderfull features in gmsh.

Regards,

Pierre
-- 
Pierre.Saramito at imag.fr
Directeur de Recherche CNRS
Laboratoire Jean Kuntzmann, Grenoble, France
http://www-ljk.imag.fr/membres/Pierre.Saramito
-------------- next part --------------
A non-text attachment was scrubbed...
Name: tetra_p1.msh
Type: model/mesh
Size: 131 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20111018/903cc339/attachment.msh>