Occ cut fails in python API when using addPlaneSurface
Hi, I'm running into an issue when trying to perform a cut between two faces created with gmsh.model.occ.addPlaneSurface. See the following script to reproduce the pb NB using Shape from Occ like Rectangle works fine but I really need to define the face with point and lines
import gmsh
gmsh.initialize()
gmsh.option.setNumber("General.Terminal", 1)
gmsh.option.setNumber("General.Verbosity", 0)
# # working
# id0 = gmsh.model.occ.addRectangle(-2, -2, 0, 4, 4)
# id1 = gmsh.model.occ.addRectangle(-1, -1, 0, 2, 2)
# cad = gmsh.model.occ.cut([(2, id0)], [(2, id1)])
pt0_id0 = gmsh.model.occ.addPoint(-2, -2, 0)
pt1_id0 = gmsh.model.occ.addPoint( 2, -2, 0)
pt2_id0 = gmsh.model.occ.addPoint( 2, 2, 0)
pt3_id0 = gmsh.model.occ.addPoint(-2, 2, 0)
l0 = gmsh.model.occ.addLine(pt0_id0, pt1_id0)
l1 = gmsh.model.occ.addLine(pt1_id0, pt2_id0)
l2 = gmsh.model.occ.addLine(pt2_id0, pt3_id0)
l3 = gmsh.model.occ.addLine(pt3_id0, pt0_id0)
cl = gmsh.model.occ.addCurveLoop([l0, l1, l2, l3])
disk = gmsh.model.occ.addPlaneSurface([cl])
pt0_id0 = gmsh.model.occ.addPoint(-1, -1, 0)
pt1_id0 = gmsh.model.occ.addPoint( 1, -1, 0)
pt2_id0 = gmsh.model.occ.addPoint( 1, 1, 0)
pt3_id0 = gmsh.model.occ.addPoint(-1, 1, 0)
l0 = gmsh.model.occ.addLine(pt0_id0, pt1_id0)
l1 = gmsh.model.occ.addLine(pt1_id0, pt2_id0)
l2 = gmsh.model.occ.addLine(pt2_id0, pt3_id0)
l3 = gmsh.model.occ.addLine(pt3_id0, pt0_id0)
disk1 = gmsh.model.occ.addPlaneSurface([cl])
cad = gmsh.model.occ.cut([(2, disk)], [(2, disk1)])
gmsh.model.occ.synchronize()
gmsh.fltk.run()
gmsh.finalize()
Edited by Christophe Trophime