Affine transformation requires at least 12 entires
Hello,
I have been trying to generate a triply periodic mesh on a unit cube following benchmarks/3d/CubePeriodic.geo
but have been getting the following error:
Error: 'outCoarse.geo', line 71 Affine transformation requires at least 12 entires
Here is the complete geo
file
// This code was created by pygmsh v6.0.4.
SetFactory("OpenCASCADE");
p0 = newp;
Point(p0) = {0.0, 0.0, 0.0, 1};
p1 = newp;
Point(p1) = {1.0, 0.0, 0.0, 1};
p2 = newp;
Point(p2) = {0.0, 1.0, 0.0, 1};
p3 = newp;
Point(p3) = {0.0, 0.0, 1.0, 1};
p4 = newp;
Point(p4) = {1.0, 0, 1, 1};
p5 = newp;
Point(p5) = {1.0, 1.0, 0, 1};
p6 = newp;
Point(p6) = {1, 1, 1, 1};
p7 = newp;
Point(p7) = {0, 1, 1, 1};
l0 = newl;
Line(l0) = {p0, p1};
l1 = newl;
Line(l1) = {p0, p2};
l2 = newl;
Line(l2) = {p0, p3};
l3 = newl;
Line(l3) = {p1, p5};
l4 = newl;
Line(l4) = {p5, p6};
l5 = newl;
Line(l5) = {p6, p4};
l6 = newl;
Line(l6) = {p4, p1};
l7 = newl;
Line(l7) = {p2, p7};
l8 = newl;
Line(l8) = {p7, p6};
l9 = newl;
Line(l9) = {p5, p2};
l10 = newl;
Line(l10) = {p3, p4};
l11 = newl;
Line(l11) = {p7, p3};
ll0 = newll;
Line Loop(ll0) = {l3, l4, l5, l6};
s0 = news;
Plane Surface(s0) = {ll0};
ll1 = newll;
Line Loop(ll1) = {l1, l7, l11, l2};
s1 = news;
Plane Surface(s1) = {ll1};
ll2 = newll;
Line Loop(ll2) = {l7, l8, l4, l9};
s2 = news;
Plane Surface(s2) = {ll2};
ll3 = newll;
Line Loop(ll3) = {l2, l10, l6, l0};
s3 = news;
Plane Surface(s3) = {ll3};
ll4 = newll;
Line Loop(ll4) = {l0, l3, l9, l1};
s4 = news;
Plane Surface(s4) = {ll4};
ll5 = newll;
Line Loop(ll5) = {l10, l5, l8, l11};
s5 = news;
Plane Surface(s5) = {ll5};
sl0 = news;
Surface Loop(sl0) = {s0,s4,s1,s5,s2,s3};
vol0 = newv;
Volume(vol0) = {sl0};
Periodic Surface { s0 } = { s1 };
Periodic Surface { s2 } = { s3 };
Periodic Surface { s5 } = { s4 };
However if I proceed to generate the mesh ignoring it, I do seem to get a reasonable mesh.
Any pointers to what could be going wrong (and if at all the error could lead to a non-matching mesh)