[Gmsh] Error transfinite volume

Amit Itagi aitagi at gmail.com
Tue Nov 4 21:52:11 CET 2008


Hi,

I am getting an error while meshing a transfinite volume. When I comment out
the transfinite volume lines, there is no error. Could you kindly take a
look ?

Here is my "geo" file.

Point(1) = {-1,-1,0,0.2};
Point(2) = {1,-1,0,0.2};
Point(3) = {1,1,0,0.2};
Point(4) = {-1,1,0,0.2};
Line(1) = {1,2};
Line(2) = {2,3};
Line(3) = {3,4};
Line(4) = {4,1};
Line Loop(5) = {4,1,2,3};
Plane Surface(6) = {5};
Point(5) = {-5,-5,0,1};
Point(6) = {5,-5,0,1};
Point(7) = {5,5,0,1};
Point(8) = {-5,5,0,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(13) = {1,10,-5,-9};
Plane Surface(14) = {13};
Line Loop(15) = {2,11,-6,-10};
Plane Surface(16) = {15};
Line Loop(17) = {3,12,-7,-11};
Plane Surface(18) = {17};
Line Loop(19) = {4,9,-8,-12};
Plane Surface(20) = {19};
Point(101) = {-1,-1,10,0.2};
Point(102) = {1,-1,10,0.2};
Point(103) = {1,1,10,0.2};
Point(104) = {-1,1,10,0.2};
Line(101) = {101,102};
Line(102) = {102,103};
Line(103) = {103,104};
Line(104) = {104,101};
Line Loop(105) = {104,101,102,103};
Plane Surface(106) = {105};
Point(105) = {-5,-5,10,1};
Point(106) = {5,-5,10,1};
Point(107) = {5,5,10,1};
Point(108) = {-5,5,10,1};
Line(105) = {105,106};
Line(106) = {106,107};
Line(107) = {107,108};
Line(108) = {108,105};
Line(109) = {101,105};
Line(110) = {102,106};
Line(111) = {103,107};
Line(112) = {104,108};
Line Loop(113) = {101,110,-105,-109};
Plane Surface(114) = {113};
Line Loop(115) = {102,111,-106,-110};
Plane Surface(116) = {115};
Line Loop(117) = {103,112,-107,-111};
Plane Surface(118) = {117};
Line Loop(119) = {104,109,-108,-112};
Plane Surface(120) = {119};
Line(200) = {1,101};
Line(201) = {2,102};
Line(202) = {3,103};
Line(203) = {4,104};
Line(204) = {5,105};
Line(205) = {6,106};
Line(206) = {7,107};
Line(207) = {8,108};
Line Loop(208) = {1,201,-101,-200};
Plane Surface(209) = {208};
Line Loop(210) = {2,202,-102,-201};
Plane Surface(211) = {210};
Line Loop(212) = {3,203,-103,-202};
Plane Surface(213) = {212};
Line Loop(214) = {4,200,-104,-203};
Plane Surface(215) = {214};
Line Loop(216) = {9,204,-109,-200};
Plane Surface(217) = {216};
Line Loop(218) = {10,205,-110,-201};
Plane Surface(219) = {218};
Line Loop(220) = {11,206,-111,-202};
Plane Surface(221) = {220};
Line Loop(222) = {12,207,-112,-203};
Plane Surface(223) = {222};
Line Loop(224) = {5,205,-105,-204};
Plane Surface(225) = {224};
Line Loop(226) = {6,206,-106,-205};
Plane Surface(227) = {226};
Line Loop(228) = {7,207,-107,-206};
Plane Surface(229) = {228};
Line Loop(230) = {8,204,-108,-207};
Plane Surface(231) = {230};
Surface Loop(232) = {209,211,213,215,6,106};
Volume(233) = {232};
Surface Loop(234) = {14,225,114,217,219,209};
Volume(235) = {234};
Surface Loop(236) = {227,16,116,219,221,211};
Volume(237) = {236};
Surface Loop(238) = {18,229,118,223,221,213};
Volume(239) = {238};
Surface Loop(240) = {20,231,120,217,223,215};
Volume(241) = {240};

Transfinite Line{1,2,3,4}=5;
Transfinite Line{5,6,7,8}=5;
Transfinite Line{101,102,103,104}=5;
Transfinite Line{9,10,11,12}=5 Using Progression 1.2;
Transfinite Line{105,106,107,108}=5;
Transfinite Line{109,110,111,112}=5 Using Progression 1.2;
Transfinite Line{200:207}=5;

Transfinite Surface{6}={1,2,3,4};
Transfinite Surface{14}={1,2,6,5};
Transfinite Surface{16}={2,3,7,6};
Transfinite Surface{18}={3,4,8,7};
Transfinite Surface{20}={4,1,5,8};
Transfinite Surface{106}={101,102,103,104};
Transfinite Surface{114}={101,102,106,105};
Transfinite Surface{116}={102,103,107,106};
Transfinite Surface{118}={103,104,108,107};
Transfinite Surface{120}={104,101,105,108};
Transfinite Surface{209}={1,2,102,101};
Transfinite Surface{211}={2,3,103,102};
Transfinite Surface{213}={3,4,104,103};
Transfinite Surface{215}={4,1,101,104};
Transfinite Surface{217}={1,5,105,101};
Transfinite Surface{219}={2,6,106,102};
Transfinite Surface{221}={3,7,107,103};
Transfinite Surface{223}={4,8,108,104};
Transfinite Surface{225}={5,6,106,105};
Transfinite Surface{227}={6,7,107,106};
Transfinite Surface{229}={7,8,108,107};
Transfinite Surface{231}={8,5,105,108};

Recombine Surface{6};
Recombine Surface{14};
Recombine Surface{16};
Recombine Surface{18};
Recombine Surface{20};
Recombine Surface{106};
Recombine Surface{114};
Recombine Surface{116};
Recombine Surface{118};
Recombine Surface{120};
Recombine Surface{209};
Recombine Surface{211};
Recombine Surface{213};
Recombine Surface{215};
Recombine Surface{217};
Recombine Surface{219};
Recombine Surface{221};
Recombine Surface{223};
Recombine Surface{225};
Recombine Surface{227};
Recombine Surface{229};Info    :
-------------------------------------------------------
Info    : Gmsh version   : 2.0.7
Info    : Build OS       : Linux
Info    : Build options  : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB MATHEVAL
ANN OCC
Info    : Build date     : Mon Apr  2 11:09:24 CEST 2007
Info    : Build host     : debian
Info    : Packager       : geuzaine
Info    : Home directory : /home/amit/
Info    : Launch date    : Tue Nov  4 15:50:00 2008
Info    : Command line   : gmsh test.geo
Info    : -------------------------------------------------------
Info    : Reading 'test.geo'
Info    : Read 'test.geo'
Info    : Meshing 1D...
Info    : Meshing curve 1
Info    : Meshing curve 2
Info    : Meshing curve 3
Info    : Meshing curve 4
Info    : Meshing curve 5
Info    : Meshing curve 6
Info    : Meshing curve 7
Info    : Meshing curve 8
Info    : Meshing curve 9
Info    : Meshing curve 10
Info    : Meshing curve 11
Info    : Meshing curve 12
Info    : Meshing curve 101
Info    : Meshing curve 102
Info    : Meshing curve 103
Info    : Meshing curve 104
Info    : Meshing curve 105
Info    : Meshing curve 106
Info    : Meshing curve 107
Info    : Meshing curve 108
Info    : Meshing curve 109
Info    : Meshing curve 110
Info    : Meshing curve 111
Info    : Meshing curve 112
Info    : Meshing curve 200
Info    : Meshing curve 201
Info    : Meshing curve 202
Info    : Meshing curve 203
Info    : Meshing curve 204
Info    : Meshing curve 205
Info    : Meshing curve 206
Info    : Meshing curve 207
Info    : Mesh 1D complete (0.008001 s)
Info    : Mesh
Info    : Meshing 2D...
Info    : Meshing surface 6 (Plane)
Info    : Meshing surface 14 (Plane)
Info    : Meshing surface 16 (Plane)
Info    : Meshing surface 18 (Plane)
Info    : Meshing surface 20 (Plane)
Info    : Meshing surface 106 (Plane)
Info    : Meshing surface 114 (Plane)
Info    : Meshing surface 116 (Plane)
Info    : Meshing surface 118 (Plane)
Info    : Meshing surface 120 (Plane)
Info    : Meshing surface 209 (Plane)
Info    : Meshing surface 211 (Plane)
Info    : Meshing surface 213 (Plane)
Info    : Meshing surface 215 (Plane)
Info    : Meshing surface 217 (Plane)
Info    : Meshing surface 219 (Plane)
Info    : Meshing surface 221 (Plane)
Info    : Meshing surface 223 (Plane)
Info    : Meshing surface 225 (Plane)
Info    : Meshing surface 227 (Plane)
Info    : Meshing surface 229 (Plane)
Info    : Meshing surface 231 (Plane)
Info    : Mesh 2D complete (0 s)
Info    : Mesh
Info    : Meshing 3D...
Info    : Meshing volume 233
Info    : Meshing volume 235
Error   : Incompatible surface 14 in transfinite volume 235
Info    : Meshing volume 237
Error   : Incompatible surface 16 in transfinite volume 237
Info    : Meshing volume 239
Error   : Incompatible surface 18 in transfinite volume 239
Info    : Meshing volume 241
Error   : Incompatible surface 20 in transfinite volume 241
Info    : Mesh 3D complete (0 s)
Info    : Mesh
Info    : 337 vertices 544 elements
Info    : Writing '/home/amit/programs/FEM/myExamples/DTR/error.txt'
Recombine Surface{231};

Transfinite Volume{233}={1,2,3,4,101,102,103,104};
Transfinite Volume{235}={1,2,5,6,101,102,105,106};
Transfinite Volume{237}={2,3,6,7,102,103,106,107};
Transfinite Volume{239}={3,4,7,8,103,104,107,108};
Transfinite Volume{241}={4,1,8,5,104,101,108,105};


I am using an older version of gmsh. Here is the error log.
Info    : -------------------------------------------------------
Info    : Gmsh version   : 2.0.7
Info    : Build OS       : Linux
Info    : Build options  : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB MATHEVAL
ANN OCC
Info    : Build date     : Mon Apr  2 11:09:24 CEST 2007
Info    : Build host     : debian
Info    : Packager       : geuzaine
Info    : Home directory : /home/amit/
Info    : Launch date    : Tue Nov  4 15:50:00 2008
Info    : Command line   : gmsh test.geo
Info    : -------------------------------------------------------
Info    : Reading 'test.geo'
Info    : Read 'test.geo'
Info    : Meshing 1D...
Info    : Meshing curve 1
Info    : Meshing curve 2
Info    : Meshing curve 3
Info    : Meshing curve 4
Info    : Meshing curve 5
Info    : Meshing curve 6
Info    : Meshing curve 7
Info    : Meshing curve 8
Info    : Meshing curve 9
Info    : Meshing curve 10
Info    : Meshing curve 11
Info    : Meshing curve 12
Info    : Meshing curve 101
Info    : Meshing curve 102
Info    : Meshing curve 103
Info    : Meshing curve 104
Info    : Meshing curve 105
Info    : Meshing curve 106
Info    : Meshing curve 107
Info    : Meshing curve 108
Info    : Meshing curve 109
Info    : Meshing curve 110
Info    : Meshing curve 111
Info    : Meshing curve 112
Info    : Meshing curve 200
Info    : Meshing curve 201
Info    : Meshing curve 202
Info    : Meshing curve 203
Info    : Meshing curve 204
Info    : Meshing curve 205
Info    : Meshing curve 206Info    :
-------------------------------------------------------
Info    : Gmsh version   : 2.0.7Info    :
-------------------------------------------------------
Info    : Gmsh version   : 2.0.7
Info    : Build OS       : Linux
Info    : Build options  : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB MATHEVAL
ANN OCC
Info    : Build date     : Mon Apr  2 11:09:24 CEST 2007
Info    : Build host     : debian
Info    : Packager       : geuzaine
Info    : Home directory : /home/amit/
Info    : Launch date    : Tue Nov  4 15:50:00 2008
Info    : Command line   : gmsh test.geo
Info    : -------------------------------------------------------
Info    : Reading 'test.geo'
Info    : Read 'test.geo'
Info    : Meshing 1D...
Info    : Meshing curve 1
Info    : Meshing curve 2
Info    : Meshing curve 3
Info    : Meshing curve 4
Info    : Meshing curve 5
Info    : Meshing curve 6
Info    : Meshing curve 7
Info    : Meshing curve 8
Info    : Meshing curve 9
Info    : Meshing curve 10
Info    : Meshing curve 11
Info    : Meshing curve 12
Info    : Meshing curve 101
Info    : Meshing curve 102
Info    : Meshing curve 103
Info    : Meshing curve 104
Info    : Meshing curve 105
Info    : Meshing curve 106
Info    : Meshing curve 107
Info    : Meshing curve 108
Info    : Meshing curve 109
Info    : Meshing curve 110
Info    : Meshing curve 111
Info    : Meshing curve 112
Info    : Meshing curve 200
Info    : Meshing curve 201
Info    : Meshing curve 202
Info    : Meshing curve 203
Info    : Meshing curve 204
Info    : Meshing curve 205
Info    : Meshing curve 206
Info    : Meshing curve 207
Info    : Mesh 1D complete (0.008001 s)
Info    : Mesh
Info    : Meshing 2D...
Info    : Meshing surface 6 (Plane)
Info    : Meshing surface 14 (Plane)
Info    : Meshing surface 16 (Plane)
Info    : Meshing surface 18 (Plane)
Info    : Meshing surface 20 (Plane)
Info    : Meshing surface 106 (Plane)
Info    : Meshing surface 114 (Plane)
Info    : Meshing surface 116 (Plane)
Info    : Meshing surface 118 (Plane)
Info    : Meshing surface 120 (Plane)
Info    : Meshing surface 209 (Plane)
Info    : Meshing surface 211 (Plane)
Info    : Meshing surface 213 (Plane)
Info    : Meshing surface 215 (Plane)
Info    : Meshing surface 217 (Plane)
Info    : Meshing surface 219 (Plane)
Info    : Meshing surface 221 (Plane)
Info    : Meshing surface 223 (Plane)
Info    : Meshing surface 225 (Plane)
Info    : Meshing surface 227 (Plane)
Info    : Meshing surface 229 (Plane)
Info    : Meshing surface 231 (Plane)
Info    : Mesh 2D complete (0 s)
Info    : Mesh
Info    : Meshing 3D...
Info    : Meshing volume 233
Info    : Meshing volume 235
Error   : Incompatible surface 14 in transfinite volume 235
Info    : Meshing volume 237
Error   : Incompatible surface 16 in transfinite volume 237
Info    : Meshing volume 239
Error   : Incompatible surface 18 in transfinite volume 239
Info    : Meshing volume 241
Error   : Incompatible surface 20 in transfinite volume 241
Info    : Mesh 3D complete (0 s)
Info    : Mesh
Info    : 337 vertices 544 elements
Info    : Writing '/home/amit/programs/FEM/myExamples/DTR/error.txt'
Info    : Build OS       : Linux
Info    : Build options  : GSL TRIANGLE NETGEN TETGEN JPEG PNG ZLIB MATHEVAL
ANN OCC
Info    : Build date     : Mon Apr  2 11:09:24 CEST 2007
Info    : Build host     : debian
Info    : Packager       : geuzaine
Info    : Home directory : /home/amit/
Info    : Launch date    : Tue Nov  4 15:50:00 2008
Info    : Command line   : gmsh test.geo
Info    : -------------------------------------------------------
Info    : Reading 'test.geo'
Info    : Read 'test.geo'
Info    : Meshing 1D...
Info    : Meshing curve 1
Info    : Meshing curve 2
Info    : Meshing curve 3
Info    : Meshing curve 4
Info    : Meshing curve 5
Info    : Meshing curve 6
Info    : Meshing curve 7
Info    : Meshing curve 8
Info    : Meshing curve 9
Info    : Meshing curve 10
Info    : Meshing curve 11
Info    : Meshing curve 12
Info    : Meshing curve 101
Info    : Meshing curve 102
Info    : Meshing curve 103
Info    : Meshing curve 104
Info    : Meshing curve 105
Info    : Meshing curve 106
Info    : Meshing curve 107
Info    : Meshing curve 108
Info    : Meshing curve 109
Info    : Meshing curve 110
Info    : Meshing curve 111
Info    : Meshing curve 112
Info    : Meshing curve 200
Info    : Meshing curve 201
Info    : Meshing curve 202
Info    : Meshing curve 203
Info    : Meshing curve 204
Info    : Meshing curve 205
Info    : Meshing curve 206
Info    : Meshing curve 207
Info    : Mesh 1D complete (0.008001 s)
Info    : Mesh
Info    : Meshing 2D...
Info    : Meshing surface 6 (Plane)
Info    : Meshing surface 14 (Plane)
Info    : Meshing surface 16 (Plane)
Info    : Meshing surface 18 (Plane)
Info    : Meshing surface 20 (Plane)
Info    : Meshing surface 106 (Plane)
Info    : Meshing surface 114 (Plane)
Info    : Meshing surface 116 (Plane)
Info    : Meshing surface 118 (Plane)
Info    : Meshing surface 120 (Plane)
Info    : Meshing surface 209 (Plane)
Info    : Meshing surface 211 (Plane)
Info    : Meshing surface 213 (Plane)
Info    : Meshing surface 215 (Plane)
Info    : Meshing surface 217 (Plane)
Info    : Meshing surface 219 (Plane)
Info    : Meshing surface 221 (Plane)
Info    : Meshing surface 223 (Plane)
Info    : Meshing surface 225 (Plane)
Info    : Meshing surface 227 (Plane)
Info    : Meshing surface 229 (Plane)
Info    : Meshing surface 231 (Plane)
Info    : Mesh 2D complete (0 s)
Info    : Mesh
Info    : Meshing 3D...
Info    : Meshing volume 233
Info    : Meshing volume 235
Error   : Incompatible surface 14 in transfinite volume 235
Info    : Meshing volume 237
Error   : Incompatible surface 16 in transfinite volume 237
Info    : Meshing volume 239
Error   : Incompatible surface 18 in transfinite volume 239
Info    : Meshing volume 241
Error   : Incompatible surface 20 in transfinite volume 241
Info    : Mesh 3D complete (0 s)
Info    : Mesh
Info    : 337 vertices 544 elements
Info    : Writing '/home/amit/programs/FEM/myExamples/DTR/error.txt'
Info    : Meshing curve 207
Info    : Mesh 1D complete (0.008001 s)
Info    : Mesh
Info    : Meshing 2D...
Info    : Meshing surface 6 (Plane)
Info    : Meshing surface 14 (Plane)
Info    : Meshing surface 16 (Plane)
Info    : Meshing surface 18 (Plane)
Info    : Meshing surface 20 (Plane)
Info    : Meshing surface 106 (Plane)
Info    : Meshing surface 114 (Plane)
Info    : Meshing surface 116 (Plane)
Info    : Meshing surface 118 (Plane)
Info    : Meshing surface 120 (Plane)
Info    : Meshing surface 209 (Plane)
Info    : Meshing surface 211 (Plane)
Info    : Meshing surface 213 (Plane)
Info    : Meshing surface 215 (Plane)
Info    : Meshing surface 217 (Plane)
Info    : Meshing surface 219 (Plane)
Info    : Meshing surface 221 (Plane)
Info    : Meshing surface 223 (Plane)
Info    : Meshing surface 225 (Plane)
Info    : Meshing surface 227 (Plane)
Info    : Meshing surface 229 (Plane)
Info    : Meshing surface 231 (Plane)
Info    : Mesh 2D complete (0 s)
Info    : Mesh
Info    : Meshing 3D...
Info    : Meshing volume 233
Info    : Meshing volume 235
Error   : Incompatible surface 14 in transfinite volume 235
Info    : Meshing volume 237
Error   : Incompatible surface 16 in transfinite volume 237
Info    : Meshing volume 239
Error   : Incompatible surface 18 in transfinite volume 239
Info    : Meshing volume 241
Error   : Incompatible surface 20 in transfinite volume 241
Info    : Mesh 3D complete (0 s)
Info    : Mesh
Info    : 337 vertices 544 elements



Thanks

Rgds,
Amit
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20081104/0f3b1f3e/attachment.html>