/* --------------------------------------------------------------------------
    This is a gmsh geometry file
   --------------------------------------------------------------------------*/

/* defining some parameters */

// resolution
res = 2e-2;

/* Width and height of domain */
L = 1; 
/* Radius of circular scatterer */
r = 0.1;

/* defining the points */
p1 = newp; Point(p1) = {-L/2,-L/2,0,res};
p2 = newp; Point(p2) = {L/2,-L/2,0,res};
p3 = newp; Point(p3) = {-L/2,L/2,0,res};
p4 = newp; Point(p4) = {L/2,L/2,0,res};

p5 = newp; Point(p5) = {0,0,0,res};
p6 = newp; Point(p6) = {r,0,0,res};
p7 = newp; Point(p7) = {-r,0,0,res};

/* defining the lines */
l1 = newl; Line(l1) = {p1,p2};
l2 = newl; Line(l2) = {p2,p4};
l3 = newl; Line(l3) = {p4,p3};
l4 = newl; Line(l4) = {p3,p1};

l5 = newl; Circle(l5) = {p6,p5,p7};
l6 = newl; Circle(l6) = {p7,p5,p6};

loop1 = newll; Line Loop(loop1) = {l1,l2,l3,l4,-l5,-l6}; 

s1 = news; Plane Surface(s1) = {loop1};

Physical Surface(1001) = {s1};
Physical Line(101) = {l5,l6};



