[Gmsh] Reproduceable gmsh 3D meshing error...

Christophe Geuzaine geuzaine at acm.caltech.edu
Wed Jan 15 02:40:52 CET 2003


"Toshiro K. Ohsumi" wrote:
> 
> Hello,
> 
>    I have a reproducable 3D mesh gmsh error.  Attached below is the model.
> With lc = 0.25, with cross.msh (attached below), the error is attached below.
> With lc = 0.375, no error messages occur.  I am using gmsh 1.36.0.  Many
> thanks,

Yep, it's a bug: the 3D algorithm isn't very robust yet. I'll keep your
file around.

Christophe

> 
>                 - Toshiro K. Ohsumi
> 
> ---------------------------- cross.msh -------------------------------------
> 
> lc = 0.25;
> Point(1) = {-0.35,-17.5,0.1,lc};
> Point(2) = {0.35,-17.5,0.1,lc};
> Point(3) = {0.35,-0.35,0.1,lc};
> Point(4) = {17.5,-0.35,0.1,lc};
> Point(5) = {17.5,0.35,0.1,lc};
> Point(6) = {0.35,0.35,0.1,lc};
> Point(7) = {0.35,17.5,0.1,lc};
> Point(8) = {-0.35,17.5,0.1,lc};
> Point(9) = {-0.35,0.35,0.1,lc};
> Point(10) = {-17.5,0.35,0.1,lc};
> Point(11) = {-17.5,-0.35,0.1,lc};
> Point(12) = {-0.35,-0.35,0.1,lc};
> Point(13) = {-0.35,-17.5,-0.1,lc};
> Point(14) = {0.35,-17.5,-0.1,lc};
> Point(15) = {0.35,-0.35,-0.1,lc};
> Point(16) = {17.5,-0.35,-0.1,lc};
> Point(17) = {17.5,0.35,-0.1,lc};
> Point(18) = {0.35,0.35,-0.1,lc};
> Point(19) = {0.35,17.5,-0.1,lc};
> Point(20) = {-0.35,17.5,-0.1,lc};
> Point(21) = {-0.35,0.35,-0.1,lc};
> Point(22) = {-17.5,0.35,-0.1,lc};
> Point(23) = {-17.5,-0.35,-0.1,lc};
> Point(24) = {-0.35,-0.35,-0.1,lc};
> Line(30) = {1,2};
> Line(31) = {2,3};
> Line(32) = {3,4};
> Line(33) = {4,5};
> Line(34) = {5,6};
> Line(35) = {6,7};
> Line(36) = {7,8};
> Line(37) = {8,9};
> Line(38) = {9,10};
> Line(39) = {10,11};
> Line(40) = {11,12};
> Line(41) = {12,1};
> Line(42) = {13,14};
> Line(43) = {14,15};
> Line(44) = {15,16};
> Line(45) = {16,17};
> Line(46) = {17,18};
> Line(47) = {18,19};
> Line(48) = {19,20};
> Line(49) = {20,21};
> Line(50) = {21,22};
> Line(51) = {22,23};
> Line(52) = {23,24};
> Line(53) = {24,13};
> Line(54) = {1,13};
> Line(55) = {2,14};
> Line(56) = {3,15};
> Line(57) = {4,16};
> Line(58) = {5,17};
> Line(59) = {6,18};
> Line(60) = {7,19};
> Line(61) = {8,20};
> Line(62) = {9,21};
> Line(63) = {10,22};
> Line(64) = {11,23};
> Line(65) = {12,24};
> Line Loop(70) = {30,31,32,33,34,35,36,37,38,39,40,41};
> Plane Surface(71) = {70};
> Line Loop(72) = {42,43,44,45,46,47,48,49,50,51,52,53};
> Plane Surface(73) = {72};
> Line Loop(74) = {30,55,-42,-54};
> Plane Surface(75) = {74};
> Line Loop(76) = {55,43,-56,-31};
> Plane Surface(77) = {76};
> Line Loop(78) = {32,57,-44,-56};
> Plane Surface(79) = {78};
> Line Loop(80) = {57,45,-58,-33};
> Plane Surface(81) = {80};
> Line Loop(82) = {58,46,-59,-34};
> Plane Surface(83) = {82};
> Line Loop(84) = {59,47,-60,-35};
> Plane Surface(85) = {84};
> Line Loop(86) = {60,48,-61,-36};
> Plane Surface(87) = {86};
> Line Loop(88) = {62,-49,-61,37};
> Plane Surface(89) = {88};
> Line Loop(90) = {-38,62,50,-63};
> Plane Surface(91) = {90};
> Line Loop(92) = {63,51,-64,-39};
> Plane Surface(93) = {92};
> Line Loop(94) = {65,-52,-64,40};
> Plane Surface(95) = {94};
> Line Loop(96) = {54,-53,-65,41};
> Plane Surface(97) = {96};
> Surface Loop(100) = {71,73,75,77,79,81,83,85,87,89,91,93,95,97};
> Volume(101) = {100};
> 
> ----------------------------- .gmsh-errors ---------------------------------
> 
> Gmsh 1.36.0
> -------------------------------------------------------
> Build OS       : CYGWIN_ME-4.90 1.3.5(0.47/3/2)
> Build date     : Thu Nov 21 02:10:15  2002
> Build host     : ELAP74
> Packager       : geuzaine
> Home directory : /cygdrive/c/DOCUME~1/TOSHIR~1.OHS/LOCALS~1/Temp/
> Launch date    : Sun Dec 15 22:16:45 2002
> Command line   : /cygdrive/c/Program Files/Gmsh/gmsh
> -------------------------------------------------------
> Geometry
> Loading '/cygdrive/c/TEMP/cross.geo'
> Loaded '/cygdrive/c/TEMP/cross.geo'
> Mesh
> Mesh 1D...
> Meshing Curve 30
> Meshing Curve 31
> Meshing Curve 32
> Meshing Curve 33
> Meshing Curve 34
> Meshing Curve 35
> Meshing Curve 36
> Meshing Curve 37
> Meshing Curve 38
> Meshing Curve 39
> Meshing Curve 40
> Meshing Curve 41
> Meshing Curve 42
> Meshing Curve 43
> Meshing Curve 44
> Meshing Curve 45
> Meshing Curve 46
> Meshing Curve 47
> Meshing Curve 48
> Meshing Curve 49
> Meshing Curve 50
> Meshing Curve 51
> Meshing Curve 52
> Meshing Curve 53
> Meshing Curve 54
> Meshing Curve 55
> Meshing Curve 56
> Meshing Curve 57
> Meshing Curve 58
> Meshing Curve 59
> Meshing Curve 60
> Meshing Curve 61
> Meshing Curve 62
> Meshing Curve 63
> Meshing Curve 64
> Meshing Curve 65
> Mesh 1D complete (0.01 s)
> Mesh 2D...
> Meshing Surface 71
> Meshing Surface 73
> Meshing Surface 75
> Meshing Surface 77
> Meshing Surface 79
> Meshing Surface 81
> Meshing Surface 83
> Meshing Surface 85
> Meshing Surface 87
> Meshing Surface 89
> Meshing Surface 91
> Meshing Surface 93
> Meshing Surface 95
> Meshing Surface 97
> Mesh 2D complete (0.601 s)
> Mesh 3D...
> Mesh 3D... (initial)
> Nod=103/2067 Elm=518
> Vol=672.28 (102 0 1)
> Nod=206/2067 Elm=1026
> Vol=672.28 (205 0 1)
> Nod=309/2067 Elm=1534
> Vol=672.28 (308 0 1)
> Nod=412/2067 Elm=2045
> Vol=672.28 (411 0 1)
> Nod=515/2067 Elm=2656
> Vol=672.28 (514 0 1)
> Nod=618/2067 Elm=3154
> Vol=672.28 (617 0 1)
> Nod=721/2067 Elm=3771
> Vol=672.28 (720 0 1)
> Nod=824/2067 Elm=4496
> Vol=672.28 (797 24 3)
> Nod=927/2067 Elm=4981
> Vol=672.28 (809 110 8)
> Nod=1030/2067 Elm=5617
> Vol=672.28 (823 194 13)
> Nod=1133/2067 Elm=6202
> Vol=672.28 (832 283 18)
> Fatal Error : 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>
> 
> _______________________________________________
> Gmsh mailing list
> Gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh

-- 
Christophe Geuzaine

Tel: (626) 395-4552    http://www.geuz.org
Fax: (626) 578-0124    mailto:geuzaine at acm.caltech.edu