// gmsh 2011-03-29  -*- c++ -*-

r = 1.0;  // cylinder radius
R = 10.0; // outer boundary radius

n1 = 10;
n2 = 10;

Point(1) = { 0, 0, 0};
Point(2) = {0, r, 0};
Point(3) = { r, 0, 0};

Point(4) = {0, R, 0};
Point(5) = { R, 0, 0};

Circle(1) = {3,1,2};
Circle(2) = {5,1,4};

Line(3) = {5,3};
Line(4) = {2,4};

Line Loop(1) = {3, 1, 4, -2};
Ruled Surface(1) = {1};
Transfinite Surface(1) = {3,5,4,2};
Recombine Surface(1);

Transfinite Line {3,4} = n1;
Transfinite Line {1,2} = n2;

Symmetry {0.0, 1.0, 0.0, 0.0} { Duplicata{Surface{1};} }

Symmetry {1, 0, 0, 0} {
  Duplicata { Surface{1, 5}; }
}
