Problem with the generation of high-order Pyramids
Dear all,
In order to check the geometric configuration of high-order elements, we tried to generate higher order elements starting from a single first order element. It worked fine for Tetrahedra, Hexaedra and Prisms but we think we might have detected some errors in the generation of high-order Pyramids.
We started from the order 1 Pyramid meshfile:
(Pyramid1.msh)
$MeshFormat 2.2 0 8 $EndMeshFormat $Nodes 5 1 -1 -1 0 2 1 -1 0 3 1 1 0 4 -1 1 0 5 0 0 1
$EndNodes $Elements 1 1 7 2 10 8 1 2 3 4 5 $EndElements
When we try to build the higher orders (using the command line "gmsh -order Pyramid1.msh"), some points are not correctly located: -The center of the quadrangular face at order 2 (Point 14) -The interior points starting from order 3 (Points 26 to 30)
I have attached the .msh files generated by gmsh for orders 2 and 3
Also, when trying to start from a valid order 2-Pyramid to generate an order 3 one, a segmentation fault is returned.
Best regards,
Hamza