src/FeaturesPlugin/Test/TestSewing_NonManifold.py

   1 # Copyright (C) 2014-2022  CEA/DEN, EDF R&D
   2 #
   3 # This library is free software; you can redistribute it and/or
   4 # modify it under the terms of the GNU Lesser General Public
   5 # License as published by the Free Software Foundation; either
   6 # version 2.1 of the License, or (at your option) any later version.
   7 #
   8 # This library is distributed in the hope that it will be useful,
   9 # but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 # Lesser General Public License for more details.
  12 #
  13 # You should have received a copy of the GNU Lesser General Public
  14 # License along with this library; if not, write to the Free Software
  15 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 #
  17 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 #
  19
  20 from salome.shaper import model
  21 from GeomAPI import GeomAPI_Shape
  22
  23 aShapeTypes = {
  24   GeomAPI_Shape.SOLID:  "GeomAPI_Shape.SOLID",
  25   GeomAPI_Shape.FACE:   "GeomAPI_Shape.FACE",
  26   GeomAPI_Shape.EDGE:   "GeomAPI_Shape.EDGE",
  27   GeomAPI_Shape.VERTEX: "GeomAPI_Shape.VERTEX"}
  28
  29 def testNbUniqueSubShapes(theFeature, theShapeType, theExpectedNbSubShapes):
  30   """ Tests number of unique feature sub-shapes of passed type for each result.
  31   :param theFeature: feature to test.
  32   :param theShapeType: shape type of sub-shapes to test.
  33   :param theExpectedNbSubShapes: list of sub-shapes numbers. Size of list should be equal to len(theFeature.results()).
  34   """
  35   aResults = theFeature.feature().results()
  36   aNbResults = len(aResults)
  37   aListSize = len(theExpectedNbSubShapes)
  38   assert (aNbResults == aListSize), "Number of results: {} not equal to list size: {}.".format(aNbResults, aListSize)
  39   for anIndex in range(0, aNbResults):
  40     aNbResultSubShapes = 0
  41     anExpectedNbSubShapes = theExpectedNbSubShapes[anIndex]
  42     aNbResultSubShapes = aResults[anIndex].shape().subShapes(theShapeType, True).size()
  43     assert (aNbResultSubShapes == anExpectedNbSubShapes), "Number of sub-shapes of type {} for result[{}]: {}. Expected: {}.".format(aShapeTypes[theShapeType], anIndex, aNbResultSubShapes, anExpectedNbSubShapes)
  44
  45 def testResults(theFeature,theModel,NbRes,NbSubRes,NbShell,NbFace,NbEdge,NbVertex):
  46   """ Tests numbers of unique sub-shapes in the results
  47   """
  48   aResults = theFeature.feature().results()
  49   aNbResults = len(aResults)
  50   assert (aNbResults == NbRes), "Number of results: {} not equal to {}}.".format(aNbResults, NbRes)
  51   theModel.testNbSubResults(theFeature, NbSubRes)
  52   testNbUniqueSubShapes(theFeature, GeomAPI_Shape.SHELL, NbShell)
  53   testNbUniqueSubShapes(theFeature, GeomAPI_Shape.FACE, NbFace)
  54   testNbUniqueSubShapes(theFeature, GeomAPI_Shape.EDGE, NbEdge)
  55   testNbUniqueSubShapes(theFeature, GeomAPI_Shape.VERTEX, NbVertex)
  56
  57 # Create document
  58 model.begin()
  59 partSet = model.moduleDocument()
  60 Part_1 = model.addPart(partSet)
  61 Part_1_doc = Part_1.document()
  62 Param_dx = model.addParameter(Part_1_doc, "dx", '10')
  63 Param_alfa = model.addParameter(Part_1_doc, "alfa", '90')
  64
  65 # =============================================================================
  66 # Test 1. Sew two shells with a single common edge
  67 # =============================================================================
  68 Box_1 = model.addBox(Part_1_doc, 10, 10, 10)
  69 Box_2 = model.addBox(Part_1_doc, 10, 10, 10)
  70 Vertex_1 = model.selection("VERTEX", "[Box_1_1/Back][Box_1_1/Left][Box_1_1/Bottom]")
  71 Vertex_2 = model.selection("VERTEX", "[Box_1_1/Front][Box_1_1/Left][Box_1_1/Top]")
  72 Translation_1 = model.addTranslation(Part_1_doc, [model.selection("SOLID", "Box_2_1")], startPoint = Vertex_1, endPoint = Vertex_2, keepSubResults = True)
  73
  74 Shell_1_objects = [model.selection("FACE", "Box_1_1/Top"), model.selection("FACE", "Box_1_1/Front")]
  75 Shell_1 = model.addShell(Part_1_doc, Shell_1_objects)
  76
  77 Shell_2_objects = [model.selection("FACE", "Translation_1_1/MF:Translated&Box_2_1/Back"), model.selection("FACE", "Translation_1_1/MF:Translated&Box_2_1/Bottom")]
  78 Shell_2 = model.addShell(Part_1_doc, Shell_2_objects)
  79 model.do()
  80
  81 Sewing_1 = model.addSewing(Part_1_doc, [model.selection("SHELL", "Shell_2_1"), model.selection("SHELL", "Shell_1_1")], 1e-07, allowNonManifold = False, alwaysCreateResult = True)
  82 model.end()
  83
  84 # sewing failed (result is a compound with the two input shells)
  85 testResults(Sewing_1, model, 1, [2], [2], [4], [14], [12])
  86
  87 # =============================================================================
  88 # Test 2. Sew the same two shells allowing non-manifold results
  89 # =============================================================================
  90 model.undo()
  91 model.begin()
  92
  93 Sewing_2 = model.addSewing(Part_1_doc, [model.selection("SHELL", "Shell_2_1"), model.selection("SHELL", "Shell_1_1")], 1e-07, allowNonManifold = True, alwaysCreateResult = True)
  94 model.end()
  95
  96 # sewing successful
  97 testResults(Sewing_2, model, 1, [0], [1], [4], [13], [10])