[Gmsh] Surface and volume mesh do not match on un-recombined transfinite mesh

[Adam Croft]

Hello,

As noted in the subject line, I am having a problem with the surface 
mesh not matching the volume mesh in un-recombined transfinite meshes. I 
have included a small example where you can see that the +x and +y faces 
are split one way by the triangles and the other way by the tetrahedra, 
the other faces are correct. I have tried orienting the faces a couple 
of different ways, but haven't had any luck. Am I doing something wrong, 
or is this a bug?

Thanks,

Adam Croft

-------------- next part --------------
sp=1.0;
xmin=-3.0;
xmax=3.0;
ymin=-2.0;
ymax=4.0;
zmin=-1.5;
zmax=4.5;
nx=2;
ny=2;
nz=2;
Point(1) = {xmin,ymin,zmin,sp};
Point(2) = {xmax,ymin,zmin,sp};
Point(3) = {xmin,ymax,zmin,sp};
Point(4) = {xmax,ymax,zmin,sp};
Point(5) = {xmin,ymin,zmax,sp};
Point(6) = {xmax,ymin,zmax,sp};
Point(7) = {xmin,ymax,zmax,sp};
Point(8) = {xmax,ymax,zmax,sp};
Line(1) = {1, 3};
Transfinite Line(1) = ny;
Line(2) = {3, 4};
Transfinite Line(2) = nx;
Line(3) = {4, 2};
Transfinite Line(3) = ny;
Line(4) = {2, 1};
Transfinite Line(4) = nx;
Line(5) = {1, 5};
Transfinite Line(5) = nz;
Line(6) = {2, 6};
Transfinite Line(6) = nz;
Line(7) = {4, 8};
Transfinite Line(7) = nz;
Line(8) = {3, 7};
Transfinite Line(8) = nz;
Line(9) = {5, 6};
Transfinite Line(9) = nx;
Line(10) = {6, 8};
Transfinite Line(10) = ny;
Line(11) = {8, 7};
Transfinite Line(11) = nx;
Line(12) = {7, 5};
Transfinite Line(12) = ny;
Line Loop(13) = {9, 10, 11, 12};
Plane Surface(14) = {13};
Transfinite Surface(14) = {5,6,8,7}; //front
Line Loop(15) = {4, 1, 2, 3};
Plane Surface(16) = {15};
Transfinite Surface(16) = {1,3,4,2}; //back
Line Loop(17) = {-6,-3,7,-10};
Plane Surface(18) = {17};
Transfinite Surface(18) = {6,2,4,8}; //right
Line Loop(19) = {-2,8,-11,-7};
Plane Surface(20) = {19};
Transfinite Surface(20) = {4,3,7,8}; //top
Line Loop(21) = {-1,5,-12,-8};
Plane Surface(22) = {21};
Transfinite Surface(22) = {1,5,7,3}; //left
Line Loop(23) = {-4,6,-9,-5};
Plane Surface(24) = {23};
Transfinite Surface(24) = {1,2,6,5}; //bottom
Surface Loop(25) = {-24, -16, -22, 20, 18, 14};
Volume(26) = {25};
Transfinite Volume(26)={1,2,4,3,5,6,8,7};
Physical Surface("bottom") = {24};
Physical Surface("top") = {20};
Physical Surface("left") = {22};
Physical Surface("right") = {18};
Physical Surface("front") = {14};
Physical Surface("back") = {16};
Physical Volume("interior") = {26};
Coherence;
Coherence;
-------------- next part --------------
A non-text attachment was scrubbed...
Name: coarseblock.png
Type: image/png
Size: 27095 bytes
Desc: not available
URL: <http://onelab.info/pipermail/gmsh/attachments/20170907/ef76776d/attachment-0001.png>