[Gmsh] Segmentation fault on extrude

Ruth V. Sabariego r.sabariego at ulg.ac.be
Tue Nov 15 22:14:33 CET 2011


You should extrude with an angle strictly smaller than Pi. 
Otherwise, there is a problem with indetermination.
Try:

// foil parameters
r1 = 2.0;
d  = 0.1;
r2 = r1 + d;

// size parameters
h1 =  4.0;
h2 = 10.0;
b = 2*r2*1.1;

// Mesh
lc = 2*d;
lv = 5;

p1 = newp; Point(p1) = {r1  , 0, h1+d, lc};

wall[] = Extrude  {0, 0, h2-h1-2*d} { Point{ p1 }; Layers{ lv }; };

wall1[] = Extrude  { {0, 0, 1}, {0, 0, 0},Pi/2 } { Line { wall[1] }; };
wall2[] = Extrude  { {0, 0, 1}, {0, 0, 0}, Pi/2 } { Line { wall1[0] }; };

wall3[] = Extrude  { {0, 0, 1}, {0, 0, 0},-Pi/2 } { Line { wall[1] }; };
wall4[] = Extrude  { {0, 0, 1}, {0, 0, 0},-Pi/2 } { Line { wall3[0] }; };


Regards,
Ruth



--
Dr. Ir. Ruth V. Sabariego
University of Liege, Electrical Engineering & Computer Science, 
Applied & Computational Electromagnetics (ACE),
phone: +32-4-3663737 - fax: +32-4-3662910 - http://ace.montefiore.ulg.ac.be/



On 15 Nov 2011, at 20:27, Anton Gladky wrote:

> Dear developers,
> 
> there is a bug on Launchpad bug-tracker [1].
> The attached script crashes. Seems, the problem is in last lines:
> 
> //Fails:
> wall2[] = Extrude  { {0, 0, 1}, {0, 0, 0},-Pi } { Line { wall[1] }; };
> //Fails:
> wall3[] = Extrude  { {0, 0, 1}, {0, 0, 0}, Pi } { Line { wall1[0] }; };
> 
> Could you, please, have a look at this script?
> 
> Thanks
> Anton
> 
> [1] https://bugs.launchpad.net/ubuntu/+source/gmsh/+bug/886607
> <atm2.geo>_______________________________________________
> gmsh mailing list
> gmsh at geuz.org
> http://www.geuz.org/mailman/listinfo/gmsh