[Gmsh] Problems with embedded crack

Fri Oct 15 17:55:50 CEST 2010

On 29/09/10 18:21, olivier jamond wrote:
> It seems to be okay for this case but not for every one. Here is one
> which seems to be much more problematic (sometimes it works but if you
> try to mesh it successive times, it will fail):

Hi Olivier - there's indeed a nasty bug in the implementation of 
embedded lines. This is now in our bug tracking database: 
https://geuz.org/trac/gmsh/ticket/96

>
> Point(1) = {80.000000,24.559024,0};
> Point(2) = {80.000000,28.000000,0};
> Point(3) = {80.000000,31.206362,0};
> Point(4) = {79.983338,31.234574,0};
> Point(5) = {79.551972,31.799626,0};
> Point(6) = {79.056802,32.309693,0};
> Point(7) = {78.504781,32.757614,0};
> Point(8) = {77.903658,33.137103,0};
> Point(9) = {77.261870,33.442831,0};
> Point(10) = {76.588427,33.670507,0};
> Point(11) = {75.892784,33.816934,0};
> Point(12) = {75.184704,33.880058,0};
> Point(13) = {74.474129,33.858992,0};
> Point(14) = {73.771033,33.754032,0};
> Point(15) = {73.085285,33.566651,0};
> Point(16) = {72.426513,33.299480,0};
> Point(17) = {71.803964,32.956269,0};
> Point(18) = {71.226377,32.541836,0};
> Point(19) = {70.701861,32.061999,0};
> Point(20) = {70.237778,31.523494,0};
> Point(21) = {69.840643,30.933880,0};
> Point(22) = {69.516031,30.301434,0};
> Point(23) = {69.268498,29.635034,0};
> Point(24) = {69.101521,28.944035,0};
> Point(25) = {69.017442,28.238137,0};
> Point(26) = {69.017442,27.527249,0};
> Point(27) = {69.101521,26.821351,0};
> Point(28) = {69.268498,26.130352,0};
> Point(29) = {69.516031,25.463952,0};
> Point(30) = {69.840643,24.831506,0};
> Point(31) = {70.237778,24.241892,0};
> Point(32) = {70.701861,23.703387,0};
> Point(33) = {71.226377,23.223550,0};
> Point(34) = {71.803964,22.809117,0};
> Point(35) = {72.426513,22.465906,0};
> Point(36) = {73.085285,22.198735,0};
> Point(37) = {73.771033,22.011354,0};
> Point(38) = {74.474129,21.906394,0};
> Point(39) = {75.184704,21.885328,0};
> Point(40) = {75.892784,21.948452,0};
> Point(41) = {76.588427,22.094880,0};
> Point(42) = {77.261870,22.322555,0};
> Point(43) = {77.903658,22.628283,0};
> Point(44) = {78.504781,23.007772,0};
> Point(45) = {79.056802,23.455693,0};
> Point(46) = {79.551972,23.965760,0};
> Point(47) = {79.983338,24.530812,0};
> Point(48) = {78.504938,28.000000,0};
> Point(49) = {78.504938,28.000000,0};
> Point(50) = {78.504938,28.000000,0};
> Point(51) = {78.466511,28.412900,0};
> Point(52) = {78.379519,28.818357,0};
> Point(53) = {78.245183,29.210680,0};
> Point(54) = {78.065389,29.584360,0};
> Point(55) = {77.842660,29.934153,0};
> Point(56) = {77.580123,30.255148,0};
> Point(57) = {77.281465,30.542839,0};
> Point(58) = {76.950876,30.793187,0};
> Point(59) = {76.592998,31.002678,0};
> Point(60) = {76.212855,31.168372,0};
> Point(61) = {75.815783,31.287941,0};
> Point(62) = {75.407356,31.359709,0};
> Point(63) = {74.993307,31.382667,0};
> Point(64) = {74.579450,31.356492,0};
> Point(65) = {74.171593,31.281554,0};
> Point(66) = {73.775462,31.158903,0};
> Point(67) = {73.396617,30.990261,0};
> Point(68) = {73.040378,30.777995,0};
> Point(69) = {72.711744,30.525086,0};
> Point(70) = {72.415330,30.235083,0};
> Point(71) = {72.155295,29.912058,0};
> Point(72) = {71.935291,29.560546,0};
> Point(73) = {71.758405,29.185479,0};
> Point(74) = {71.627121,28.792125,0};
> Point(75) = {71.543282,28.386004,0};
> Point(76) = {71.508065,27.972818,0};
> Point(77) = {71.521964,27.558366,0};
> Point(78) = {71.584783,27.148468,0};
> Point(79) = {71.695642,26.748876,0};
> Point(80) = {71.852983,26.365201,0};
> Point(81) = {72.054599,26.002828,0};
> Point(82) = {72.297659,25.666844,0};
> Point(83) = {72.578750,25.361966,0};
> Point(84) = {72.893928,25.092473,0};
> Point(85) = {73.238766,24.862149,0};
> Point(86) = {73.608426,24.674226,0};
> Point(87) = {73.997718,24.531344,0};
> Point(88) = {74.401176,24.435507,0};
> Point(89) = {74.813137,24.388061,0};
> Point(90) = {75.227818,24.389672,0};
> Point(91) = {75.639398,24.440317,0};
> Point(92) = {76.042100,24.539286,0};
> Point(93) = {76.430269,24.685189,0};
> Point(94) = {76.798457,24.875978,0};
> Point(95) = {77.141496,25.108974,0};
> Point(96) = {77.454571,25.380908,0};
> Point(97) = {77.733285,25.687961,0};
> Point(98) = {77.973726,26.025823,0};
> Point(99) = {78.172520,26.389752,0};
> Point(100) = {78.326876,26.774638,0};
> Point(101) = {78.434627,27.175079,0};
> Point(102) = {78.494260,27.585453,0};
> Point(103) = {80.000000,28.000000,0};
> Point(104) = {78.504938,28.000000,0};
> Point(105) = {78.504938,28.000000,0};
> Point(106) = {76.000000,28.000000,0};
> Point(107) = {75.006904,27.882693,0};
> Line(1) = {1,2};
> Line(2) = {2,3};
> Line(3) = {3,4};
> Line(4) = {4,5};
> Line(5) = {5,6};
> Line(6) = {6,7};
> Line(7) = {7,8};
> Line(8) = {8,9};
> Line(9) = {9,10};
> Line(10) = {10,11};
> Line(11) = {11,12};
> Line(12) = {12,13};
> Line(13) = {13,14};
> Line(14) = {14,15};
> Line(15) = {15,16};
> Line(16) = {16,17};
> Line(17) = {17,18};
> Line(18) = {18,19};
> Line(19) = {19,20};
> Line(20) = {20,21};
> Line(21) = {21,22};
> Line(22) = {22,23};
> Line(23) = {23,24};
> Line(24) = {24,25};
> Line(25) = {25,26};
> Line(26) = {26,27};
> Line(27) = {27,28};
> Line(28) = {28,29};
> Line(29) = {29,30};
> Line(30) = {30,31};
> Line(31) = {31,32};
> Line(32) = {32,33};
> Line(33) = {33,34};
> Line(34) = {34,35};
> Line(35) = {35,36};
> Line(36) = {36,37};
> Line(37) = {37,38};
> Line(38) = {38,39};
> Line(39) = {39,40};
> Line(40) = {40,41};
> Line(41) = {41,42};
> Line(42) = {42,43};
> Line(43) = {43,44};
> Line(44) = {44,45};
> Line(45) = {45,46};
> Line(46) = {46,47};
> Line(47) = {47,1};
> Line(48) = {48,49};
> Line(49) = {49,50};
> Line(50) = {50,51};
> Line(51) = {51,52};
> Line(52) = {52,53};
> Line(53) = {53,54};
> Line(54) = {54,55};
> Line(55) = {55,56};
> Line(56) = {56,57};
> Line(57) = {57,58};
> Line(58) = {58,59};
> Line(59) = {59,60};
> Line(60) = {60,61};
> Line(61) = {61,62};
> Line(62) = {62,63};
> Line(63) = {63,64};
> Line(64) = {64,65};
> Line(65) = {65,66};
> Line(66) = {66,67};
> Line(67) = {67,68};
> Line(68) = {68,69};
> Line(69) = {69,70};
> Line(70) = {70,71};
> Line(71) = {71,72};
> Line(72) = {72,73};
> Line(73) = {73,74};
> Line(74) = {74,75};
> Line(75) = {75,76};
> Line(76) = {76,77};
> Line(77) = {77,78};
> Line(78) = {78,79};
> Line(79) = {79,80};
> Line(80) = {80,81};
> Line(81) = {81,82};
> Line(82) = {82,83};
> Line(83) = {83,84};
> Line(84) = {84,85};
> Line(85) = {85,86};
> Line(86) = {86,87};
> Line(87) = {87,88};
> Line(88) = {88,89};
> Line(89) = {89,90};
> Line(90) = {90,91};
> Line(91) = {91,92};
> Line(92) = {92,93};
> Line(93) = {93,94};
> Line(94) = {94,95};
> Line(95) = {95,96};
> Line(96) = {96,97};
> Line(97) = {97,98};
> Line(98) = {98,99};
> Line(99) = {99,100};
> Line(100) = {100,101};
> Line(101) = {101,102};
> Line(102) = {102,48};
> Line(103) = {2,48};
> Line(104) = {48,48};
> Line(105) = {48,106};
> Line(106) = {106,107};
> Line Loop(107) =
> {1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47};
>
> Plane Surface(108) = {107};
> Recombine Surface{108};
> Line{48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106}
> In Surface{108};
>
> Field[1]=MathEval;
> Field[1].F ="min(0.025+( 1-Exp( -((x-75.0069)^2+(y-27.8827)^2)/((1)^2) )
> )*(0.975))";
> Background Field = 1;
>
> Mesh.CharacteristicLengthExtendFromBoundary=0;
>
>
>
> On 09/28/2010 10:57 AM, Ruth V. Sabariego wrote:
>> Indeed, it seems there is a problem with Lines 52 and 80.
>> Changing the following option in your geo file, seems to work:
>>
>> Mesh.CharacteristicLengthExtendFromBoundary = 0;
>> // Extend computation of mesh element sizes from the boundaries into
>> the surfaces/volumes
>>
>> By default is set to 1.
>>
>> Regards,
>> Ruth
>>
>>
>> On 27/09/10 15:40, olivier jamond wrote:
>>> Hi, thank you for helping!
>>>
>>> I have just tried your solution. It behaves better but there is still
>>> some problematic cases... Here is one:
>>>
>>> Point(1) = {0.000000,0.000000,0};
>>> Point(2) = {80.000000,0.000000,0};
>>> Point(3) = {80.000000,28.000000,0};
>>> Point(4) = {80.000000,40.000000,0};
>>> Point(5) = {0.000000,40.000000,0};
>>> Point(6) = {0.000000,12.000000,0};
>>> Point(7) = {20.000000,28.000000,0};
>>> Point(8) = {18.730028,32.325127,0};
>>> Point(9) = {15.323320,35.277056,0};
>>> Point(10) = {10.861481,35.918572,0};
>>> Point(11) = {6.761114,34.045997,0};
>>> Point(12) = {4.324056,30.253860,0};
>>> Point(13) = {4.324056,25.746140,0};
>>> Point(14) = {6.761114,21.954003,0};
>>> Point(15) = {10.861481,20.081428,0};
>>> Point(16) = {15.323320,20.722944,0};
>>> Point(17) = {18.730028,23.674873,0};
>>> Point(18) = {76.000000,12.000000,0};
>>> Point(19) = {74.730028,16.325127,0};
>>> Point(20) = {71.323320,19.277056,0};
>>> Point(21) = {66.861481,19.918572,0};
>>> Point(22) = {62.761114,18.045997,0};
>>> Point(23) = {60.324056,14.253860,0};
>>> Point(24) = {60.324056,9.746140,0};
>>> Point(25) = {62.761114,5.954003,0};
>>> Point(26) = {66.861481,4.081428,0};
>>> Point(27) = {71.323320,4.722944,0};
>>> Point(28) = {74.730028,7.674873,0};
>>> Point(29) = {0.000000,12.000000,0};
>>> Point(30) = {4.000000,12.000000,0};
>>> Point(31) = {4.998921,12.046438,0};
>>> Point(32) = {5.976553,12.256763,0};
>>> Point(33) = {6.959056,12.443009,0};
>>> Point(34) = {7.917472,12.728383,0};
>>> Point(35) = {8.881660,12.993604,0};
>>> Point(36) = {9.847070,13.254340,0};
>>> Point(37) = {10.791857,13.582024,0};
>>> Point(38) = {11.732240,13.922143,0};
>>> Point(39) = {12.683953,14.229132,0};
>>> Point(40) = {13.637635,14.529946,0};
>>> Point(41) = {14.621307,14.709919,0};
>>> Point(42) = {15.612789,14.840166,0};
>>> Point(43) = {16.598456,15.008865,0};
>>> Point(44) = {17.595157,15.090027,0};
>>> Point(45) = {18.594970,15.109361,0};
>>> Point(46) = {19.594568,15.081007,0};
>>> Point(47) = {20.593809,15.042039,0};
>>> Point(48) = {21.587684,14.931535,0};
>>> Point(49) = {22.583660,14.841912,0};
>>> Point(50) = {23.579780,14.753904,0};
>>> Point(51) = {24.570180,14.615672,0};
>>> Point(52) = {25.556851,14.452945,0};
>>> Point(53) = {26.539498,14.267462,0};
>>> Point(54) = {27.538870,14.302900,0};
>>> Point(55) = {28.538787,14.315780,0};
>>> Point(56) = {29.538208,14.281758,0};
>>> Point(57) = {80.000000,28.000000,0};
>>> Point(58) = {76.000000,28.000000,0};
>>> Point(59) = {75.004838,27.901747,0};
>>> Point(60) = {74.022794,27.713099,0};
>>> Point(61) = {73.057275,27.452766,0};
>>> Point(62) = {72.100706,27.161260,0};
>>> Point(63) = {71.108506,27.036600,0};
>>> Point(64) = {70.143832,26.773154,0};
>>> Point(65) = {69.183306,26.494963,0};
>>> Point(66) = {68.253875,26.125967,0};
>>> Point(67) = {67.303473,25.814946,0};
>>> Point(68) = {66.338865,25.551256,0};
>>> Point(69) = {65.359379,25.349745,0};
>>> Point(70) = {64.389750,25.105165,0};
>>> Point(71) = {63.405932,24.925991,0};
>>> Point(72) = {62.419237,24.763411,0};
>>> Point(73) = {61.421902,24.690451,0};
>>> Point(74) = {60.422787,24.732517,0};
>>> Point(75) = {59.424031,24.782376,0};
>>> Point(76) = {58.441564,24.968812,0};
>>> Point(77) = {57.462291,25.171360,0};
>>> Point(78) = {56.476254,25.337884,0};
>>> Point(79) = {55.482089,25.445756,0};
>>> Point(80) = {54.482485,25.473869,0};
>>> Point(81) = {53.483407,25.430930,0};
>>> Point(82) = {52.483763,25.404241,0};
>>> Point(83) = {51.485758,25.467367,0};
>>> Point(84) = {50.487634,25.528594,0};
>>> Line(1) = {1,2};
>>> Line(2) = {2,3};
>>> Line(3) = {3,4};
>>> Line(4) = {4,5};
>>> Line(5) = {5,6};
>>> Line(6) = {6,1};
>>> Line(7) = {7,8};
>>> Line(8) = {8,9};
>>> Line(9) = {9,10};
>>> Line(10) = {10,11};
>>> Line(11) = {11,12};
>>> Line(12) = {12,13};
>>> Line(13) = {13,14};
>>> Line(14) = {14,15};
>>> Line(15) = {15,16};
>>> Line(16) = {16,17};
>>> Line(17) = {17,7};
>>> Line(18) = {18,19};
>>> Line(19) = {19,20};
>>> Line(20) = {20,21};
>>> Line(21) = {21,22};
>>> Line(22) = {22,23};
>>> Line(23) = {23,24};
>>> Line(24) = {24,25};
>>> Line(25) = {25,26};
>>> Line(26) = {26,27};
>>> Line(27) = {27,28};
>>> Line(28) = {28,18};
>>> Line(29) = {6,30};
>>> Line(30) = {30,31};
>>> Line(31) = {31,32};
>>> Line(32) = {32,33};
>>> Line(33) = {33,34};
>>> Line(34) = {34,35};
>>> Line(35) = {35,36};
>>> Line(36) = {36,37};
>>> Line(37) = {37,38};
>>> Line(38) = {38,39};
>>> Line(39) = {39,40};
>>> Line(40) = {40,41};
>>> Line(41) = {41,42};
>>> Line(42) = {42,43};
>>> Line(43) = {43,44};
>>> Line(44) = {44,45};
>>> Line(45) = {45,46};
>>> Line(46) = {46,47};
>>> Line(47) = {47,48};
>>> Line(48) = {48,49};
>>> Line(49) = {49,50};
>>> Line(50) = {50,51};
>>> Line(51) = {51,52};
>>> Line(52) = {52,53};
>>> Line(53) = {53,54};
>>> Line(54) = {54,55};
>>> Line(55) = {55,56};
>>> Line(56) = {3,58};
>>> Line(57) = {58,59};
>>> Line(58) = {59,60};
>>> Line(59) = {60,61};
>>> Line(60) = {61,62};
>>> Line(61) = {62,63};
>>> Line(62) = {63,64};
>>> Line(63) = {64,65};
>>> Line(64) = {65,66};
>>> Line(65) = {66,67};
>>> Line(66) = {67,68};
>>> Line(67) = {68,69};
>>> Line(68) = {69,70};
>>> Line(69) = {70,71};
>>> Line(70) = {71,72};
>>> Line(71) = {72,73};
>>> Line(72) = {73,74};
>>> Line(73) = {74,75};
>>> Line(74) = {75,76};
>>> Line(75) = {76,77};
>>> Line(76) = {77,78};
>>> Line(77) = {78,79};
>>> Line(78) = {79,80};
>>> Line(79) = {80,81};
>>> Line(80) = {81,82};
>>> Line(81) = {82,83};
>>> Line(82) = {83,84};
>>> Line Loop(83) = {1,2,3,4,5,6};
>>> Line Loop(84) = {7,8,9,10,11,12,13,14,15,16,17};
>>> Line Loop(85) = {18,19,20,21,22,23,24,25,26,27,28};
>>> Plane Surface(86) = {83,-84,-85};
>>> Line{29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82}
>>> In Surface{86};
>>>
>>> Field[1]=MathEval;
>>> Field[1].F ="min(0.25+( 1-Exp(
>>> -((x-29.5382)^2+(y-14.2818)^2)/((12.5)^2) ) )*(4.75),0.25+( 1-Exp(
>>> -((x-50.4876)^2+(y-25.5286)^2)/((12.5)^2) ) )*(4.75))";
>>> Background Field = 1;
>>>
>>> Thank you again for your help!
>>>
>>> Olivier
>>>
>>>> The problem is that Line(81) cuts the Surfaces of your problem.
>>>> You have to split in three. That is:
>>>> Line(81) = {40,80};
>>>> Line{81} In Surface{41};
>>>> Line(82) = {1,40};
>>>> Line(83) = {79,1};
>>>>
>>>> Regards,
>>>> Ruth
>>>>
>>>>
>>>> On 21/09/10 09:17, olivier jamond wrote:
>>>> >/ Hi,
>>>> />/
>>>> />/ I am a PhD student and I am working with evolving cracks in
>>>> />/ finite-elements modelization. I am dealing with a finite-elements
>>>> />/ patch located in the vicinity of the crack tip. Its mesh has to be
>>>> />/ regenerated at each step of the propagation in order to fit the
>>>> crack.
>>>> />/ So I generate automatically a .geo files and call gmsh at each
>>>> step of
>>>> />/ the propagation. Some time ago, I asked in this mailling list
>>>> how to
>>>> />/ introduce cracks with gmsh. I had been advised to use the command
>>>> />/ "Line{x} In Surface{y}".
>>>> />/
>>>> />/ But in some cases, some elements do not respect this constraint and
>>>> />/ cross the line (with the 3 2D mesh algorithms), which is really
>>>> />/ problematic in my developments. Here is such a case:
>>>> />/
>>>> />/ Point(1) = {2.200000,5.000000,0};
>>>> />/ Point(2) = {2.223310,4.711260,0};
>>>> />/ Point(3) = {2.292634,4.429998,0};
>>>> />/ Point(4) = {2.406179,4.163498,0};
>>>> />/ Point(5) = {2.561003,3.918664,0};
>>>> />/ Point(6) = {2.753096,3.701836,0};
>>>> />/ Point(7) = {2.977483,3.518629,0};
>>>> />/ Point(8) = {3.228353,3.373789,0};
>>>> />/ Point(9) = {3.499209,3.271067,0};
>>>> />/ Point(10) = {3.783034,3.213124,0};
>>>> />/ Point(11) = {4.072479,3.201460,0};
>>>> />/ Point(12) = {4.360046,3.236377,0};
>>>> />/ Point(13) = {4.638289,3.316971,0};
>>>> />/ Point(14) = {4.900000,3.441154,0};
>>>> />/ Point(15) = {5.138402,3.605711,0};
>>>> />/ Point(16) = {5.347319,3.806379,0};
>>>> />/ Point(17) = {5.521342,4.037962,0};
>>>> />/ Point(18) = {5.655963,4.294460,0};
>>>> />/ Point(19) = {5.747695,4.569232,0};
>>>> />/ Point(20) = {5.794163,4.855160,0};
>>>> />/ Point(21) = {5.794163,5.144840,0};
>>>> />/ Point(22) = {5.747695,5.430768,0};
>>>> />/ Point(23) = {5.655963,5.705540,0};
>>>> />/ Point(24) = {5.521342,5.962038,0};
>>>> />/ Point(25) = {5.347319,6.193621,0};
>>>> />/ Point(26) = {5.138402,6.394289,0};
>>>> />/ Point(27) = {4.900000,6.558846,0};
>>>> />/ Point(28) = {4.638289,6.683029,0};
>>>> />/ Point(29) = {4.360046,6.763623,0};
>>>> />/ Point(30) = {4.072479,6.798540,0};
>>>> />/ Point(31) = {3.783034,6.786876,0};
>>>> />/ Point(32) = {3.499209,6.728933,0};
>>>> />/ Point(33) = {3.228353,6.626211,0};
>>>> />/ Point(34) = {2.977483,6.481371,0};
>>>> />/ Point(35) = {2.753096,6.298164,0};
>>>> />/ Point(36) = {2.561003,6.081336,0};
>>>> />/ Point(37) = {2.406179,5.836502,0};
>>>> />/ Point(38) = {2.292634,5.570002,0};
>>>> />/ Point(39) = {2.223310,5.288740,0};
>>>> />/ Line(1) = {1,2};
>>>> />/ Line(2) = {2,3};
>>>> />/ Line(3) = {3,4};
>>>> />/ Line(4) = {4,5};
>>>> />/ Line(5) = {5,6};
>>>> />/ Line(6) = {6,7};
>>>> />/ Line(7) = {7,8};
>>>> />/ Line(8) = {8,9};
>>>> />/ Line(9) = {9,10};
>>>> />/ Line(10) = {10,11};
>>>> />/ Line(11) = {11,12};
>>>> />/ Line(12) = {12,13};
>>>> />/ Line(13) = {13,14};
>>>> />/ Line(14) = {14,15};
>>>> />/ Line(15) = {15,16};
>>>> />/ Line(16) = {16,17};
>>>> />/ Line(17) = {17,18};
>>>> />/ Line(18) = {18,19};
>>>> />/ Line(19) = {19,20};
>>>> />/ Line(20) = {20,21};
>>>> />/ Line(21) = {21,22};
>>>> />/ Line(22) = {22,23};
>>>> />/ Line(23) = {23,24};
>>>> />/ Line(24) = {24,25};
>>>> />/ Line(25) = {25,26};
>>>> />/ Line(26) = {26,27};
>>>> />/ Line(27) = {27,28};
>>>> />/ Line(28) = {28,29};
>>>> />/ Line(29) = {29,30};
>>>> />/ Line(30) = {30,31};
>>>> />/ Line(31) = {31,32};
>>>> />/ Line(32) = {32,33};
>>>> />/ Line(33) = {33,34};
>>>> />/ Line(34) = {34,35};
>>>> />/ Line(35) = {35,36};
>>>> />/ Line(36) = {36,37};
>>>> />/ Line(37) = {37,38};
>>>> />/ Line(38) = {38,39};
>>>> />/ Line(39) = {39,1};
>>>> />/ Line Loop(40) =
>>>> />/
>>>> {1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39};
>>>>
>>>> />/ Plane Surface(41) = {40};
>>>> />/ //Recombine Surface{41};
>>>> />/ Point(40) = {2.600000,5.000000,0};
>>>> />/ Point(41) = {2.618130,4.775424,0};
>>>> />/ Point(42) = {2.672049,4.556665,0};
>>>> />/ Point(43) = {2.760362,4.349388,0};
>>>> />/ Point(44) = {2.880780,4.158961,0};
>>>> />/ Point(45) = {3.030186,3.990317,0};
>>>> />/ Point(46) = {3.204709,3.847823,0};
>>>> />/ Point(47) = {3.399830,3.735169,0};
>>>> />/ Point(48) = {3.610496,3.655275,0};
>>>> />/ Point(49) = {3.831249,3.610208,0};
>>>> />/ Point(50) = {4.056372,3.601135,0};
>>>> />/ Point(51) = {4.280036,3.628293,0};
>>>> />/ Point(52) = {4.496447,3.690977,0};
>>>> />/ Point(53) = {4.700000,3.787564,0};
>>>> />/ Point(54) = {4.885424,3.915553,0};
>>>> />/ Point(55) = {5.047915,4.071628,0};
>>>> />/ Point(56) = {5.183266,4.251748,0};
>>>> />/ Point(57) = {5.287971,4.451247,0};
>>>> />/ Point(58) = {5.359319,4.664958,0};
>>>> />/ Point(59) = {5.395460,4.887347,0};
>>>> />/ Point(60) = {5.395460,5.112653,0};
>>>> />/ Point(61) = {5.359319,5.335042,0};
>>>> />/ Point(62) = {5.287971,5.548753,0};
>>>> />/ Point(63) = {5.183266,5.748252,0};
>>>> />/ Point(64) = {5.047915,5.928372,0};
>>>> />/ Point(65) = {4.885424,6.084447,0};
>>>> />/ Point(66) = {4.700000,6.212436,0};
>>>> />/ Point(67) = {4.496447,6.309023,0};
>>>> />/ Point(68) = {4.280036,6.371707,0};
>>>> />/ Point(69) = {4.056372,6.398865,0};
>>>> />/ Point(70) = {3.831249,6.389792,0};
>>>> />/ Point(71) = {3.610496,6.344725,0};
>>>> />/ Point(72) = {3.399830,6.264831,0};
>>>> />/ Point(73) = {3.204709,6.152177,0};
>>>> />/ Point(74) = {3.030186,6.009683,0};
>>>> />/ Point(75) = {2.880780,5.841039,0};
>>>> />/ Point(76) = {2.760362,5.650612,0};
>>>> />/ Point(77) = {2.672049,5.443335,0};
>>>> />/ Point(78) = {2.618130,5.224576,0};
>>>> />/ Line(42) = {40,41};
>>>> />/ Line(43) = {41,42};
>>>> />/ Line(44) = {42,43};
>>>> />/ Line(45) = {43,44};
>>>> />/ Line(46) = {44,45};
>>>> />/ Line(47) = {45,46};
>>>> />/ Line(48) = {46,47};
>>>> />/ Line(49) = {47,48};
>>>> />/ Line(50) = {48,49};
>>>> />/ Line(51) = {49,50};
>>>> />/ Line(52) = {50,51};
>>>> />/ Line(53) = {51,52};
>>>> />/ Line(54) = {52,53};
>>>> />/ Line(55) = {53,54};
>>>> />/ Line(56) = {54,55};
>>>> />/ Line(57) = {55,56};
>>>> />/ Line(58) = {56,57};
>>>> />/ Line(59) = {57,58};
>>>> />/ Line(60) = {58,59};
>>>> />/ Line(61) = {59,60};
>>>> />/ Line(62) = {60,61};
>>>> />/ Line(63) = {61,62};
>>>> />/ Line(64) = {62,63};
>>>> />/ Line(65) = {63,64};
>>>> />/ Line(66) = {64,65};
>>>> />/ Line(67) = {65,66};
>>>> />/ Line(68) = {66,67};
>>>> />/ Line(69) = {67,68};
>>>> />/ Line(70) = {68,69};
>>>> />/ Line(71) = {69,70};
>>>> />/ Line(72) = {70,71};
>>>> />/ Line(73) = {71,72};
>>>> />/ Line(74) = {72,73};
>>>> />/ Line(75) = {73,74};
>>>> />/ Line(76) = {74,75};
>>>> />/ Line(77) = {75,76};
>>>> />/ Line(78) = {76,77};
>>>> />/ Line(79) = {77,78};
>>>> />/ Line(80) = {78,40};
>>>> />/
>>>> Line{42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80}
>>>>
>>>> />/ In Surface{41};
>>>> />/ Point(79) = {-1.000000,5.000000,0};
>>>> />/ Point(80) = {4.000000,5.000000,0};
>>>> />/ Line(81) = {79,80};
>>>> />/ Line{81} In Surface{41};
>>>> />/
>>>> />/ Field[1]=MathEval;
>>>> />/ Field[1].F ="0.080000+( 1-Exp(
>>>> />/ -((x-4.000000)^2+(y-5.000000)^2)/((1.600000)^2) )
>>>> )*(0.400000-0.080000)";
>>>> />/ Background Field = 1;
>>>> />/
>>>> />/ Is there something wrong with my input file or any kind of bug
>>>> (I am
>>>> />/ using gmsh 2.4.2) ?
>>>> />/
>>>> />/ Thanks,
>>>> />/
>>>> />/ Best regards,
>>>> />/
>>>> />/ Olivier Jamond
>>>> />/
>>>> />/ _______________________________________________
>>>> />/ gmsh mailing list
>>>> />/ gmsh at geuz.org <http://www.geuz.org/mailman/listinfo/gmsh>
>>>> />/ http://www.geuz.org/mailman/listinfo/gmsh
>>>> />/
>>>> /
>>>> --
>>>> Dr. Ir. Ruth V. Sabariego
>>>> University of Liege, Dept. of Electrical Engineering& Computer Science,
>>>> Applied& Computational Electromagnetics (ACE),
>>>> phone: +32-4-3663737 - fax: +32-4-3662910
>>>> -http://ace.montefiore.ulg.ac.be/
>>>
>>>
>>> _______________________________________________
>>> gmsh mailing list
>>> gmsh at geuz.org
>>> http://www.geuz.org/mailman/listinfo/gmsh
>>>
>>
>
>
> _______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh
>

-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine