Discovered bug in nodes numbering
Hey guys, how are you?
I discovered a minor bug in nodes numbering.
I created next mesh:
// Gmsh project created on Thu Jan 14 13:35:09 2021
boxSize = 5.0;
p = 1;
Point(p) = {-boxSize, -boxSize, 0.0, 1.0}; p += 1;
Point(p) = { boxSize, -boxSize, 0.0, 1.0}; p += 1;
Point(p) = { boxSize, boxSize, 0.0, 1.0}; p += 1;
Point(p) = {-boxSize, boxSize, 0.0, 1.0}; p += 1;
l = 1;
Line(l) = {1, 2}; l += 1;
Line(l) = {2, 3}; l += 1;
Line(l) = {3, 4}; l += 1;
Line(l) = {4, 1}; l += 1;
Curve Loop(1) = {1, 2, 3, 4};
Plane Surface(1) = {1};
Transfinite Curve {1, 2, 3, 4} = 3 Using Progression 1;
Periodic Line {3} = {1} Translate {0, 2*boxSize, 0};
Periodic Line {2} = {4} Translate {2*boxSize, 0, 0};
Physical Curve("Periodic1to3", 1) = {1, 3};
Physical Curve("Periodic2to4", 2) = {2, 4};
Physical Surface("Fluid", 3) = {1};
I created 2D mesh using front2d algorithm.
In .msh file, in elements, I got next two elements:
13 3 7 12
14 6 3 12
but when I presented cells numbers and nodes numbers, I got that the nodes swapped between this two elements, as you can see in the attached picture:
The same issue happens when I the "Transfinite Curve" setting to 5, this is what I've got:
42 3 11 29
43 10 3 29
and as before, you can see in attached picture, the nodes between these two cells are swapped:
Hope that you will find the problem.
Thank you in advance! :)
Edited by Christophe Geuzaine