Problem with Cut method and output resulting map
Hi Christophe,
I'm using gmsh 4.7.1 and I have a little problem. I have to create a parametric hole in a surface and then read the resulting edges. So, my workaround is: create the surface(2, tag1), create the hole(2, tag2), cut the surface with the hole. And so far so good. The problem arises when the hole divides the surface into 2 surfaces. Using gmsh.model.occ.cut everything works but the outMap objectDimTags and toolDimTags return only (2, tag1) and (2, tag2). And the second half of the surface is lost. Is there any way to fix it?
An example is easy to create.
// SetFactory("OpenCASCADE"); cl__1 = 1e+22;
Point(1) = {0, 0, 0, cl__1}; Point(2) = {100, 0, 0, cl__1}; Point(3) = {100, 100, 0, cl__1}; Point(4) = {0, 100, 0, cl__1}; Line(1) = {1, 2}; Line(2) = {2, 3}; Line(3) = {3, 4}; Line(4) = {4, 1};
Line Loop(1) = {1, 2, 3, 4}; Plane Surface(1) = {1};
Point(5) = {0, 30, 0, cl__1}; Point(6) = {100, 30, 0, cl__1}; Point(7) = {100, 50, 0, cl__1}; Point(8) = {0, 50, 0, cl__1}; Line(5) = {5, 6}; Line(6) = {6, 7}; Line(7) = {7, 8}; Line(8) = {8, 5};
Line Loop(2) = {5, 6, 7, 8}; Plane Surface(2) = {2};
BooleanDifference{ Surface{1}; Delete; }{ Surface{2}; Delete; } //
Thanks a lot
Gabriele