[Gmsh] GMSH fails to mesh cube with second order triangles

Michael Carley m.j.carley at bath.ac.uk
Fri Dec 11 12:34:03 CET 2009


Hello all,
           I am trying to mesh a cube surface with second order triangles. 
The first order mesh is fine, but the second order hangs. I am using 
the downloaded Linux binary of gmsh 2.4.2 and my input geometry is:

lc = 1.0 ;

Point(1) = {-1,-1,-1,lc} ;
Point(2) = { 1,-1,-1,lc} ;
Point(3) = { 1, 1,-1,lc} ;
Point(4) = {-1, 1,-1,lc} ;
Point(5) = {-1,-1, 1,lc} ;
Point(6) = { 1,-1, 1,lc} ;
Point(7) = { 1, 1, 1,lc} ;
Point(8) = {-1, 1, 1,lc} ;
Line(1) = {4, 8};
Line(2) = {8, 7};
Line(3) = {7, 3};
Line(4) = {3, 4};
Line(5) = {1, 1};
Line(6) = {1, 1};
Line(7) = {1, 5};
Line(8) = {5, 6};
Line(9) = {6, 2};
Line(10) = {2, 1};
Line(11) = {1, 4};
Line(12) = {8, 5};
Line(13) = {6, 7};
Line(14) = {3, 2};
Line Loop(15) = {4, 1, 2, 3};
Plane Surface(16) = {15};
Line Loop(17) = {4, -11, -10, -14};
Plane Surface(18) = {-17};
Line Loop(19) = {14, -9, 13, 3};
Plane Surface(20) = {-19};
Line Loop(21) = {2, -13, -8, -12};
Plane Surface(22) = {-21};
Line Loop(23) = {1, 12, -7, 11};
Plane Surface(24) = {-23};
Line Loop(25) = {8, 9, 10, 7};
Plane Surface(26) = {-25};


-- 
`To tell the truth, let us be honest at least, it is some considerable
 	  time since I last knew what I was talking about.'
No MS attachments: http://www.gnu.org/philosophy/no-word-attachments.html
Home page: http://people.bath.ac.uk/ensmjc/