New section test (#159)

mpoudot · web-flow · commit f92e745157a7 · 2025-08-08T16:35:20.000+02:00
* Remove inoperative code for section of surfaces

* New section test

* test_sectionbyplane enhanced
diff --git a/src/Core/Geom/CommandSection.cpp b/src/Core/Geom/CommandSection.cpp
@@ -44,14 +44,13 @@ CommandSection::~CommandSection()
 void CommandSection::validate()
 {
     for(unsigned int i=0;i<m_entities.size();i++){
-        if(m_entities[i]->getDim()<2)
+        if(m_entities[i]->getDim()<=2)
         throw TkUtil::Exception(TkUtil::UTF8String ("Seuls des volumes géométriques peuvent être découpés", TkUtil::Charset::UTF_8));
     }
 
     int dim_tool = m_tool->getDim();
     if(dim_tool!=2)
         throw TkUtil::Exception(TkUtil::UTF8String ("L'outil de découpe doit être une surface", TkUtil::Charset::UTF_8));
-
 }
 /*----------------------------------------------------------------------------*/
 void CommandSection::internalSpecificPreExecute()
diff --git a/src/Core/Geom/GeomSectionImplementation.cpp b/src/Core/Geom/GeomSectionImplementation.cpp
@@ -12,6 +12,7 @@
 #include <BRepAlgoAPI_BuilderAlgo.hxx>
 #include <BRepAlgoAPI_Fuse.hxx>
 #include <BRepAlgoAPI_Common.hxx>
+#include <BRepAlgoAPI_Cut.hxx>
 #include <BRepFeat_Gluer.hxx>
 #include <TopTools_ListOfShape.hxx>
 /*----------------------------------------------------------------------------*/
@@ -67,10 +68,8 @@ void GeomSectionImplementation::perform(std::vector<GeomEntity*>& res)
     if(dim_ref==3){
     	sectionVolumes(res);
     }
-    else if(dim_ref==2)
-        sectionSurfaces(res);
     else
-        throw TkUtil::Exception (TkUtil::UTF8String ("La section géométrique ne s'applique qu'aux volumes ou surfaces", TkUtil::Charset::UTF_8));
+        throw TkUtil::Exception (TkUtil::UTF8String ("La section géométrique ne s'applique qu'aux volumes", TkUtil::Charset::UTF_8));
 
 }
 /*----------------------------------------------------------------------------*/
@@ -124,61 +123,6 @@ void GeomSectionImplementation::sectionVolumes(std::vector<GeomEntity*>& res)
     res.insert(res.end(), entities_new.begin(), entities_new.end());
 }
 /*----------------------------------------------------------------------------*/
-void GeomSectionImplementation::sectionSurfaces(std::vector<GeomEntity*>& res)
-{
-    //========================================================================
-    //1 - Recuperation de la surface de coupe
-    //========================================================================
-    TopoDS_Face tool_face = checkSurface(m_tool);
-
-    //========================================================================
-    // 2 - Conservation du plan de coupe réduit aux objets coupés
-    //========================================================================
-    // On commence par fusionner toutes les entites a couper
-
-    GeomEntity* e1 = m_init_entities[0];
-//    TopoDS_Shape s_fuse;
-//    getOCCShape(e1, s_fuse);
-//    for(unsigned int i=1;i<m_init_entities.size();i++){
-//        GeomEntity* e2 = m_init_entities[i];
-//        TopoDS_Shape s2;
-//        getOCCShape(e2, s2);
-//
-//        BRepAlgoAPI_Fuse fuse_operator(s_fuse,s2);
-//        if(fuse_operator.IsDone())
-//            s_fuse = fuse_operator.Shape();
-//        else
-//            throw TkUtil::Exception (TkUtil::UTF8String ("Problème OCC lors de l'union avant coupe", TkUtil::Charset::UTF_8);
-//    }
-//    // On recupere l'intersection de la surface wf et de l'union des entites à
-//    // couper
-//    BRepAlgoAPI_Common common_operator(s_fuse,tool_face);
-//    if(common_operator.IsDone())
-//        m_restricted_section_tool = common_operator.Shape();
-//    else
-//        throw TkUtil::Exception (TkUtil::UTF8String ("Problème OCC lors de l'intersection avant coupe", TkUtil::Charset::UTF_8);
-
-    //========================================================================
-    //3 - Decoupe des entités
-    //========================================================================
-    TopoDS_Shape s;
-    BRepAlgoAPI_BuilderAlgo splitter;
-    TopTools_ListOfShape list_of_arguments;
-    for (unsigned int i = 0; i < m_init_entities.size(); i++) {
-        TopoDS_Shape si = checkSurface(m_init_entities[i]);
-        list_of_arguments.Append(si);
-    }
-    list_of_arguments.Append(tool_face);
-    splitter.Build();
-    s = splitter.Shape();
-
-
-    //createGeomEntitiesBut(s,tool_shape);
-    createGeomEntities(s,true);
-
-    res.insert(res.end(), entities_new.begin(), entities_new.end());
-}
-/*----------------------------------------------------------------------------*/
 TopoDS_Face GeomSectionImplementation::
 checkSurface(GeomEntity* ge) const
 {
diff --git a/src/Core/protected/Geom/GeomSectionImplementation.h b/src/Core/protected/Geom/GeomSectionImplementation.h
@@ -68,11 +68,6 @@ class GeomSectionImplementation: public GeomModificationBaseClass{
      */
     void sectionVolumes(std::vector<GeomEntity*>& res);
 
-    /*------------------------------------------------------------------------*/
-    /** \brief  réalisation de la coupe par m_tool
-     */
-    void sectionSurfaces(std::vector<GeomEntity*>& res);
-
 private:
     TopoDS_Face checkSurface(GeomEntity* ge) const;
     TopoDS_Shape checkVolume(GeomEntity* ge) const;
diff --git a/test_link/test_align_vertices.py b/test_link/test_align_vertices.py
@@ -17,6 +17,9 @@ def test_align_vertices(capfd):
     tm.alignVertices(Mgx3D.Point(0, 0, 1), Mgx3D.Point(0, 0, 0), ["Som0008"])
     assertPoint(tm, 0, 0, 0.5)
 
+    out, err = capfd.readouterr()
+    assert len(err) == 0
+
 def assertPoint(tm, x, y, z):
     p = tm.getCoord("Som0008")
     assert p.getX() == x
diff --git a/test_link/test_section.py b/test_link/test_section.py
@@ -0,0 +1,103 @@
+import os
+import pyMagix3D as Mgx3D
+
+def test_section_box1(capfd):
+    ctx = Mgx3D.getStdContext()
+    gm = ctx.getGeomManager()
+    tm = ctx.getTopoManager()
+    mm = ctx.getMeshManager()
+    ctx.clearSession() # Clean the session after the previous test
+
+    # Création de la boite Vol0000
+    gm.newBox (Mgx3D.Point(0, 0, 0), Mgx3D.Point(1, 1, 1), "A")
+    # Création de la boite Vol0001
+    gm.newBox (Mgx3D.Point(1, 0, 0), Mgx3D.Point(2, 1, 1), "B")
+    # Collage entre Vol0000 Vol0001
+    gm.glue (["Vol0000","Vol0001"])
+    # Création d'un bloc topologique structuré sur une géométrie (Vol0000)
+    ctx.getTopoManager().newStructuredTopoOnGeometry ("Vol0000")
+    # Création d'un bloc topologique structuré sur une géométrie (Vol0001)
+    ctx.getTopoManager().newStructuredTopoOnGeometry ("Vol0001")
+
+    # Création du sommet Pt0016 sur la courbe Crb0001
+    gm.newVertex ("Crb0001", 5.000000e-01, "C")
+    # Création du sommet Pt0017 sur la courbe Crb0003
+    gm.newVertex ("Crb0003", 5.000000e-01, "C")
+    # Création du sommet Pt0018 sur la courbe Crb0007
+    gm.newVertex ("Crb0007", 5.000000e-01, "C")
+    # Création du sommet Pt0019 sur la courbe Crb0005
+    gm.newVertex ("Crb0005", 5.000000e-01, "C")
+
+    # Création du segment Crb0024
+    gm.newSegment("Pt0016", "Pt0017", "C")
+    # Création du segment Crb0025
+    gm.newSegment("Pt0018", "Pt0017", "C")
+    # Création du segment Crb0026
+    gm.newSegment("Pt0019", "Pt0018", "C")
+    # Création du segment Crb0027
+    gm.newSegment("Pt0016", "Pt0019", "C")
+    # Création de la surface Surf0012
+    gm.newPlanarSurface (["Crb0024","Crb0025","Crb0026","Crb0027"], "C")
+
+    # Découpage suivant Ar0008 des blocs Bl0000
+    tm.splitBlocks (["Bl0000"],"Ar0008", .5)
+    # Section par un plan entre géométries avec topologies
+    gm.section (["Vol0000"], "Surf0012")
+    # Création du maillage pour tous les blocs
+    mm.newAllBlocksMesh()
+    # Sauvegarde du maillage (mli)
+    filename = "section.mli2"
+    mm.writeMli(filename)
+
+    assert os.path.exists(filename)
+    assert os.path.getsize(filename) > 0
+    assert gm.getNbVolumes() == 3
+    assert gm.getNbSurfaces() == 17
+    assert tm.getNbBlocks() == 3
+    assert tm.getNbFaces() == 16
+    assert mm.getNbFaces() == 1200
+
+    out, err = capfd.readouterr()
+    assert len(err) == 0
+
+
+def test_section_box2(capfd):
+    ctx = Mgx3D.getStdContext()
+    gm = ctx.getGeomManager()
+    ctx.clearSession() # Clean the session after the previous test
+
+    # Création de la boite Vol0000
+    gm.newBox (Mgx3D.Point(0, 0, 0), Mgx3D.Point(1, 1, 1), "A")
+    # Création de la boite Vol0001
+    gm.newBox (Mgx3D.Point(1, 0, 0), Mgx3D.Point(2, 1, 1), "B")
+    # Collage entre Vol0000 Vol0001
+    gm.glue (["Vol0000","Vol0001"])
+    # Création de la boite Vol0002
+    gm.newBox (Mgx3D.Point(0, .9, 0), Mgx3D.Point(2, 1.3, 1), "C")
+    # Section entre géométries
+    gm.section (["Vol0000","Vol0001"], "Surf0014")
+
+    assert gm.getNbVolumes() == 5
+    assert gm.getNbSurfaces() == 26
+
+    out, err = capfd.readouterr()
+    assert len(err) == 0
+
+
+def test_section_box_by_cylinder(capfd):
+    ctx = Mgx3D.getStdContext()
+    gm = ctx.getGeomManager()
+    ctx.clearSession() # Clean the session after the previous test
+
+    # Création de la boite Vol0000
+    gm.newBox (Mgx3D.Point(-.1, -.25, -.1), Mgx3D.Point(.1, .25, .1), "A")
+    # Création du cylindre Vol0001
+    gm.newCylinder (Mgx3D.Point(-.25, 0, 0), 0.05, Mgx3D.Vector(.5, 0, 0), 3.600000e+02, "B")
+    # Section suivant Surf0006 des entités Vol0000
+    gm.section (["Vol0000"], "Surf0006")
+
+    assert gm.getNbVolumes() == 3
+    assert gm.getNbSurfaces() == 12
+
+    out, err = capfd.readouterr()
+    assert len(err) == 0
diff --git a/test_link/test_sectionbyplane.py b/test_link/test_sectionbyplane.py
@@ -4,6 +4,21 @@
 # pour laquelle un certain nombres de bug a été trouvé
 # (cf. Issue#208)
 
+def test_cylinder(capfd):
+    ctx = Mgx3D.getStdContext()
+    ctx.clearSession() # Clean the session after the previous test
+    gm = ctx.getGeomManager()
+
+    # Création du cylindre Vol0000
+    gm.newCylinder (Mgx3D.Point(0, 0, 0), .1, Mgx3D.Vector(0, 0, 1), 180, "CYL")
+    # Section par un plan de Vol0000 suivant  [  [ 0, 0, .31833] ,  [ 0, 0, 1] ]
+    gm.sectionByPlane (["Vol0000"], Mgx3D.Plane(Mgx3D.Point(0, 0, .31833), Mgx3D.Vector(0, 0, 1)), "XY")
+
+    assert gm.getNbVolumes() == 2
+    assert gm.getNbSurfaces() == 9
+    out, err = capfd.readouterr()
+    assert len(err) == 0
+
 def test_3_boxes():
     ctx = Mgx3D.getStdContext()
     ctx.clearSession() # Clean the session after the previous test

Original file line number	Diff line number	Diff line change
`@@ -44,14 +44,13 @@ CommandSection::~CommandSection()`
`44`	`44`	`void CommandSection::validate()`
`45`	`45`	`{`
`46`	`46`	`for(unsigned int i=0;i<m_entities.size();i++){`
`47`		`- if(m_entities[i]->getDim()<2)`
	`47`	`+ if(m_entities[i]->getDim()<=2)`
`48`	`48`	`throw TkUtil::Exception(TkUtil::UTF8String ("Seuls des volumes géométriques peuvent être découpés", TkUtil::Charset::UTF_8));`
`49`	`49`	`}`
`50`	`50`
`51`	`51`	`int dim_tool = m_tool->getDim();`
`52`	`52`	`if(dim_tool!=2)`
`53`	`53`	`throw TkUtil::Exception(TkUtil::UTF8String ("L'outil de découpe doit être une surface", TkUtil::Charset::UTF_8));`
`54`		`-`
`55`	`54`	`}`
`56`	`55`	`/----------------------------------------------------------------------------/`
`57`	`56`	`void CommandSection::internalSpecificPreExecute()`