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.string("selection_type").setValue("topological_naming")
  67   model.do()
  68   model.end()
  69   pass
  70
  71 def getTmpFileName(ext):
  72   tempdir = tempfile.gettempdir()
  73   tmp_file = tempfile.NamedTemporaryFile(suffix=".%s"%ext, prefix='shaper_', dir=tempdir, delete=False)
  74   tmp_filename = tmp_file.name
  75   return tmp_filename
  76
  77 def testGroupsAndFieldsExportToGEOM():
  78   model.begin()
  79   partSet = model.moduleDocument()
  80   Part_1 = model.addPart(partSet)
  81   Part_1_doc = Part_1.document()
  82   Box_1 = model.addBox(Part_1_doc, 10, 10, 10)
  83   Box_2 = model.addBox(Part_1_doc, 20, 20, 20)
  84   Translation_1 = model.addTranslation(Part_1_doc, [model.selection("SOLID", "Box_1_1")], model.selection("EDGE", "PartSet/OX"), -10)
  85   Translation_1.result().setName("Translation_1_1")
  86   Partition_1 = model.addPartition(Part_1_doc, [model.selection("SOLID", "Translation_1_1"), model.selection("SOLID", "Box_2_1")])
  87   Partition_1.result().setName("Partition_1_1")
  88   Box_3 = model.addBox(Part_1_doc, 10, 10, 10)
  89   Translation_2 = model.addTranslation(Part_1_doc, [model.selection("SOLID", "Box_3_1")], model.selection("EDGE", "PartSet/OX"), 20)
  90   Translation_2.result().setName("Translation_2_1")
  91   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")]
  92   Group_1 = model.addGroup(Part_1_doc, Group_1_objects)
  93   Group_1.result().setName("faces_top")
  94   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]")]
  95   Group_2 = model.addGroup(Part_1_doc, Group_2_objects)
  96   Group_2.result().setName("edges_x")
  97   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]")]
  98   Group_3 = model.addGroup(Part_1_doc, Group_3_objects)
  99   Group_3.result().setName("vertices_bottom")
 100   Group_4 = model.addGroup(Part_1_doc, [model.selection("SOLID", "Partition_1_1_1"), model.selection("SOLID", "Translation_2_1")])
 101   Group_4.result().setName("solids_small")
 102   Field_1_objects = [model.selection("SOLID", "Partition_1_1_1"), model.selection("SOLID", "Partition_1_1_2"), model.selection("SOLID", "Translation_2_1")]
 103   Field_1 = model.addField(Part_1_doc, 1, "DOUBLE", 1, ["Comp 1"], Field_1_objects)
 104   Field_1.result().setName("Field_solids")
 105   Field_1.addStep(0, 0, [[0], [1], [2], [3]])
 106   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")]
 107   Field_2 = model.addField(Part_1_doc, 1, "DOUBLE", 1, ["Comp 1"], Field_2_objects)
 108   Field_2.result().setName("Field_faces")
 109   Field_2.addStep(0, 0, [[0], [1], [1], [1], [2], [2], [2]])
 110   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]")]
 111   Field_3 = model.addField(Part_1_doc, 1, "DOUBLE", 1, ["Comp 1"], Field_3_objects)
 112   Field_3.result().setName("Field_edges")
 113   Field_3.addStep(0, 0, [[0], [1], [1], [1], [2], [2], [2], [2], [3], [3], [3], [3]])
 114   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]")]
 115   Field_4 = model.addField(Part_1_doc, 1, "DOUBLE", 1, ["Comp 1"], Field_4_objects)
 116   Field_4.result().setName("Field_vertices")
 117   Field_4.addStep(0, 0, [[0], [1], [1], [1], [1], [1], [1], [1], [1], [2], [2], [2], [2], [3], [3], [3], [3]])
 118   model.exportToGEOM(Part_1_doc)
 119   model.do()
 120   model.end()
 121
 122 # check the groups content by the coordinates of a point on its sub-shapes
 123 def checkGroupByCoords(group, coords, tolerance=1e-7):
 124   for coord in coords:
 125     p = geompy.MakeVertex(*coord)
 126     assert geompy.MinDistance(group, p) < tolerance
 127   pass
 128
 129 ## Check the result imported in GEOM
 130 def checkResultInGEOM():
 131   geomObject_1 = getLastGEOMShape()
 132
 133   group_1_GEOM = getSubObject(geomObject_1, 1)
 134   assert group_1_GEOM.GetName() == 'faces_top'
 135   assert geompy.NumberOfFaces(group_1_GEOM) == 3
 136
 137   # coordinates of the barycenters of the faces of Group_1
 138   coords_1 = [[-5, 5, 10], [10, 10, 20], [25, 5, 10]]
 139   checkGroupByCoords(group_1_GEOM, coords_1)
 140
 141   group_2_GEOM = getSubObject(geomObject_1, 2)
 142   assert group_2_GEOM.GetName() == 'edges_x'
 143   assert geompy.NumberOfEdges(group_2_GEOM) == 3
 144
 145   # coordinates of the barycenters of the edges of Group_2
 146   coords_2 = [[-5, 0, 0], [10, 0, 0], [25, 0, 0]]
 147   checkGroupByCoords(group_2_GEOM, coords_2)
 148
 149   group_3_GEOM = getSubObject(geomObject_1, 3)
 150   assert group_3_GEOM.GetName() == 'vertices_bottom'
 151   assert geompy.NumberOfSubShapes(group_3_GEOM, geompy.ShapeType["VERTEX"]) == 8
 152
 153   # coordinates of the points of of Group_3
 154   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]]
 155   checkGroupByCoords(group_3_GEOM, coords_3)
 156
 157   group_4_GEOM = getSubObject(geomObject_1, 4)
 158   assert group_4_GEOM.GetName() == 'solids_small'
 159   assert geompy.NumberOfSolids(group_4_GEOM) == 2
 160
 161   # coordinates of the barycenters of the solids of Group_4
 162   coords_4 = [[-5, 5, 5], [25, 5, 5]]
 163   checkGroupByCoords(group_4_GEOM, coords_4)
 164
 165   field_1_GEOM = getSubObject(geomObject_1, 5)
 166   assert field_1_GEOM.GetName() == 'Field_solids'
 167   assert field_1_GEOM.GetStep(0).GetValues() == [1.0, 2.0, 3.0]
 168
 169   field_2_GEOM = getSubObject(geomObject_1, 6)
 170   assert field_2_GEOM.GetName() == 'Field_faces'
 171   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]
 172
 173   field_3_GEOM = getSubObject(geomObject_1, 7)
 174   assert field_3_GEOM.GetName() == 'Field_edges'
 175   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]
 176
 177   field_4_GEOM = getSubObject(geomObject_1, 8)
 178   assert field_4_GEOM.GetName() == 'Field_vertices'
 179   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]
 180
 181   pass
 182
 183 def checkDump():
 184   # Dump the salome study (only CORBA modules, SHAPER dump is not in it)
 185   tempdir = tempfile.gettempdir()
 186   dumpFileGeomBase = "dump_test_geom"
 187   dumpFileGeom = os.path.join(tempdir, "%s.py"%dumpFileGeomBase)
 188   salome.myStudy.DumpStudy(tempdir, dumpFileGeomBase, True, False)
 189
 190   # Dump SHAPER
 191   dumpFileShaper = os.path.join(tempdir, "dump_test_shaper.py")
 192   dumpShaper(dumpFileShaper)
 193
 194   # Load SHAPER dump
 195   exec(compile(open(dumpFileShaper).read(), dumpFileShaper, 'exec'))
 196
 197   # Load GEOM dump
 198   exec(compile(open(dumpFileGeom).read(), dumpFileGeom, 'exec'))
 199
 200   # Clean files
 201   files = [dumpFileGeom, dumpFileShaper]
 202   for f in files:
 203     os.remove(f)
 204
 205   pass
 206
 207 if __name__ == '__main__':
 208   # create 3 boxes with groups and fields
 209   testGroupsAndFieldsExportToGEOM()
 210   # check the result in GEOM
 211   checkResultInGEOM()
 212   # check that dump produces no error and can be reloaded
 213   checkDump()
 214   # check the result of the dump
 215   checkResultInGEOM()