src/ConnectorAPI/Test/TestExportToGEOMAllGroupsAndFields.py

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