src/ConnectorAPI/Test/TestExportToGEOMAllGroupsAndFields.py

   1 # Copyright (C) 2014-2019  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 # -*- coding: utf-8 -*-
  20
  21 import os
  22 import tempfile
  23
  24 import sys
  25 import salome
  26
  27 salome.salome_init()
  28 theStudy = salome.myStudy
  29
  30 from salome.shaper import model
  31
  32 from salome.geom import geomBuilder
  33 geompy = geomBuilder.New()
  34
  35 ## Get the last object published in the GEOM section of the object browser
  36 def getLastGEOMShape():
  37   sb = salome.myStudy.NewBuilder()
  38   comp = salome.myStudy.FindComponent("GEOM")
  39   obj = None
  40   if comp:
  41     iterator = salome.myStudy.NewChildIterator( comp )
  42     sobj = None
  43     while iterator.More():
  44       sobj = iterator.Value()
  45       iterator.Next()
  46     if sobj:
  47       obj = sobj.GetObject()
  48   else:
  49     raise Exception("GEOM component not found.")
  50   return obj
  51
  52 ## Get the sub-object i of an object in the object browser
  53 # Numerotation starts at 1
  54 def getSubObject(obj, i):
  55   ok, sub_sobj = salome.ObjectToSObject(obj).FindSubObject(i)
  56   if not ok:
  57     raise Exception("No child found at %i for %s"%(i, obj.GetName()))
  58   sub_obj = sub_sobj.GetObject()
  59   return sub_obj
  60
  61 def dumpShaper(fileName):
  62   model.begin()
  63   dump=model.moduleDocument().addFeature("Dump")
  64   dump.string("file_path").setValue(fileName)
  65   dump.string("file_format").setValue("py")
  66   dump.boolean("topological_naming").setValue(True)
  67   dump.boolean("geometric_selection").setValue(False)
  68   dump.boolean("weak_naming").setValue(False)
  69   model.do()
  70   model.end()
  71   pass
  72
  73 def getTmpFileName(ext):
  74   tempdir = tempfile.gettempdir()
  75   tmp_file = tempfile.NamedTemporaryFile(suffix=".%s"%ext, prefix='shaper_', dir=tempdir, delete=False)
  76   tmp_filename = tmp_file.name
  77   return tmp_filename
  78
  79 def testGroupsAndFieldsExportToGEOM():
  80   model.begin()
  81   partSet = model.moduleDocument()
  82   Part_1 = model.addPart(partSet)
  83   Part_1_doc = Part_1.document()
  84   Box_1 = model.addBox(Part_1_doc, 10, 10, 10)
  85   Box_2 = model.addBox(Part_1_doc, 20, 20, 20)
  86   Translation_1 = model.addTranslation(Part_1_doc, [model.selection("SOLID", "Box_1_1")], model.selection("EDGE", "PartSet/OX"), -10)
  87   Translation_1.result().setName("Translation_1_1")
  88   Partition_1 = model.addPartition(Part_1_doc, [model.selection("SOLID", "Translation_1_1"), model.selection("SOLID", "Box_2_1")])
  89   Partition_1.result().setName("Partition_1_1")
  90   Box_3 = model.addBox(Part_1_doc, 10, 10, 10)
  91   Translation_2 = model.addTranslation(Part_1_doc, [model.selection("SOLID", "Box_3_1")], model.selection("EDGE", "PartSet/OX"), 20)
  92   Translation_2.result().setName("Translation_2_1")
  93   Group_1_objects = [model.selection("FACE", "Translation_1_1/MF:Translated&Box_1_1/Top"), model.selection("FACE", "Box_2_1/Top"), model.selection("FACE", "Translation_2_1/MF:Translated&Box_3_1/Top")]
  94   Group_1 = model.addGroup(Part_1_doc, Group_1_objects)
  95   Group_1.result().setName("faces_top")
  96   Group_2_objects = [model.selection("EDGE", "[Partition_1_1_1/Modified_Face&Box_1_1/Left][Partition_1_1_1/Modified_Face&Box_1_1/Bottom]"), model.selection("EDGE", "[Partition_1_1_2/Modified_Face&Box_2_1/Left][Partition_1_1_2/Modified_Face&Box_2_1/Bottom]"), model.selection("EDGE", "[Translation_2_1/MF:Translated&Box_3_1/Left][Translation_2_1/MF:Translated&Box_3_1/Bottom]")]
  97   Group_2 = model.addGroup(Part_1_doc, Group_2_objects)
  98   Group_2.result().setName("edges_x")
  99   Group_3_objects = [model.selection("VERTEX", "[Translation_2_1/MF:Translated&Box_3_1/Front][Translation_2_1/MF:Translated&Box_3_1/Left][Translation_2_1/MF:Translated&Box_3_1/Bottom]"), model.selection("VERTEX", "[Translation_2_1/MF:Translated&Box_3_1/Front][Translation_2_1/MF:Translated&Box_3_1/Right][Translation_2_1/MF:Translated&Box_3_1/Bottom]"), model.selection("VERTEX", "[Translation_2_1/MF:Translated&Box_3_1/Back][Translation_2_1/MF:Translated&Box_3_1/Right][Translation_2_1/MF:Translated&Box_3_1/Bottom]"), model.selection("VERTEX", "[Box_2_1/Front][Box_2_1/Right][Partition_1_1_2/Modified_Face&Box_2_1/Bottom]"), model.selection("VERTEX", "[(Partition_1_1_2/Modified_Face&Box_2_1/Back&Box_1_1/Front)(Partition_1_1_2/Modified_Face&Box_2_1/Left)(Box_2_1/Right)(Partition_1_1_2/Modified_Face&Box_2_1/Bottom)(Box_2_1/Top)][Box_2_1/Right][Partition_1_1_2/Modified_Face&Box_2_1/Bottom]"), model.selection("VERTEX", "[Partition_1_1_1/Modified_Face&Box_2_1/Back&Box_1_1/Front][Translation_1_1/MF:Translated&Box_1_1/Right][Partition_1_1_1/Modified_Face&Box_1_1/Bottom]"), model.selection("VERTEX", "[Translation_1_1/MF:Translated&Box_1_1/Back][Translation_1_1/MF:Translated&Box_1_1/Right][Partition_1_1_1/Modified_Face&Box_1_1/Bottom]"), model.selection("VERTEX", "[Translation_1_1/MF:Translated&Box_1_1/Back][Partition_1_1_1/Modified_Face&Box_1_1/Left][Partition_1_1_1/Modified_Face&Box_1_1/Bottom]")]
 100   Group_3 = model.addGroup(Part_1_doc, Group_3_objects)
 101   Group_3.result().setName("vertices_bottom")
 102   Group_4 = model.addGroup(Part_1_doc, [model.selection("SOLID", "Partition_1_1_1"), model.selection("SOLID", "Translation_2_1")])
 103   Group_4.result().setName("solids_small")
 104   Field_1_objects = [model.selection("SOLID", "Partition_1_1_1"), model.selection("SOLID", "Partition_1_1_2"), model.selection("SOLID", "Translation_2_1")]
 105   Field_1 = model.addField(Part_1_doc, 1, "DOUBLE", 1, ["Comp 1"], Field_1_objects)
 106   Field_1.result().setName("Field_solids")
 107   Field_1.addStep(0, 0, [[0], [1], [2], [3]])
 108   Field_2_objects = [model.selection("FACE", "Partition_1_1_1/Modified_Face&Box_1_1/Bottom"), model.selection("FACE", "Partition_1_1_2/Modified_Face&Box_2_1/Bottom"), model.selection("FACE", "Translation_2_1/MF:Translated&Box_3_1/Bottom"), model.selection("FACE", "Translation_1_1/MF:Translated&Box_1_1/Top"), model.selection("FACE", "Box_2_1/Top"), model.selection("FACE", "Translation_2_1/MF:Translated&Box_3_1/Top")]
 109   Field_2 = model.addField(Part_1_doc, 1, "DOUBLE", 1, ["Comp 1"], Field_2_objects)
 110   Field_2.result().setName("Field_faces")
 111   Field_2.addStep(0, 0, [[0], [1], [1], [1], [2], [2], [2]])
 112   Field_3_objects = [model.selection("EDGE", "[Translation_2_1/MF:Translated&Box_3_1/Front][Translation_2_1/MF:Translated&Box_3_1/Left]"), model.selection("EDGE", "[Translation_2_1/MF:Translated&Box_3_1/Left][Translation_2_1/MF:Translated&Box_3_1/Bottom]"), model.selection("EDGE", "[Partition_1_1_2/Modified_Face&Box_2_1/Left][Partition_1_1_2/Modified_Face&Box_2_1/Bottom]"), model.selection("EDGE", "[Partition_1_1_1/Modified_Face&Box_1_1/Left][Partition_1_1_1/Modified_Face&Box_1_1/Bottom]"), model.selection("EDGE", "[Translation_1_1/MF:Translated&Box_1_1/Back][Partition_1_1_1/Modified_Face&Box_1_1/Left]"), model.selection("EDGE", "[Translation_2_1/MF:Translated&Box_3_1/Front][Translation_2_1/MF:Translated&Box_3_1/Right]"), model.selection("EDGE", "[Translation_1_1/MF:Translated&Box_1_1/Back][Translation_1_1/MF:Translated&Box_1_1/Right]"), model.selection("EDGE", "[Box_2_1/Front][Box_2_1/Right]"), model.selection("EDGE", "[Box_2_1/Right][Box_2_1/Top]"), model.selection("EDGE", "[(Partition_1_1_2/Modified_Face&Box_2_1/Back&Box_1_1/Front)(Partition_1_1_2/Modified_Face&Box_2_1/Left)(Box_2_1/Right)(Partition_1_1_2/Modified_Face&Box_2_1/Bottom)(Box_2_1/Top)][Box_2_1/Right]"), model.selection("EDGE", "[Box_2_1/Right][Partition_1_1_2/Modified_Face&Box_2_1/Bottom]")]
 113   Field_3 = model.addField(Part_1_doc, 1, "DOUBLE", 1, ["Comp 1"], Field_3_objects)
 114   Field_3.result().setName("Field_edges")
 115   Field_3.addStep(0, 0, [[0], [1], [1], [1], [2], [2], [2], [2], [3], [3], [3], [3]])
 116   Field_4_objects = [model.selection("VERTEX", "[Translation_1_1/MF:Translated&Box_1_1/Back][Partition_1_1_1/Modified_Face&Box_1_1/Left][Partition_1_1_1/Modified_Face&Box_1_1/Bottom]"), model.selection("VERTEX", "[Translation_1_1/MF:Translated&Box_1_1/Back][Translation_1_1/MF:Translated&Box_1_1/Right][Partition_1_1_1/Modified_Face&Box_1_1/Bottom]"), model.selection("VERTEX", "[Partition_1_1_1/Modified_Face&Box_2_1/Back&Box_1_1/Front][Translation_1_1/MF:Translated&Box_1_1/Right][Partition_1_1_1/Modified_Face&Box_1_1/Bottom]"), model.selection("VERTEX", "[(Partition_1_1_2/Modified_Face&Box_2_1/Back&Box_1_1/Front)(Partition_1_1_2/Modified_Face&Box_2_1/Left)(Box_2_1/Right)(Partition_1_1_2/Modified_Face&Box_2_1/Bottom)(Box_2_1/Top)][Box_2_1/Right][Partition_1_1_2/Modified_Face&Box_2_1/Bottom]"), model.selection("VERTEX", "[Box_2_1/Front][Box_2_1/Right][Partition_1_1_2/Modified_Face&Box_2_1/Bottom]"), model.selection("VERTEX", "[Translation_2_1/MF:Translated&Box_3_1/Back][Translation_2_1/MF:Translated&Box_3_1/Right][Translation_2_1/MF:Translated&Box_3_1/Bottom]"), model.selection("VERTEX", "[Translation_2_1/MF:Translated&Box_3_1/Front][Translation_2_1/MF:Translated&Box_3_1/Right][Translation_2_1/MF:Translated&Box_3_1/Bottom]"), model.selection("VERTEX", "[Translation_2_1/MF:Translated&Box_3_1/Front][Translation_2_1/MF:Translated&Box_3_1/Left][Translation_2_1/MF:Translated&Box_3_1/Bottom]"), model.selection("VERTEX", "[Translation_1_1/MF:Translated&Box_1_1/Back][Partition_1_1_1/Modified_Face&Box_1_1/Left][Translation_1_1/MF:Translated&Box_1_1/Top]"), model.selection("VERTEX", "[Translation_1_1/MF:Translated&Box_1_1/Back][Translation_1_1/MF:Translated&Box_1_1/Right][Translation_1_1/MF:Translated&Box_1_1/Top]"), model.selection("VERTEX", "[(Partition_1_1_2/Modified_Face&Box_2_1/Back&Box_1_1/Front)(Partition_1_1_2/Modified_Face&Box_2_1/Left)(Box_2_1/Right)(Partition_1_1_2/Modified_Face&Box_2_1/Bottom)(Box_2_1/Top)][Partition_1_1_2/Modified_Face&Box_2_1/Left][Box_2_1/Top]"), model.selection("VERTEX", "[(Partition_1_1_2/Modified_Face&Box_2_1/Back&Box_1_1/Front)(Partition_1_1_2/Modified_Face&Box_2_1/Left)(Box_2_1/Right)(Partition_1_1_2/Modified_Face&Box_2_1/Bottom)(Box_2_1/Top)][Box_2_1/Right][Box_2_1/Top]"), model.selection("VERTEX", "[Box_2_1/Front][Partition_1_1_2/Modified_Face&Box_2_1/Left][Box_2_1/Top]"), model.selection("VERTEX", "[Box_2_1/Front][Box_2_1/Right][Box_2_1/Top]"), model.selection("VERTEX", "[Translation_2_1/MF:Translated&Box_3_1/Front][Translation_2_1/MF:Translated&Box_3_1/Left][Translation_2_1/MF:Translated&Box_3_1/Top]"), model.selection("VERTEX", "[Translation_2_1/MF:Translated&Box_3_1/Front][Translation_2_1/MF:Translated&Box_3_1/Right][Translation_2_1/MF:Translated&Box_3_1/Top]")]
 117   Field_4 = model.addField(Part_1_doc, 1, "DOUBLE", 1, ["Comp 1"], Field_4_objects)
 118   Field_4.result().setName("Field_vertices")
 119   Field_4.addStep(0, 0, [[0], [1], [1], [1], [1], [1], [1], [1], [1], [2], [2], [2], [2], [3], [3], [3], [3]])
 120   model.exportToGEOM(Part_1_doc)
 121   model.do()
 122   model.end()
 123
 124 # check the groups content by the coordinates of a point on its sub-shapes
 125 def checkGroupByCoords(group, coords, tolerance=1e-7):
 126   for coord in coords:
 127     p = geompy.MakeVertex(*coord)
 128     assert geompy.MinDistance(group, p) < tolerance
 129   pass
 130
 131 ## Check the result imported in GEOM
 132 def checkResultInGEOM():
 133   geomObject_1 = getLastGEOMShape()
 134
 135   group_1_GEOM = getSubObject(geomObject_1, 1)
 136   assert group_1_GEOM.GetName() == 'faces_top'
 137   assert geompy.NumberOfFaces(group_1_GEOM) == 3
 138
 139   # coordinates of the barycenters of the faces of Group_1
 140   coords_1 = [[-5, 5, 10], [10, 10, 20], [25, 5, 10]]
 141   checkGroupByCoords(group_1_GEOM, coords_1)
 142
 143   group_2_GEOM = getSubObject(geomObject_1, 2)
 144   assert group_2_GEOM.GetName() == 'edges_x'
 145   assert geompy.NumberOfEdges(group_2_GEOM) == 3
 146
 147   # coordinates of the barycenters of the edges of Group_2
 148   coords_2 = [[-5, 0, 0], [10, 0, 0], [25, 0, 0]]
 149   checkGroupByCoords(group_2_GEOM, coords_2)
 150
 151   group_3_GEOM = getSubObject(geomObject_1, 3)
 152   assert group_3_GEOM.GetName() == 'vertices_bottom'
 153   assert geompy.NumberOfSubShapes(group_3_GEOM, geompy.ShapeType["VERTEX"]) == 8
 154
 155   # coordinates of the points of of Group_3
 156   coords_3 = [[-10, 0, 0], [-10, 10, 0], [0, 10, 0], [0, 20, 0], [20, 20, 0], [20, 10, 0], [30, 10, 0], [30, 0, 0]]
 157   checkGroupByCoords(group_3_GEOM, coords_3)
 158
 159   group_4_GEOM = getSubObject(geomObject_1, 4)
 160   assert group_4_GEOM.GetName() == 'solids_small'
 161   assert geompy.NumberOfSolids(group_4_GEOM) == 2
 162
 163   # coordinates of the barycenters of the solids of Group_4
 164   coords_4 = [[-5, 5, 5], [25, 5, 5]]
 165   checkGroupByCoords(group_4_GEOM, coords_4)
 166
 167   field_1_GEOM = getSubObject(geomObject_1, 5)
 168   assert field_1_GEOM.GetName() == 'Field_solids'
 169   assert field_1_GEOM.GetStep(0).GetValues() == [1.0, 2.0, 3.0]
 170
 171   field_2_GEOM = getSubObject(geomObject_1, 6)
 172   assert field_2_GEOM.GetName() == 'Field_faces'
 173   assert field_2_GEOM.GetStep(0).GetValues() == [0.0, 0.0, 0.0, 0.0, 1.0, 2.0, 0.0, 0.0, 0.0, 0.0, 1.0, 2.0, 0.0, 0.0, 0.0, 0.0, 1.0, 2.0]
 174
 175   field_3_GEOM = getSubObject(geomObject_1, 7)
 176   assert field_3_GEOM.GetName() == 'Field_edges'
 177   assert field_3_GEOM.GetStep(0).GetValues() == [2.0, 0.0, 2.0, 0.0, 0.0, 0.0, 0.0, 0.0, 2.0, 0.0, 0.0, 0.0, 0.0, 0.0, 3.0, 0.0, 0.0, 0.0, 3.0, 0.0, 1.0, 0.0, 3.0, 3.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 2.0, 0.0, 1.0, 0.0, 0.0, 0.0]
 178
 179   field_4_GEOM = getSubObject(geomObject_1, 8)
 180   assert field_4_GEOM.GetName() == 'Field_vertices'
 181   assert field_4_GEOM.GetStep(0).GetValues() == [2.0, 1.0, 2.0, 1.0, 0.0, 0.0, 0.0, 1.0, 2.0, 2.0, 1.0, 3.0, 0.0, 3.0, 1.0, 0.0, 0.0, 0.0, 1.0, 3.0, 1.0, 3.0, 1.0]
 182
 183   pass
 184
 185 def checkDump():
 186   # Dump the salome study (only CORBA modules, SHAPER dump is not in it)
 187   tempdir = tempfile.gettempdir()
 188   dumpFileGeomBase = "dump_test_geom"
 189   dumpFileGeom = os.path.join(tempdir, "%s.py"%dumpFileGeomBase)
 190   salome.myStudy.DumpStudy(tempdir, dumpFileGeomBase, True, False)
 191
 192   # Dump SHAPER
 193   dumpFileShaper = os.path.join(tempdir, "dump_test_shaper.py")
 194   dumpShaper(dumpFileShaper)
 195
 196   # Load SHAPER dump
 197   exec(compile(open(dumpFileShaper).read(), dumpFileShaper, 'exec'))
 198
 199   # Load GEOM dump
 200   exec(compile(open(dumpFileGeom).read(), dumpFileGeom, 'exec'))
 201
 202   # Clean files
 203   files = [dumpFileGeom, dumpFileShaper]
 204   for f in files:
 205     os.remove(f)
 206
 207   pass
 208
 209 if __name__ == '__main__':
 210   # create 3 boxes with groups and fields
 211   testGroupsAndFieldsExportToGEOM()
 212   # check the result in GEOM
 213   checkResultInGEOM()
 214   # check that dump produces no error and can be reloaded
 215   checkDump()
 216   # check the result of the dump
 217   checkResultInGEOM()