[Gmsh] Meshing challenge: extrusion of compound surfaces

[François Beaubert]

Hi all,

I'm discovering gmsh and it's new "Compound" feature.
I'm wondering if it's possible to extrude a compound of surfaces using a 
revolution extrusion.

I've tried many ways but I'm always getting this error message:

"Error   : Impossible to extrude entity 1000 (of type 307)"

Is it a bug a limitation or do I something wrong ?

Please find enclosed the small example I use to make my test:

Any help will be welcome
Thanks !

Francois

---------------------- geo file starts here 
-------------------------------------

Point(1) = {0, 0, 0};
Point(2) = {10, 0, 0};
Point(3) = {10, 0, 10};
Point(4) = {0, 0, 10};
Point(5) = {10, 0, 20};
Point(6) = {0, 0, 20};
Point(7) = {10, 0, 30};
Point(8) = {0, 0, 30};

Line(1) = {1, 4};
Line(2) = {2, 3};
Line(3) = {1, 2};
Line(4) = {4, 3};
Line(5) = {6, 5};
Line(6) = {3, 5};
Line(7) = {4, 6};
Line(8) = {8, 7};
Line(9) = {5, 7};
Line(10)= {6, 8};

Line Loop(1) = {3, 2, -4, -1};
Plane Surface(1) = {1};

Line Loop(2) = {4, 6, -5, -7};
Plane Surface(2) = {2};

Line Loop(3) = {5, 9, -8, -10};
Plane Surface(3) = {3};

Transfinite Line {1, 2, 3, 4, 5, 6, 7, 8, 9, 10} = 3 Using Progression 1;

Transfinite Surface{1} = {1,2,3,4};
Recombine Surface{1};

Transfinite Surface{2} = {4,3,5,6};
Recombine Surface{2};

Transfinite Surface{3}={6,5,7,8};
Recombine Surface{3};

Compound Surface(1000) = {1,2,3} Boundary {{},{},{},{}};

// revolution[] = Extrude { {0,0,1} , {0,0,0} , Pi/2 } {
//   Surface{1,2,3}; Layers{{2,1,2}, {0.25,0.75,1}}; Recombine;
// };

revolution[] = Extrude { {0,0,1} , {0,0,0} , Pi/2 } {
   Surface{1000}; Layers{{2,1,2}, {0.25,0.75,1}}; Recombine;
};

---------------------- geo file ends here 
-------------------------------------

-------------- next part --------------
A non-text attachment was scrubbed...
Name: francois_beaubert.vcf
Type: text/x-vcard
Size: 390 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/gmsh/attachments/20110120/90a071d9/attachment.vcf>