Package: gmsh
Version: 2.8.4+dfsg-1+b1
Severity: normal
Merging a .geo file and then translating the volume in it results in a
wrong mesh.
I've attached an example: the unit cube (0,0,0) to (1,1,1) is
translated along the z-axis. The translated cube is defined by the
corners (0,0,1) and (1,1,2). However the resulting mesh has a node on
the z=0 plane (where z should be 1) and another on the z=1 plane that
belongs on the z=2 plane.
Run "gmsh -3 cube1.geo" or choose Mesh->3D in the GUI to reproduce the
problem. One has to rotate the cube to see the problem. Or look at
these nodes in the .msh file:
9 0.5 0.5 0
10 0.5 0.5 1
The problem does not happen when using Include instead of Merge or
when including the contents of unit-cube.geo directly.
Ansgar
-- System Information:
Debian Release: jessie/sid
APT prefers testing
APT policy: (500, 'testing'), (100, 'unstable'), (1, 'experimental')
Architecture: amd64 (x86_64)
Kernel: Linux 3.14-1-amd64 (SMP w/8 CPU cores)
Locale: LANG=C.UTF-8, LC_CTYPE=C.UTF-8 (charmap=UTF-8)
Shell: /bin/sh linked to /bin/dash
Versions of packages gmsh depends on:
ii libann0 1.1.2+doc-5
ii libblas3 [libblas.so.3] 1.2.20110419-7
ii libc6 2.19-1
ii libfltk-gl1.3 1.3.2-5
ii libfltk-images1.3 1.3.2-5
ii libfltk1.3 1.3.2-5
ii libgcc1 1:4.9.0-5
ii libgl1-mesa-glx [libgl1] 10.1.4-1
ii libgl2ps0 1.3.8-1
ii libglu1-mesa [libglu1] 9.0.0-2
ii libjpeg8 8d-2
ii liblapack3 [liblapack.so.3] 3.5.0-2
ii libmed1 3.0.6-5
ii liboce-foundation8 0.15-4
ii liboce-modeling8 0.15-4
ii libopenmpi1.6 1.6.5-8
ii libpng12-0 1.2.50-1
ii libstdc++6 4.9.0-5
ii libtet1.5 1.5.0-3
ii libx11-6 2:1.6.2-2
ii mpi-default-bin 1.0.2+nmu1
Versions of packages gmsh recommends:
pn gmsh-doc <none>
gmsh suggests no packages.
-- no debconf information
Merge "unit-cube.geo";
Translate{0,0,1}{
Volume{1};
}
Mesh.CharacteristicLengthMin = 1;
Point(1) = {0, 0, 0};
Point(2) = {1, 0, 0};
Point(3) = {1, 1, 0};
Point(4) = {0, 1, 0};
Point(5) = {0, 0, 1};
Point(6) = {1, 0, 1};
Point(7) = {1, 1, 1};
Point(8) = {0, 1, 1};
Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,1};
Line(5) = {5,6};
Line(6) = {6,7};
Line(7) = {7,8};
Line(8) = {8,5};
Line(9) = {1,5};
Line(10) = {2,6};
Line(11) = {3,7};
Line(12) = {4,8};
Line Loop(1) = {1,2,3,4};
Line Loop(2) = {5,6,7,8};
Line Loop(3) = {1,10,-5,-9};
Line Loop(4) = {2,11,-6,-10};
Line Loop(5) = {3,12,-7,-11};
Line Loop(6) = {4,9,-8,-12};
Plane Surface(1) = {1};
Plane Surface(2) = {2};
Plane Surface(3) = {3};
Plane Surface(4) = {4};
Plane Surface(5) = {5};
Plane Surface(6) = {6};
Surface Loop(1) = {1,2,3,4,5,6};
Volume (1) = {1};
$MeshFormat
2.2 0 8
$EndMeshFormat
$Nodes
14
1 0 0 1
2 1 0 1
3 1 1 1
4 0 1 1
5 0 0 2
6 1 0 2
7 1 1 2
8 0 1 2
9 0.5 0.5 0
10 0.5 0.5 1
11 0.5 0 1.5
12 1 0.5 1.5
13 0.5 1 1.5
14 0 0.5 1.5
$EndNodes
$Elements
64
1 15 2 0 1 1
2 15 2 0 2 2
3 15 2 0 3 3
4 15 2 0 4 4
5 15 2 0 5 5
6 15 2 0 6 6
7 15 2 0 7 7
8 15 2 0 8 8
9 1 2 0 1 1 2
10 1 2 0 2 2 3
11 1 2 0 3 3 4
12 1 2 0 4 4 1
13 1 2 0 5 5 6
14 1 2 0 6 6 7
15 1 2 0 7 7 8
16 1 2 0 8 8 5
17 1 2 0 9 1 5
18 1 2 0 10 2 6
19 1 2 0 11 3 7
20 1 2 0 12 4 8
21 2 2 0 1 2 3 9
22 2 2 0 1 1 2 9
23 2 2 0 1 3 4 9
24 2 2 0 1 1 9 4
25 2 2 0 2 6 10 5
26 2 2 0 2 10 8 5
27 2 2 0 2 7 10 6
28 2 2 0 2 8 10 7
29 2 2 0 3 6 11 2
30 2 2 0 3 2 11 1
31 2 2 0 3 11 6 5
32 2 2 0 3 11 5 1
33 2 2 0 4 3 12 2
34 2 2 0 4 12 6 2
35 2 2 0 4 7 12 3
36 2 2 0 4 12 7 6
37 2 2 0 5 4 13 3
38 2 2 0 5 13 7 3
39 2 2 0 5 8 13 4
40 2 2 0 5 13 8 7
41 2 2 0 6 4 14 8
42 2 2 0 6 1 14 4
43 2 2 0 6 5 8 14
44 2 2 0 6 1 5 14
45 4 2 0 1 10 7 12 3
46 4 2 0 1 10 1 5 14
47 4 2 0 1 10 6 2 12
48 4 2 0 1 4 10 8 14
49 4 2 0 1 10 7 3 13
50 4 2 0 1 10 6 11 2
51 4 2 0 1 10 4 8 13
52 4 2 0 1 1 10 5 11
53 4 2 0 1 1 10 2 9
54 4 2 0 1 1 4 10 9
55 4 2 0 1 4 3 10 9
56 4 2 0 1 10 3 2 9
57 4 2 0 1 11 1 10 2
58 4 2 0 1 1 4 14 10
59 4 2 0 1 3 12 10 2
60 4 2 0 1 13 4 3 10
61 4 2 0 1 8 10 13 7
62 4 2 0 1 6 5 10 11
63 4 2 0 1 7 6 10 12
64 4 2 0 1 10 5 8 14
$EndElements
