hi to everyone,

I would like to ask please;

I created a square with a hole using the commands “square” and then “trunc”.

I wrote the following code in FF++ for the square domain with a hole :

//parameters

real x0=-7.,y0=-7.;//left bottom node of the omega square

real xF=7.,yF=7.;//right top node of the omega square

int n=7;

int N =2^n;//number of partition in each direction

real xT=-2.1,yT=-2.1;//left bottom node of TRAC square

real xTF=2.1,yTF=2.1;//right top node of TRAC square

// domain of the problem

int[int] outer=[1,1,1,1], bTRAC=[2,2,2,2];

mesh Th = square(N,N,[x0+(xF-x0)*x ,y0+(yF-y0)*y],label=outer, region=1);

mesh ThTRAC = square(N,N,[xT+(xTF-xT)*x,yT+(yTF-yT)*y],label=bTRAC,region=2);

mesh Thwithole=trunc(Th,((x<xT | x>xTF|y<yT) | (x<xT | x>xTF| y>yTF)));

I know that if add the following lines:

int outerboundary=1;

mesh Thwithole=trunc(Th,((x<xT | x>xTF|y<yT) | (x<xT | x>xTF| y>yTF)),label=outerboundary);

Then FF++ would assign the same label for each of the boundaries.

But, I would like to assign different labels for each of the boundaries, i.e., for the inner boundary and for the outer boundary.

Thanks in advance,

Mordechai.