[Gmsh] bug report & question
Magdalena Stolarska
stolarsk at math.umn.edu
Wed Nov 16 18:49:43 CET 2005
Hi. I've come across a fatal error in gmsh that results from a
segmentation violation. The .gmsh-error file and the .geo file that I
used are attached to this email. In order to compute the 3D mesh, I used
the isotropic algorithm. The netgen algorithm was causing gmsh to abort
without any warning and without producing a .gmsh-errors file. Is there
anything wrong with the geo file that would cause this to happen? Any
insight would be very helpful. Thanks in advance for your response.
Magda Stolarska
-------------- next part --------------
h1 = 0.25;
h2 = 0.4;
Point(1) = {0.0,0.0,0.0,h1};
Point(2) = {0.5,0.0,0.0,h1};
Point(3) = {-0.5,0.0,0.0,h1};
Point(4) = {0,0.5,0.0,h1};
Point(5) = {0,-0.5,0.0,h1};
Point(6) = {0,0,5,h2};
Point(7) = {0.5,0,5,h2};
Point(8) = {-0.5,0,5,h2};
Point(9) = {0,0.5,5,h2};
Point(10) = {0,-0.5,5,h2};
Point(11) = {0,0,.2,h1};
Point(12) = {0.2,0,0,h1};
Point(13) = {-0.2,0,0,h1};
Point(14) = {0,0.2,0,h1};
Point(15) = {0,-0.2,0,h1};
Circle(1) = {4,1,3};
Circle(2) = {4,1,2};
Circle(3) = {3,1,5};
Circle(4) = {2,1,5};
Circle(5) = {14,1,13};
Circle(6) = {14,1,12};
Circle(7) = {12,1,15};
Circle(8) = {13,1,15};
Circle(9) = {14,1,11};
Circle(10) = {12,1,11};
Circle(11) = {15,1,11};
Circle(12) = {13,1,11};
Circle(13) = {9,6,8};
Circle(14) = {9,6,7};
Circle(15) = {7,6,10};
Circle(16) = {10,6,8};
Line(17) = {7,2};
Line(18) = {10,5};
Line(19) = {8,3};
Line(20) = {9,4};
Line Loop(21) = {14,15,16,-13};
Plane Surface(22) = {21};
Line Loop(23) = {2,4,-3,-1};
Line Loop(24) = {5,8,-7,-6};
Plane Surface(25) = {23,24};
Plane Surface(26) = {24};
Line Loop(27) = {6,10,-9};
Ruled Surface(28) = {27};
Line Loop(29) = {5,12,-9};
Ruled Surface(30) = {29};
Line Loop(31) = {12,-11,-8};
Ruled Surface(32) = {31};
Line Loop(33) = {11,-10,7};
Ruled Surface(34) = {33};
Line Loop(35) = {17,4,-18,-15};
Ruled Surface(36) = {35};
Line Loop(37) = {17,-2,-20,14};
Ruled Surface(38) = {37};
Line Loop(39) = {13,19,-1,-20};
Ruled Surface(40) = {39};
Line Loop(41) = {16,19,3,-18};
Ruled Surface(42) = {41};
Surface Loop(43) = {26,-30,32,34,28};
//Surface Loop(43) = {-26,30,-32,-34,-28};
Volume(44) = {43};
Surface Loop(45) = {22,-38,36,-25,-42,40,28,34,32,-30};
//Surface Loop(45) = {-22,38,-36,25,42,-40,-28,-34,-32,30};
Volume(46) = {45};
Physical Surface(47) = {22};
Physical Surface(48) = {38,36,40,42};
Physical Surface(49) = {25,26};
Physical Volume(50) = {44};
Physical Volume(51) = {46};
-------------- next part --------------
Info : Plane Surface 22 (4 cross points)
Info : Plane Surface 25 (8 cross points)
Info : Plane Surface 26 (4 cross points)
Info : Mesh 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 13
Info : Meshing curve 14
Info : Meshing curve 15
Info : Meshing curve 16
Info : Meshing curve 17
Info : Meshing curve 18
Info : Meshing curve 19
Info : Meshing curve 20
Info : Mesh 1D complete (0.002999 s)
Info : Mesh 2D...
Info : Meshing surface 22
Info : Meshing surface 25
Info : Meshing surface 26
Info : Meshing surface 28
Info : Meshing surface 30
Info : Meshing surface 32
Info : Meshing surface 34
Info : Meshing surface 36
Info : Meshing surface 38
Info : Meshing surface 40
Info : Meshing surface 42
Info : Mesh 2D complete (0.085987 s)
Info : Mesh 3D...
Info : Mesh 3D... (initial)
Info : Nod=35/713 Elm=160
Info : Vol=13.72 (33 1 1)
Info : Nod=70/713 Elm=357
Info : Vol=13.72 (59 10 1)
Info : Nod=105/713 Elm=570
Info : Vol=13.72 (82 22 1)
Info : Nod=140/713 Elm=799
Info : Vol=13.72 (108 31 1)
Info : Nod=175/713 Elm=1061
Info : Vol=13.72 (122 52 1)
Info : Nod=210/713 Elm=1301
Info : Vol=13.72 (146 63 1)
Info : Nod=245/713 Elm=1540
Info : Vol=13.72 (164 80 1)
Info : Nod=280/713 Elm=1797
Info : Vol=13.72 (187 92 1)
Info : Nod=315/713 Elm=2027
Info : Vol=13.72 (215 99 1)
Info : Nod=350/713 Elm=2281
Info : Vol=13.72 (250 99 1)
Info : Nod=385/713 Elm=2564
Info : Vol=13.72 (284 100 1)
Info : Nod=420/713 Elm=2764
Info : Vol=13.72 (306 113 1)
Info : Nod=455/713 Elm=3078
Info : Vol=13.72 (330 124 1)
Info : Nod=490/713 Elm=3319
Info : Vol=13.72 (350 139 1)
Info : Nod=525/713 Elm=3620
Info : Vol=13.72 (374 150 1)
Info : Nod=560/713 Elm=3863
Info : Vol=13.72 (394 165 1)
Info : Nod=595/713 Elm=4159
Info : Vol=13.72 (416 178 1)
Info : Nod=630/713 Elm=4465
Info : Vol=13.72 (441 188 1)
Info : Nod=665/713 Elm=4792
Info : Vol=13.72 (466 198 1)
Info : Nod=700/713 Elm=5048
Info : Vol=13.72 (499 200 1)
Info : Boundary edges recovery
Info : Boundary faces recovery
Info : Volume recovery
Info : Surface 22
Info : Surface 25
Info : Surface 26
Info : Surface 36
Info : Surface 38
Info : Surface 40
Info : Surface 42
Info : Surface 22
Info : Surface 25
Info : Surface 28
Info : Surface 30
Info : Surface 32
Info : Surface 34
Info : Surface 36
Info : Surface 38
Info : Surface 40
Info : Surface 42
Info : Surface 26
Info : Surface 28
Info : Surface 30
Info : Surface 32
Info : Surface 34
Info : Mesh 3D... (final)
Info : Nod=735 Elm=1963
Info : Vol(3.82351) Conv(2.12209->0.8)
Info : Nod=759 Elm=2140
Info : Vol(3.82351) Conv(1.844->0.8)
Info : Nod=783 Elm=2318
Info : Vol(3.82351) Conv(1.6142->0.8)
Info : Nod=807 Elm=2498
Info : Vol(3.82351) Conv(1.52943->0.8)
Info : Nod=831 Elm=2683
Info : Vol(3.82351) Conv(1.43614->0.8)
Info : Nod=855 Elm=2844
Info : Vol(3.82351) Conv(1.36667->0.8)
Info : Nod=879 Elm=3028
Info : Vol(3.82351) Conv(1.32745->0.8)
Info : Nod=903 Elm=3193
Info : Vol(3.82351) Conv(1.29444->0.8)
Info : Nod=927 Elm=3371
Info : Vol(3.82351) Conv(1.25236->0.8)
Info : Nod=951 Elm=3539
Info : Vol(3.82351) Conv(1.19972->0.8)
Info : Nod=975 Elm=3715
Info : Vol(3.82351) Conv(1.17294->0.8)
Info : Nod=999 Elm=3883
Info : Vol(3.82351) Conv(1.1521->0.8)
Info : Nod=1023 Elm=4049
Info : Vol(3.82351) Conv(1.12663->0.8)
Info : Nod=1047 Elm=4189
Info : Vol(3.82351) Conv(1.10436->0.8)
Info : Nod=1071 Elm=4357
Info : Vol(3.82351) Conv(1.08356->0.8)
Info : Nod=1095 Elm=4503
Info : Vol(3.82351) Conv(1.069->0.8)
Info : Nod=1119 Elm=4637
Info : Vol(3.82351) Conv(1.05405->0.8)
Info : Nod=1143 Elm=4782
Info : Vol(3.82351) Conv(1.03843->0.8)
Info : Nod=1167 Elm=4929
Info : Vol(3.82351) Conv(1.02492->0.8)
Info : Nod=1191 Elm=5096
Info : Vol(3.82351) Conv(1.00304->0.8)
Info : Nod=1215 Elm=5246
Info : Vol(3.82351) Conv(0.987187->0.8)
Info : Nod=1239 Elm=5404
Info : Vol(3.82351) Conv(0.975215->0.8)
Info : Nod=1263 Elm=5553
Info : Vol(3.82351) Conv(0.962466->0.8)
Info : Nod=1287 Elm=5720
Info : Vol(3.82351) Conv(0.952105->0.8)
Info : Nod=1311 Elm=5872
Info : Vol(3.82351) Conv(0.942153->0.8)
Info : Nod=1335 Elm=6017
Info : Vol(3.82351) Conv(0.932487->0.8)
Info : Nod=1359 Elm=6160
Info : Vol(3.82351) Conv(0.92357->0.8)
Info : Nod=1383 Elm=6314
Info : Vol(3.82351) Conv(0.909862->0.8)
Info : Nod=1407 Elm=6482
Info : Vol(3.82351) Conv(0.903469->0.8)
Info : Nod=1431 Elm=6632
Info : Vol(3.82351) Conv(0.897597->0.8)
Info : Nod=1455 Elm=6771
Info : Vol(3.82351) Conv(0.891553->0.8)
Info : Nod=1478 Elm=6928
Info : Vol(3.82351) Conv(0.885092->0.8)
Info : Nod=1502 Elm=7091
Info : Vol(3.82351) Conv(0.878992->0.8)
Info : Nod=1526 Elm=7241
Info : Vol(3.82351) Conv(0.873392->0.8)
Info : Nod=1550 Elm=7414
Info : Vol(3.82351) Conv(0.868922->0.8)
Info : Nod=1573 Elm=7551
Info : Vol(3.82351) Conv(0.863861->0.8)
Info : Nod=1597 Elm=7705
Info : Vol(3.82351) Conv(0.857172->0.8)
Info : Nod=1621 Elm=7848
Info : Vol(3.82351) Conv(0.851277->0.8)
Info : Nod=1642 Elm=7991
Info : Vol(3.82351) Conv(0.845036->0.8)
Info : Nod=1664 Elm=8131
Info : Vol(3.82351) Conv(0.839822->0.8)
Info : Nod=1688 Elm=8267
Info : Vol(3.82351) Conv(0.835383->0.8)
Info : Nod=1712 Elm=8417
Info : Vol(3.82351) Conv(0.830876->0.8)
Info : Nod=1736 Elm=8571
Info : Vol(3.82351) Conv(0.823484->0.8)
Info : Nod=1759 Elm=8725
Info : Vol(3.82351) Conv(0.819886->0.8)
Info : Nod=1781 Elm=8866
Info : Vol(3.82351) Conv(0.81363->0.8)
Info : Nod=1802 Elm=9004
Info : Vol(3.82351) Conv(0.808183->0.8)
Info : Nod=1825 Elm=9148
Info : Vol(3.82351) Conv(0.802846->0.8)
Info : Nod=1847 Elm=9287
Info : Vol(3.82351) Conv(0.800283->0.8)
Info : Mesh 3D complete (1.60576 s)
Info : Optimize volume mesh...
Info : Optimizing volume 44
Info : Optimizing volume 46
Info : Optimize volume mesh complete (6.36703 s)
Info : Plane Surface 22 (4 cross points)
Info : Plane Surface 25 (8 cross points)
Info : Plane Surface 26 (4 cross points)
Info : Mesh 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 13
Info : Meshing curve 14
Info : Meshing curve 15
Info : Meshing curve 16
Info : Meshing curve 17
Info : Meshing curve 18
Info : Meshing curve 19
Info : Meshing curve 20
Info : Mesh 1D complete (0.002 s)
Info : Mesh 2D...
Info : Meshing surface 22
Info : Meshing surface 25
Info : Meshing surface 26
Info : Meshing surface 28
Info : Meshing surface 30
Info : Meshing surface 32
Info : Meshing surface 34
Info : Meshing surface 36
Info : Meshing surface 38
Info : Meshing surface 40
Info : Meshing surface 42
Info : Mesh 2D complete (0.020997 s)
Info : Mesh 3D...
Info : Mesh 3D... (initial)
Info : Nod=7/158 Elm=32
Info : Vol=13.72 (6 0 1)
Info : Nod=14/158 Elm=60
Info : Vol=13.72 (13 0 1)
Info : Nod=21/158 Elm=99
Info : Vol=13.72 (19 1 1)
Info : Nod=28/158 Elm=131
Info : Vol=13.72 (26 1 1)
Info : Nod=35/158 Elm=169
Info : Vol=13.72 (30 4 1)
Info : Nod=42/158 Elm=211
Info : Vol=13.72 (34 7 1)
Info : Nod=49/158 Elm=255
Info : Vol=13.72 (39 9 1)
Info : Nod=56/158 Elm=303
Info : Vol=13.72 (46 9 1)
Info : Nod=63/158 Elm=341
Info : Vol=13.72 (53 9 1)
Info : Nod=70/158 Elm=393
Info : Vol=13.72 (60 9 1)
Info : Nod=77/158 Elm=433
Info : Vol=13.72 (66 10 1)
Info : Nod=84/158 Elm=475
Info : Vol=13.72 (73 10 1)
Info : Nod=91/158 Elm=522
Info : Vol=13.72 (80 10 1)
Info : Nod=98/158 Elm=568
Info : Vol=13.72 (87 10 1)
Info : Nod=105/158 Elm=626
Info : Vol=13.72 (93 11 1)
Info : Nod=112/158 Elm=671
Info : Vol=13.72 (100 11 1)
Info : Nod=119/158 Elm=712
Info : Vol=13.72 (104 14 1)
Info : Nod=126/158 Elm=762
Info : Vol=13.72 (109 16 1)
Info : Nod=133/158 Elm=818
Info : Vol=13.72 (114 18 1)
Info : Nod=140/158 Elm=874
Info : Vol=13.72 (121 18 1)
Info : Nod=147/158 Elm=918
Info : Vol=13.72 (128 18 1)
Info : Nod=154/158 Elm=968
Info : Vol=13.72 (135 18 1)
Info : Boundary edges recovery
Info : Boundary faces recovery
Info : Volume recovery
Info : Surface 22
Info : Surface 25
Info : Surface 26
Info : Surface 36
Info : Surface 38
Info : Surface 40
Info : Surface 42
Info : Surface 26
Info : Surface 28
Info : Surface 30
Info : Surface 32
Info : Surface 34
Info : Surface 22
Info : Surface 25
Info : Surface 28
Info : Surface 30
Info : Surface 32
Info : Surface 34
Info : Surface 36
Info : Surface 38
Info : Surface 40
Info : Surface 42
Info : Mesh 3D... (final)
Info : Nod=162 Elm=423
Info : Vol(3.57464) Conv(1.52652->0.8)
Info : Nod=168 Elm=435
Info : Vol(3.57464) Conv(1.187->0.8)
Info : Nod=174 Elm=485
Info : Vol(3.57464) Conv(1.11365->0.8)
Info : Nod=180 Elm=540
Info : Vol(3.57464) Conv(1.03598->0.8)
Info : Nod=186 Elm=578
Info : Vol(3.57464) Conv(1.00313->0.8)
Info : Nod=192 Elm=622
Info : Vol(3.57464) Conv(0.968589->0.8)
Info : Nod=198 Elm=665
Info : Vol(3.57464) Conv(0.95392->0.8)
Info : Nod=204 Elm=707
Info : Vol(3.57464) Conv(0.937352->0.8)
Info : Nod=210 Elm=734
Info : Vol(3.57464) Conv(0.923661->0.8)
Info : Nod=216 Elm=783
Info : Vol(3.57464) Conv(0.91226->0.8)
Info : Nod=222 Elm=813
Info : Vol(3.57464) Conv(0.887802->0.8)
Info : Nod=228 Elm=842
Info : Vol(3.57464) Conv(0.870903->0.8)
Info : Nod=234 Elm=877
Info : Vol(3.57464) Conv(0.851815->0.8)
Info : Nod=240 Elm=908
Info : Vol(3.57464) Conv(0.843642->0.8)
Info : Nod=246 Elm=938
Info : Vol(3.57464) Conv(0.828406->0.8)
Info : Nod=252 Elm=974
Info : Vol(3.57464) Conv(0.819225->0.8)
Info : Nod=254 Elm=986
Info : Vol(3.57464) Conv(0.806347->0.8)
Info : Mesh 3D complete (0.201969 s)
Info : Optimize volume mesh...
Info : Optimizing volume 44
Fatal : Segmentation violation (invalid memory reference)
: ------------------------------------------------------
: You have discovered a bug in Gmsh! You may report it
: by e-mail (together with any helpful data permitting to
: reproduce it) to <gmsh at geuz.org>